FLIGHT  CONTROL  SYSTEM 
RECONFIGURATION  DESIGN  USING 
QUANTITATIVE  FEEDBACK  THEORY 

THESIS 


AFIT/GE/ENG/84D-15  Phillip  B.  Arnold 

Major  •  USAF 


DEPARTMENT  OF  THE  AIR  FORCE 
AIR  UNIVERSITY 

AIR  FORCE  INSTITUTE  OF  TECHNOLOGY 


tEGT  rfiymmanl  limf  gpprmtT 
for  public  ratoaw  tad  Mbf  fei 
dlitrtbtitloo  b  i  liiitoi.  <  0«4||| 


DT1C 

iELECTE 


SAPR01B86 

X  E 


Wright-Patterson  Air  Force  Base,  Ohio 


AFIT/GE/ENG/84D-15 


FLIGHT  CONTROL  SYSTEM 
RECONFIGURATION  DESIGN  USING 
QUANTITATIVE  FEEDBACK  THEORY 

THESIS 


AFIT/GE/ENG/84D-15  Phillip  B.  Arnold 

Major  •  USAF 


Approved  for  public  release;  distribution  unlimited 


DTIC 

ELETCTE 


APR  0  1 1985 


AFIT/GE/ENG/84D-15 


FLIGHT  CONTROL  SYSTEM  RECONFIGURATION  DESIGN 
USING  QUANTITATIVE  FEEDBACK  THEORY 

THESIS 

Presented  to  the  Faculty  of  the  School  of  Engineering 
of  the  Air  Force  Institute  of  Technology 
Air  University 

in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of 
Master  of  Science 


Phillip  B.  Arnold,  B.S.E.E. 

Major  USAF 

Graduate  Electrical  Engineering 
December  1984 

Approved  for  public  release;  distribution  limited 


Preface 


I  extend  my  sincere  thanks  to  my  thesis  advisors 
Professors  Constantine  H.  Houpis,  Isaac  Horowitz,  and 
John  J.  D'Azzo  for  their  cooperation  and  patience  through¬ 
out  this  thesis  effort.  Their  help  contributed  immensely 
to  the  successful  completion  of  this  thesis.  A  special 
acknowledgement  is  given  to  Mr.  Yin-Kuei  Liao  of  the 
Chung-Shan  Institute  of  Science,  Republic  of  China.  His 
expert  assistance  given  so  freely  and  tirelessly  is  greatly 
appreciated. 

My  fellow  classmates  in  the  flights  controls 
sequence  provided  many  new  insights  and  fresh  ideas  during 
countless  informal  discussions.  For  their  contributions, 

I  extend  my  sincere  gratitude  to  Major  Terry  Courtheyn, 
Captain  Harvey  Russell  of  the  Canadian  Armed  Forces,  and 
Lieutenant  Robert  Eslinger. 

My  family  has  been  very  understanding  throughout 
this  demanding  ordeal.  Their  encouragement  and  enthusiasm 
sustained  me  during  the  long  days  and  nights.  To  my  wife 
Marybeth,  and  my  daughters  Carey  and  Becky,  I  give  my  love. 


Table  of  Contents 


Page 

Preface .  ii 

List  of  Figures .  v 

List  of  Tables . xiii 

Abstract . xix 

I.  Introduction  .  1 

1 . 1  Background .  2 

1.2  Problem .  3 

1.3  Assumptions  .  4 

1 . 4  Approach .  6 

1.5  Presentation  .  7 

II.  Aircraft .  8 


II.  1  Introduction .  8 

11. 2  Aerodynamic  Model  .  8 

11. 3  Summary .  15 

III.  Reconfiguration  Theory  .  16 

•  III.l  Introduction  .  16 

111. 2  Reconfiguration  Theory  and  System 

Model . 16 

111. 3  Equivalent  Plant  Development  .  19 

III.  4  Summary .  22 

IV.  Compensator  Designs  . .  23 


IV. 1  Introduction  .  , 

IV. 2  Plant  Equations  . 

IV. 3  Specifications  . 

IV. 4  Loop  Bound  Design,  L2q  . 
IV. 5  Loop  Bound  Design,  L. 

IV. 6  Stability  Design  for  Three 
Simultaneous  Failures  .  . 
IV.  7  Summary . 


23 

23 

25 

29 

40 

45 

51 


4  / 


« 


Page 


V.  Results . .  53 

V.l  Introduction .  53 

V.2  Simulation .  53 

V.2.1  Pitch  Rate  Command  Response  .  .  54 

V.2. 2  Roll  Rate  Command  Response  ...  67 

V.  3  Saturation .  80 

V.4  QFT  Insights .  86 

*  V.5  Control  Surface  "Hardover"  .  87 

V.6  Summary .  92 


VI.  Conclusions  and  Recommendations 


94 


VI.  1  Discussion .  94 

VI.  2  Conclusions .  95 

VI .  3  Recommendations .  96 

Appendix  As  Flight  Parameters  and 

Aerodynamic  Derivatives  .  97 

Appendix  B:  Reconfiguration  Theory  .  102 

Appendix  C:  Equivalent  Plant  Transfer 

Function  Matrices  .  131 

Appendix  D:  Pitch  Rate  and  Roll  Rate 

Response  Models  .  152 

Appendix  E:  Templates  of  and  Q22 . 156 

Appendix  F:  Simulation  Programs  .  163 

Appendix  G:  Control  Surface  Deflections  .  189 

Bibliography . ■ . 208 

Vita . 209 


Accession  Fc 

NTIS  GRA&I 
DTIC  TAB 
Unonnoujieed 
Justif icatic 

)r 

□ 

in  _ 

By 

Distribution/ 
Availability  Codes 

Dist 

Avail 

Spec 

and/or 

ial 

"A 


Vv’,1 


■ .  -  >  a  v*  y.vvv.v, 

liaiaiaMimin Mtiatoi 


iv 


Figure 

1-1 


Flight  Envelope  Design  Points 


3- 1 

4- 1 
4-2 

4-3 

4-4 

4-5 

4-6 

4-7 

4-8 

4-9 

4-10 

4-11 

4-12 

4- 13 

5- 1 

5-2 (a) 
5-2 (b) 


Plant  and  Compensation  Model  Signal 

Flow  Diagram  . 

Pitch  Rate  Tolerances  in  the  Time  Domain  .  . 

Pitch  Rate  Tolerances  in  the  Frequency 
Domain  . 

Roll  Rate  Tolerances  in  the  Time  Domain  .  . 

Roll  Rate  Tolerances  in  the  Frequency 
Domain  . 

Equivalent  S1SO  Systems  . 

Template  of  Qi;l  at  u>  =  0.1  Radians /Second 

Template  of  Q22  at  w  =  0.1  Radians/Second 

Nichol's  Chart  Plot  of  L2q  and  Bounds  .  .  . 

Bode  Plot  of  L„  . 

2o 

Nichol's  Chart  Plot  of  L1q,  Llo‘* ,  and  Bounds 

Bode  Plot  of  . . 

L2q  and  Bounds  for  Stability  Only 

(On!  Elevator  Only) . 

L^q  and  Bounds  for  Stability  Only 

(On!  Elevator  Only) . 

Bode  Plot  Comparison  of  and  the 
Approximation  for  . . 

Pitch  Rate  Responses,  Flight  Condition  One  . 

Pitch  Rate  Responses,  Flight  Condition  Two  . 


Page 


Figure 


5-2 (c)  Pitch  Rate  Responses,  Flight  Condition 

Three . 57 

5-2 (d)  Pitch  Rate  Responses,  Flight  Condition 

Four  . 57 

5-3 (a)  Roll  Rate  Cross-Coupling  Responses, 

.  Flight  Condition  One  .  61 

5-3 (b)  Roll  Rate  Cross-Coupling  Responses, 

Flight  Condition  Two . 61 

5-3 (c)  Roll  Rate  Cross-Coupling  Responses, 

Flight  Condition  Three  .  62 

5-3  (d)  Roll  Rate  Cross-Coupling  Responses, 

Flight  Condition  Four . 62 

5-4 (a)  Elevator  and  Flaperon  Deflections  for 
Pitch  Rate,  CSC  Mode  One  and  Flight 
Condition  One . 64 

5-4 (b)  Elevator  and  Flaperon  Deflections  for 
Pitch  Rate,  CSC  Mode  Two  and  Flight 
Condition  One . 64 

5-4 (c)  Elevator  and  Flaperon  Deflections  for 
Pitch  Rate,  CSC  Mode  Three  and  Flight 
Condition  One . 65 

5-4j(d)  Elevator  and  Flaperon  Deflections  for 
Pitch  Rate,  CSC  Mode  Four  and  Flight 
Condition  One . 65 

5-4 (e)  Elevator  and  Flaperon  Deflections  for 
Pitch  Rate,  CSC  Mode  Five  and  Flight 
Condition  One . 66 

5-4 (f)  Elevator  and  Flaperon  Deflections  for 
Pitch  Rate,  CSC  Mode  Six  and  Flight 
Condition  One . 66 


5-5 (a)  Roll  Rate  Responses,  Flight  Condition  One  .  68 
5-5 (b)  Roll  Rate  Responses,  Flight  Condition  Two  .  68 
5-5 (c)  Roll  Rate  Responses,  Flight  Condition  Three.  69 
5-5 (d)  Roll  Rate  Responses,  Flight  Condition  Four  .  69 


Figure  Page 

5-6  Bode  Plot  Comparison  of  G-  and  the 

Approximation  for  G ^ . 72 

5-7  (a)  Pitch  Rate  Cross-Coupling  Responses, 

Flight  Condition  One . 74 

5-7 (b)  Pitch  Rate  Cross-Coupling  Responses, 

Flight  Condition  Two . 74 

5-7  (c)  Pitch  Rate  Cross-Coupling  Responses, 

Flight  Condition  Three  .  75 

5-7 (d)  Pitch  Rate  Cross-Coupling  Responses, 

Flight  Condition  Four . 75 

5-8  (a)  Elevator  and  Flaperon  Deflections  for 
Roll  Rate,  CSC  Mode  One  and  Flight 
Condition  One . 77 

5-8  (b)  Elevator  and  Flaperon  Deflections  for 
Roll  Rate,  CSC  Mode  Two  and  Flight 
Condition  One  .....  .  77 

5-8  (c)  Elevator  and  Flaperon  Deflections  for 
Roll  Rate,  CSC  Mode  Three  and  Flight 
Condition  One  .....  .  78 

5-8 (d)  Elevator  and  Flaperon  Deflections  for 
Roll  Rate,  CSC  Mode  Four  and  Flight 
Condition  One . 78 

5-8  (e)  Elevator  and  Flaperon  Deflections  for 
Roll  Rate,  CSC  Mode  Five  and  Flight 
Condition  One . 79 

5-8  (f)  Elevator  and  Flaperon  Deflections  for 
Roll  Rate,  CSC  Mode  Six  and  Flight 
Condition  One  . . 79 

5-9  Elevator  and  Flaperon  Deflections  With 

Saturation,  CSC  Mode  One  and  Flight 
Condition  One . 81 

5-10  Roll  Rate  Response  With  and  Without 

Saturation,  CSC  Mode  One  and  Flight 
Condition  One . 81 

5-11  Roll  Rate  Response  With  and  Without 
Saturation,  CSC  Mode  Two  and  Flight 
Condition  One . 82 


vii 


Page 


Figure 


5-12  (a) 

Pitch  Rate  Cross-coupling  Response 

With  and  Without  Saturation,  CSC  Mode 
Two  and  Flight  Condition  One  .... 

82 

•  •  • 

5-12 (b) 

Elevator  and  Flaperon  Deflections  for 
CSC  Mode  Two  and  Flight  Condition  One 

83 

•  •  • 

_  J 

5-13 

Elevator  and  Flaperon  Deflections  for 
Roll  Rate  with  2  =  0.5,  CSC  Mode  One 

and  Flight  Condition  One 

...  84 

5-14 

Elevator  and  Flaperon  Deflections  for 
Pitch  Rate  with  a  Ramped  Input,  0.5 
sec  to  Peak  Value  (CSC  Mode  One  and 
Flight  Condition  One)  . 

.  .  .  85 

5-15 

Elevator  and  Flaperon  Deflections  for 
Roll  Rate  with  a  Ramped  Input,  1  sec 
to  Peak  Value  (CSC  Mode  One  and  Flight 
Condition  One)  . 

.  .  .  85 

5-16 

Template  of  at  u  =  1.0  r/s. 

Includes  =0  . 

.  .  .  88 

5-17 

Bounds  on  L  Resulting  from  t 11  Demands 

.  .  .  89 

“  .  - 

5-18  (a) 

Pitch  Rate  Response  from  Elevator 
Hardover,  Flight  Condition  Two  .... 

.  .  .  90 

5-18 (b) 

Roll  Rate  Response  from  Elevator 
Hardover,  Flight  Condition  Two  .... 

.  .  .  91 

•  >-  -»>  J 

5-19  (a) 

Pitch  Rate  Response  from  Flaperon 
Hardover,  Flight  Condition  Two,.  .  .  . 

.  .  .  91 

V*".**" 

5-19 (b) 

Roll  Rate  Response  from  Flaperon 
Hardover,  Flight  Condition  Two  .... 

.  .  .  92 

- 

B-l 

A  4-input,  2-output,  2-mode  structure 

.  .  .  128 

B-2 

A  5-input,  3-output,  2-mode  structure 

.  .  .  129 

B-3 

A  10-input,  6-output,  2-mode  structure 

.  .  .  130 

E-l(a) 

Template  of  at  10  -  0.5  Rad/Sec  .  . 

.  .  .  157 

E-l(b) 

Template  of  Q22  at  to  =  0.5  Rad/Sec  .  . 

.  .  .  157 

-1 

viii 

■  .  ■  .  ~ 4 

Page 

158 


Figure 

E-2  (a)  Template  at  Qn  at  u  =  4.0  Rad/Sec  .... 

E-2  (b)  Template  of  Q22  at  oj  =  4.0  Rad/Sec  ....  158 

E-3  (a)  Template  of  Q.^  at  oj  =  8  Rad/Sec . 159 

E-3  (b)  Template  of  Q22  at  oj  =  8  Rad/Sec . 159 

E-4 (a)  Template  of  Q1]L  at  oj  =10  Rad/Sec . 160 

E-4(b)  Template  of  Q22  at  cj  =10  Rad/Sec . 160 

E-5(a)  Template  of  at  oj  =20  Rad/Sec .  161 


E-5  (b)  Template  of  Q22  at  oj  =20  Rad/Sec .  161 

E-6  (a)  Template  of  Q-^  at  oj  =50  Rad/Sec  .  .  .  .  .  162 

E-6  (b)  Template  of  Q22  at  oj  =50  Rad/Sec . 162 


G-l(a)  Control  Surface  Deflections  for  Pitch 
Rate  Command,  Flight  Condition  2 — 

CSC  Mode  1 .  190 

G-l  (b)  Control  Surface  Deflections  for  Pitch 
Rate  Command,  Flight  Condition  2 — 

CSC  Mode  2 .  190 

G-l  (c)  Control  Surface  Deflections  for  Pitch 
Rate  Command,  Flight  Condition  2-- 
CSC  Mode  3 .  191 

G-l  (d)  Control  Surface  Deflections  for  Pitch 
Rate  Command,  Flight  Condition  2-- 
CSC  Mode  4 . 191 

G-l (e)  Control  Surface  Deflections  for  Pitch 
Rate  Command,  Flight  Condition  2— 

CSC  Mode  5 .  192 

G-l  (f )  Control  Surface  Deflections  for  Pitch 

Rate  Command,  Flight  Condition  2 — 

CSC  Mode  6 .  192 

G-2  (a)  Control  Surface  Deflections  for  Pitch 
Rate  Command,  Flight  Condition  3 — 

CSC  Mode  1 .  193 

G-2 (b)  Control  Surface  Deflections  for  Pitch 
Rate  Command,  Flight  Condition  3-- 
CSC  Mode  2 .  193 


significantly  impacts  LCC  as  actuators  are  50  percent  of 
FCS  LCC.  Correspondingly,  simpler  actuators  reduce  main¬ 
tenance  actions  required  (4). 

Quantitative  feedback  theory  is  ideally  suited  for 
the  FCS  reconfiguration  problem.  In  QFT  the  design  is 
accomplished  a  priori  so  that  the  system  output  is  in  the 
acceptable  set  of  outputs  for  all  possible  members  of  the 
plant  set.  Surface  failures  simply  increase  the  number  of 
equations  in  the  equivalent  plant  set,  as  is  seen  in 
Chapter  III.  The  transparency  of  QFT  enables  the  designer 
to  readily  see  the  extra  "cost  of  feedback"  for  this  enlarge¬ 
ment  of  the  plant  equation  set  by  including  surface  failures 
(5) .  A  QFT  design  approach  also  reduces  the  need  for  com¬ 
plex  fault  identification  and  isolation  (FDI)  schemes. 

I . 2  Problem 

The  main  purpose  of  this  thesis  effort  is  to 
investigate  the  usefulness  of  QFT  in  designing  for  FCS 
reconfiguration.  The  minimum  phase  plant  equations  are 
based  upon  linearized  equations  of  motion  for  the  AFTI/F-16. 
Non-minimum  phase  plant  equations  are  not  covered  in  this 
thesis.  Using  QFT,  control  laws  are  developed  to  achieve 
specified  performance  tolerances  for  several  flight  con¬ 
trol  surface  failures  over  a  large  portion  of  the  flight 
envelope.  Pitch  rate  (q)  and  roll  rate  (p)  are  the  two 
output  variables  to  be  controlled.  Control  surfaces  are 
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In  this  thesis,  QFT  is  applied  to  a  FCS  reconfiguration 
problem  in  an  attempt  to  gain  insight  into  its  utility  as 
a  design  tool. 

The  remainder  of  the  first  chapter  presents  per¬ 
tinent  background  information,  a  statement  and  scope  of  the 
problem,  the  assumptions  used,  the  approach  taken,  and  the 
sequence  of  presentation  for  the  thesis. 

I . 1  Background 

Several  promising  benefits  of  inherent  reconfigura¬ 
tion  have  been  detailed  in  Self- Repairing  Flight  Control 
System  Reliability  and  Maintainability,  Program  Plan  (4)  . 
The  primary  benefits  to  be  derived  from  inherent  recon¬ 
figuration  include  improvements  in  reliability  and  maintain 
ability  (R&M) ,  survivability,  and  reduced  life  cycle  costs 
(LCC) . 

Improvements  in  R&M  are  expected  as  a  result  of 
fewer  and  less  complex  FCS  components.  Survivability  will 
be  improved  by  taking  advantage  of  existing  redundant  con¬ 
trol  surfaces,  thereby  reducing  the  dependency  on  any  one 
set  of  control  surfaces  for  primary  control  about  an  axis. 
The  reduction  in  criticality  of  the  control  surfaces  allows 
a  reduction  in  the  redundancy  level  at  the  actuator.  The 
actuator  no  longer  becomes  the  redundancy  integrator  and 
the  actuator  then  can  be  made  simple  and  less  complex. 
Simplex  or  dual  actuators  can  be  contemplated  which 
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FLIGHT  CONTROL  SYSTEM  RECONFIGURATION  DESIGN 


USING  QUANTITATIVE  FEEDBACK  THEORY 


I.  Introduction 


Design  planning  for  the  next  generation  fighter 
aircraft  is  presently  in  progress.  The  aircraft,  known 
as  the  Advanced  Tactical  Fighter  (ATF) ,  will  have  a  flight 
control  system  (FCS)  capable  of  taking  advantage  of  the 
inherent  redundancy  existing  in  redundant  flight  control 
surfaces.  This  means  that  no  single  loss  or  failure  of 
a  FCS  component  will  result  in  loss  of  aircraft  control. 

A  digital  flight  control  system  will  redistribute  control 
authority  to  the  remaining  components  to  generate  the  neces¬ 
sary  aerodynamic  forces  and  moments.  This  reconfiguration 
of  control  authority  by  an  aircraft  with  inherent  redun¬ 
dancy  gives  rise  to  the  term  "inherent  reconfiguration." 

Methods  of  developing  reconfiguration  schemes  are 
now  being  investigated  by  the  Flight  Dynamics  Laboratory 
(1) .  One  promising  approach  is  quantitative  feedback  theory 
(QFT)  as  developed  and  successfully  applied  to  a  number  of 
difficult  problems  by  Dr.  Isaac  Horowitz  (2;  3) .  Quantita¬ 
tive  feedback  theory  is  ideally  suited  to  develop  robust 
control  laws  for  systems  with  highly  uncertain  parameters. 
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Fixed-compensators  and  prefilters  are  synthesized 
for  the  pitch  rate  and  roll  rate  channels,  and  a  digital 
simulation  is  conducted  at  four  points  in  the  flight 
envelope  with  up  to  two  simultaneous  surface  failures,  as 
well  as,  for  the  unimpaired  aircraft.  Results  of  the  simu¬ 
lation  show  that  the  compensators  provide  robust  control 
for  each  control  surface  configuration  at  each  flight  condi 
tion.  Compensator  designs  are  also  synthesized  to  provide 
only  stability  when  three  of  the  four  surfaces  are  failed. 

The  theory  developed  by  Dr.  Horowitz  is  found  to 
be  effective  in  designing  for  robustness  despite  large 
uncertainty  arising  from  control  surface  failures,  and 
plant  parameter  variation.  Flight  control  research  using 
quantitative  feedback  theory  should  be  continued.  The 
fixed-compensator  design  approach  has  potential  to  signifi¬ 
cantly  reduce  the  complexity  of  a  reconf igurable  flight 
control  system.  A  computer-aided  design  package  for  quanti 
tative  feedback  theory  should  be  developed  to  facilitate 
design  work. 
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Abstract 


Quantitative  theory  developed  by  Dr.  Isaac  Horowitz 
of  the  Weizmann  Institute  of  Science,  Israel  is  used  to 
develop  control  laws  for  the  AFTI/F-16  with  a  reconfigur- 
able  flight  control  system.  Compensators  are  synthesized 
to  control  pitch  rate  and  roll  rate  through  individually 
controlled  elevators  and  flaperons.  Robust  control  of 
these  variables  is  required  over  a  large  portion  of  the 
flight  envelope  despite  flight  control  surface  failures. 

Linearized  aerodynamic  data  are  used  to  develop 
the  aircraft  model  in  state-variable  format.  The  longi¬ 
tudinal  and  lateral-directional  equations  are  coupled  in 
the  control  matrix.  Individual  control  of  the  elevators 
and  flaperons  is  obtained  by  dividing  the  dimensionalized 
control  derivatives  for  a  control  surface  pair  in  half  and 
assigning  each  surface  of  the  pair  one-half  of  the  total 
derivative  value.  The  system  with  individually  controlled 
surfaces  represents  a  four  input-two  output  system  which  is 
transformed  into  an  equivalent  two  input-two  output  system 
for  each  control  surface  configuration  and  flight  condition 
Quantitative  feedback  theory  is  then  applied  to  the  equiva¬ 
lent  systems. 


8  —  Pitch  angle 

<p  —  Roll  angle 

u. .  —  Control  surface  constants 

ID 

a)  —  Frequency 


xviii 


State  vector 


Dimensional  variation  of  x- force  with  pitch  angle 

Dimensional  variation  of  x-force  with  forward 
velocity 

Dimensional  variation  of  x-force  with  angle-of- 
attack 

Dimensional  variation  of  x-force  with  pitch  rate 

Dimensional  variation  of  x-force  with  elevator  (6..) 
and  flaperon  (S2) 

Dimensional  variation  of  y-force  with  roll  angle 

Dimensional  variation  of  y-force  with  sideslip 
angle 

Dimensional  variation  of  y-force  with  roll  rate 

Dimensional  variation  of  y-force  with  yaw  rate 

Dimensional  variation  of  y-force  with  elevator  (6..) 
or  flaperon 

Dimensional  variation  of  z-force  with  pitch  angle 

Dimensional  variation  of  z-force  with  forward 
velocity 

Dimensional  variation  of  z-force  with  angle-of- 
attack 

Dimensional  variation  of  z-force  with  pitch  rate 

Dimensional  variation  of  z-force  with  elevator  (5,) 
or  flaperon  (62) 

Elevator  input 

Flaperon  input 

Rudder  input 

Canard  input 


Peak  value 


M  —  Dimensional  variation  of  pitching  moment  with 
a  pitch  angle 

M  —  Dimensional  variation  of  pitching  moment  with 
u  forward  velocity 


M 

a 


—  Dimensional  variation  of  pitching  moment  with 
angle-of- attack 


Dimensional  variation  of  pitching  moment  with 
pitch  rate 


M.  —  Dimensional  variation  of  pitching  moment  with 
elevator  (6^  or  flaperon  (&2) 

Nft  —  Dimensional  variation  of  yawing  moment  with  side- 
**  slip  angle 


Dimensional  variation  of  yawing  moment  with 
roll  rate 


N  —  Dimensional  variation  of  yawing  moment  with 
r  yaw  rate 

N,.  —  Dimensional  variation  of  yawing  moment  with 

0  elevator  (6^)  or  flaperon  (<$2) 

p  —  Roll  rate 

Pp  —  Equivalent  plant  matrix 

P„  —  Basic  plant  equations 

q  —  Pitch  rate 

Q.  .  —  Plant  equations  for  equivalent  single  input-single 

output  systems 

r  —  Yaw  rate 


s  —  Wing  reference  area 

T  —  System  transfer  function  matrix  ( ^ ) 

t  —  Settling  time 

tr  —  Rise  time 

u  —  Forward  velocity 


xvi 


a  —  Angle-of-attack 

A  #  —  Plant  matrix 

B  —  Control  matrix 

—  Performance  specification  bounds 
b  —  Wing  span 

0  —  Sideslip  angle 

C  —  Wing  aerodynamic  chord 
Ft  —  Feet 

lbs  —  Pounds 

g^  —  Pitch  rate  compensator  transfer  function 

g2  —  Roll  rate  compensator  transfer  function 

f. .  —  Prefilter 

ID 

I  —  Moment  of  inertia  about  x-axis 
xx 

I  —  Moment  of  inertia  about  y-axis 

yy 

I  —  Moment  of  inertia  about  z-axis 
zz 

I  —  Product  of  inertia  in  xyz  system 

xz 

Lg  —  Dimensional  variation  of  rolling  moment  with 
p  sideslip  angle 

L  —  Dimensional  variation  of  rolling  moment  with 
"  roll  rate 


L  —  Dimensional  variation  of  rolling  moment  with 
yaw  rate 

—  Dimensional  variation  of  rolling  moment  with 
elevator  (6^)  or  flaperon  (<S2) 

--  Pitch  rate  nominal  loop  transmission 

L2q  —  Roll  rate  nominal  loop  transmission 


Table 


Equivalent  Plant  Equations  for  the  Single 
Input-Single  Output  Systems  for  Flight 
Condition  Two  . 


C-6 


C-7 

C-8  • 


Equivalent  Plant  Equations  for  the  Single 
Input-Single  Output  Systems  for  Flight 
Condition  Three  . 

Equivalent  Plant  Equations  for  the  Single 
Input-Single  Output  Systems  for  Flight 
Condition  Four  . 


List  of  Tables 


Table 

2-1  AFTI/F-16  Aircraft  Data  . 

5-1  Figure  of  Merit  Tolerances  for  Pitch  Rate 

5-2  Pitch  Rate  Response  Figures  of  Merit  .  . 

5-3  Estimated  Roll  Angles  with  Pitch 

Rate  Commanded  . 

5-4  Surface  Position  and  Rate  Limits  .... 

5-5  Figure  of  Merit  Tolerances  for  Roll  Rate 

5-6  Roll  Rate  Response  Figures  of  Merit  .  . 

5-7  Estimated  Pitch  Angle  with  Roll 

Rate  Commanded  . 

A-l  Aircraft  Data  for  Flight  Condition  One  . 

A-2  Aircraft  Data  for  Flight  Condition  Two  . 

A-3  *  Aircraft  Data  for  Flight  Condition  Three 

A-4  Aircraft  Data  for  Flight  Condition  Four 

C— 1  Equivalent  Plant  Matrices  for  Flight 

Condition  One  . 

C-2  Equivalent  Plant  Matrices  for  Flight 

Condition  Two  . 

C-3  Equivalent  Plant  Matrices  for  Flight 

Condition  Three  .  . 

C-4  Equivalent  Plant  Matrices  for  Flight 

Condition  Four  . 

C-5  Equivalent  Plant  Equations  for  the  Single 
Input-Single  Output  Systems  for  Flight 
Condition  One  . 


Figure 

G-6 (d) 

Control  Surface  Deflections  for  Roll 
Rate  Command,  Flight  Condition  4-- 
CSC  Mode  4  . 

G-6  (e) 

Control  Surface  Deflections  for  Roll 
Rate  Command,  Flight  Condition  4-- 
CSC  Mode  5  . 

G-6 (f ) 

Control  Surface  Deflections  for  Roll 
Rate  Command,  Flight  Condition  4 — 
CSC  Mode  6  . 

Page 


Figure 


G-4(d)  Control  Surface  Def lectio.  ^  for  Roll 
Rate  Command,  Flight  Condition  2 — 

CSC  Mode  4 . 200 

G-4  (e)  Control  Surface  Deflections  for  Roll 
Rate  Command,  Flight  Condition  2— 

CSC  Mode  5 . 201 

G-4(f)  Control  Surface  Deflections  for  Roll 
Rate  Command,  Flight  Condition  2-- 
CSC  Mode  6 . 201 

G-5  (a)  Control  Surface  Deflections  for  Roll 
Rate  Command,  Flight  Condition  3 — 

CSC  Mode  1 . 202 

G-5(b)  Control  Surface  Deflections  for  Roll 
Rate  Command,  Flight  Condition  3 — 

CSC  Mode  2 . 202 

G-5(c)  Control  Surface  Deflections  for  Roll 
Rate  Command,  Flight  Condition  3 — 

CSC  Mode  3 . 203 

G-5(d)  Control  Surface  Deflections  for  Roll 
Rate  Command,  Flight  Condition  3~ 

CSC  Mode  4 . 203 

G-5(e)  Control  Surface  Deflections  for  Roll 
Rate  Command,  Flight  Condition  3 — 

CSC  Mode  5 . 204 

G-5(f)  Control  Surface  Deflections  for  Roll 
Rate  Command,  Flight  Condition  3 — 

CSC  Mode  6 . • .  ’ . 204 

G-6(a)  Control  Surface  Deflections  for  Roll 
Rate  Command,  Flight  Condition  4 — 

CSC  Mode  1 . 205 

G-6(b)  Control  Surface  Deflections  for  Roll 
Rate  Command,  Flight  Condition  4 — 

CSC  Mode  2 . 205 

G-6  (c)  Control  Surface  Deflections  for  Roll 
Rate  Command,  Flight  Condition  4-- 
CSC  Mode  3 . 206 


Figure 
G-2  (c) 

G-2  (d) 

G-2 (e) 

G-2 (f ) 

G-3  (a) 

G-3 (b) 

G-3  (c) 

G-  3  ( d ) 

G-3 (e) 

G-3 (f ) 

G-4  (a) 

G-4 (b) 

G-4 (c) 


Control  Surface  Deflections  for  Pitch 
Rate  Command,  Flight  Condition  3-- 
CSC  Mode  3  . 

Control  Surface  Deflections  for  Pitch 
Rate  Command,  Flight  Condition  3 — 

CSC  Mode  4  . 

Control  Surface  Deflections  for  Pitch 
Rate  Command,  Flight  Condition  3-- 
CSC  Mode  5  . 

Control  Surface  Deflections  for  Pitch 
Rate  Command,  Flight  Condition  3 — 

CSC  Mode  6  . 

Control  Surface  Deflections  for  Pitch 
Rate  Command,  Flight  Condition  4 — 

CSC  Mode  1  . 

Control  Surface  Deflections  for  Pitch 
Rate  Command,  Flight  Condition  4-- 
CSC  Mode  2  . 

Control  Surface  Deflections  for  Pitch 
Rate  Command,  Flight  Condition  4 — 

CSC  Mode  3  . 

Control  Surface  Deflections  for  Pitch 
Rate  Command,  Flight  Condition  4-- 
CSC  Mode  4  . 

Control  Surface  Deflections  for  Pitch 
Rate  Command,  Flight  Condition  4-- 
CSC  Mode  5 . . 

Control  Surface  Deflections  for  Pitch 
Rate  Command,  Flight  Condition  4 — 

CSC  Mode  6  . 

Control  Surface  Deflections  for  Roll 
Rate  Command,  Flight  Condition  2 — 

CSC  Mode  1  . 

Control  Surface  Deflections  for  Roll 
Rate  Command,  Flight  Condition  2 — 

CSC  Mode  2  . 

Control  Surface  Deflections  for  Roll 
Rate  Command,  Flight  Condition  2-- 
CSC  Mode  3  . 


x 


two  flaperons  and  two  horizontal  tails  each  one  indi¬ 
vidually  controlled.  Control  surface  configurations  are 
as  follows: 


Mode  1  -  All  four  surfaces  are  operating  normally. 

Mode  2  -  One  horizontal  tail  is  failed. 

'  Mode  3  -  One  flaperon  is  failed. 

Mode  4  -  One  horizontal  tail  and  one  flaperon  are 
failed  (same  side) . 

Mode  5  -  One  horizontal  tail  and  one  flaperon  are 
failed  (opposite  sides) . 

Mode  6  -  Both  flaperons  are  failed. 

Throughout  this  thesis,  modes  1  through  6  are  referred  to 
as  control  surface  configuration  (CSC)  modes,  and  a  case 
refers  to  a  particular  flight  condition  and  CSC  mode. 

Also,  for  the  purposes  of  this  study,  a  failed  surface  is 
defined  to  be  fixed  in  the  neutral  reference  position. 

A  large  degree  of  parameter  uncertainty  is 
incorporated  into  the  design  by  choosing  four  design  points 
representing  a  large  portion  of  the  flight  envelope.  These 
design  points  are  seen  in  Figure  1-1. 


I . 3  Assumptions 

The  aircraft  equations  of  motion  were  developed  by 
Mr.  Finley  Barfield  using  the  assumptions  listed  below  (6). 

1.  The  aircraft  is  a  rigid  body,  and  mass  is 
constant. 
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Fig.  1-1.  Flight  Envelope  Design  Points 

2.  Thrust  is  not  changing. 

3.  The  earth's  surface  is  an  inertial  reference 

frame . 

4.  The  atmosphere  is  assumed  fixed  with  respect 
to  the  earth. 

5.  The  equations  of  motion  can  be  decoupled  into 
longitudinal  and  a  lateral-directional  set. 

6.  Linearization  about  an  operating  condition  is 
acceptable  for  point  designs. 

7.  Aerodynamics  are  fixed  for  Mach  number  and 
altitude . 

A  detailed  explanation  of  each  of  the  above  assumptions  is 
presented  in  Reference  6. 
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In  addition  to  the  above  assumptions,  the  lateral 


and  longitudinal  control  derivatives  for  the  elevators 
and  flaperons  are  halved  in  order  to  obtain  control  coeffi¬ 
cients  for  the  individual  control  surfaces.  This  approach 
is  assumed  suitably  accurate  for  a  preliminary  design  study 
of  this  nature. 

I. 4  Approach 

The  first  step  in  the  design  process  is  to  derive 
the  basic  plant  transfer  functions.  These  are  obtained 
from  the  aircraft  equations  in  state-variable  form,  with 
the  longitudinal  and  lateral  control  elements  coupled  in 
the  input  matrix.  From  the  system  block  diagram,  equiva¬ 
lent  plant  equations  are  derived  for  each  of  the  previously 
described  surface  configurations.  As  described  in  Chapter 
III,  the  four  input- two  output  system  is  transformed  into 
an  equivalent  two-by-two  system.  The  resulting  two-by-two 
equivalent  plant  transfer  function  matrices,  £e,  are  then 
inverted  to  yield  reciprocal  expressions,  P  =  (1/Q. .), 
for  equivalent  single  input-single  output  (SISO)  systems. 
Next,  output  specifications  are  formulated. 

The  pitch  rate  and  roll  rate  specifications  are 
first  developed  in  the  time  domain,  and  then  transformed 
into  the  frequency  domain  via  the  Laplace  transformation. 
Designing  in  the  frequency  domain  avoids  the  more  difficult 
convolution  operations  of  time  functions.  The  design  speci¬ 
fications  serve  mainly  to  demonstrate  the  application  of 


V 


the  technique.  They  are  based  upon  what  the  author  feels 
are  reasonable  from  his  personal  experience  gained  from 
flying  fighter  aircraft,  and  serve  primarily  as  a  basis 
for  demonstrating  QFT.  Next,  the  compensators  are  synthe¬ 
sized  according  to  the  techniques  of  Reference  3. 

'  The  longitudinal  and  lateral  controller  designs 
are  incorporated  into  the  complete  system,  and  the  design 
is  verified  by  simulations  at  each  flight  condition  for 
each  control  surface  configuration.  Control  surface  deflec¬ 
tions  are  also  examined  for  position  and  rate  saturation. 

I. 5  Presentation 

This  thesis  consists  of  six  chapters.  Chapter  II 
contains  the  aircraft  equations  of  motion  in  state  vari¬ 
able  format  that  are  used  to  derive  the  basic  plant  transfer 
function  equations.  Chapter  III  presents  the  system  model 
used  in  the  simulation  and  the  quantitative  feedback  synthe¬ 
sis  theory  applied  to  the  reconfiguration  problem.  The 
actual  design  of  the  control  laws  is  presented  in  Chapter 
IV.  Chapter  V  displays  the  output  and  control  surface 
responses,  and  evaluations  of  such.  Chapter  VI  contains 


the  conclusions  and  recommendations. 


II .  Aircraft 


11. 1  Introduction 

In  this  chapter  the  decoupled  longitudinal  and 
lateral  aircraft  equations  of  motion  in  state-variable 
format  are  coupled  together  in  the  input  matrix  allowing 
control  of  individual  elevators  and  flaperons.  From  these 
coupled  aircraft  equations,  the  basic  plant  transfer  func¬ 
tions  of  interest  are  computed.  The  linearized,  aero¬ 
dynamic  coefficients  used  to  develop  the  equations  of  motion 
are  developed  by  the  General  Dynamics  Corporation,  and  are 
taken  from  Reference  6.  A  description  of  the  AFTI/F-16  and 
its  flight  control  system  is  also  given  in  Reference  6. 

The  aerodynamic  coefficients  used  in  the  equations  are 
listed  in  Appendix  A  for  four  flight  conditions.  Table  2-1 
presents  pertinent  aircraft  data  for  the  AFTI/F-16. 

11. 2  Aerodynamic  Model 

The  equations  of  motion  are  written  in  state-space 
form,  x  =  Ax  +  Bu.  They  are  a  set  of  first-order  differen¬ 
tial  equations  which  represent  the  aircraft  dynamics.  The 
elements  of  the  A  and  B  matrices  are  dimensionalized  body 
axis  coefficients.  Coefficients  in  the  stability  axis 
system  are  converted  to  the  body  axis  system  using  a  com¬ 
puter  program  called  Conversion  and  Transformation  (6) . 


TABLE  2-1 


AFTI/F-16  AIRCRAFT  DATA 

Aircraft  Parameters 

Wing  Mean  Aerodynamic  Chord,  c(ft)  =  11.32 

2 

Wing  Reference  Area,  s(ft  )  =  300.0 
Wing  Span,  b(ft)  =  30.0 
Weight,  w(lbs)  =  21,018.0 

Moment  of  Inertias 

I  (slug-ft2)  =  10,033.4 

I  (slug-ft2)  =  53,876.3 

I  (slug-ft2)  =  61,278.4 

z  z 

I  (slug-ft2)  =  282.132 

X  z 

The  body  axis  is  chosen  as  the  design  axis  for  two  primary 
reasons.  First,  the  sensors  providing  feedback  signals 
actually  measure  body  accelerations  and  rates.  Secondly, 
the  pilot  desires  to  control  accelerations  and  rates  that 
he  actually  feels  which  are  body  accelerations  and  rates  (6) 
Using  the  assumptions  of  Chapter  I,  and  assuming 
that  coupling  can  be  ignored  between  the  longitudinal  and 
lateral  axes,  two  sets  of  linear,  first-order  differential 
equations  result.  The  longitudinal  set  of  equations  has 
the  following  form: 


Using  the  computer  program  Total,  transfer  func¬ 
tion  expressions  for  the  desired  output- to- input  rela¬ 
tionships  are  computed  (7) .  These  transfer  functions 
are  referred  to  as  the  basic  plant  equations.  They  are 
identified  by  a  two-digit  subscript,  and  a  superscript 
(a  or  b) .  The  first  digic  in  the  subscript  denotes  the 
output  while  the  second  digit  denotes  the  control  surface 
input  that  generates  the  output.  Pitch  and  roll  rates  are 
identified  by  a  1  or  a  2 ,  respectively.  Likewise,  the 
elevators  and  flaperons  are  identified  by  a  1  or  a  2, 
respectively.  A  superscript  "a"  or  "b"  denotes  a  right 
or  left  control  surface,  respectively.  For  example, 

P12a  identifies  the  basic  plant  equation  whose  output  is 
pitch  rate  generated  by  the  right  flaperon.  Combining 
the  longitudinal  and  lateral  matrices,  and  dividing  con¬ 
trol  between  individual  elevator  and  flaperon  surfaces 
results  in  a  set  of  equations  with  coupling  between  the 
longitudinal  and  lateral  axes  in  the  control  matrix 
(Eq.  2-3). 
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Plant  expressions  are  computed  for  each  of  the 


flight  conditions,  and  are  listed  next.  Due  to  symmetry. 


Flight  Condition  I.  (Mach  =  0.2  at  30  ft) 


„  a  _  q  -1.118s(s+0. 01822) (s+0. 4568) 

11  6  a  ”  (s-0.3633) (s+1.3) (s+0. 07683+j0. 2065) 


P10a=-3-=  .1209s (s+0. 06537) (s+0. 2589)  (2-5) 

62a  (s-0.3633) (s+1.3) (s+0. 07683+ jO. 2065) 


,  a  _  _2_  _ -2 ,239s  (s+0 . 205+ j0  ■  853) _ 

22  6a  (s+0. 1041) (s+0. 6835) (s+0. 2741±jl. 909) 

2 


.a  p  -2.142s (s+0. 3017±il. 562) 

21  “  .  a  (s+0. 1041) (s+0. 6835) (s+0. 2741+ jl. 909) 


Flight  Condition  2  (Mach  =  0.6  at  30,000  ft) 


a  -2. 931s (s+0. 01004) (s+0. 5502) 

P11  (s-1.167) (s+2.028) (s+0. 006472+ jO. 07803) 


,  a  _  -0.1059s (s+0. 006697) (s+1.861) 

12  “  (s-1.167)  (s+2.028) (s+0. 006472+j0. 07803) 


(2-10) 


p  a  _  -8. 723s (s+0. 2191 jl. 607) _ 

22  (s+0. 07795)  (s+0.8265)  (s+0.211+jl. 953) 


a  _  -6.792s (s+0. 2442±j 2. 101) 

21  (s+0. 07795)  (s+0.8265)  (s+0. 211± jl. 953) 


Flight  Condition  3  (Mach  =  0.9  at  20,000  ft) 


p  a_  -12.03s (s+0. 01262) (s+1. 51) 

11  (s-0.9645) (s+3.223) (s+0. 007553+ jO. 05384) 

a  -3. 236s (s+0. 01254)  (s+1. 646) 

12  (s-0.9645) (s+3.223) (s+0. 0075531 jO. 05384) 

p  a  -25. 53s (s+0. 3541+j2. 927) 

22  (s+0. 02719)  (s+2.697) (s+0. 391± j2. 962) 

„  a  _  -25. 36s (s+0. 3749±i3. 578) 

21  (s+0. 02719)  (s+2.697)  (s+0.391±j2.962) 


Flight  Condition  4_  (Mach  =  1.6  at  30,000  ft) 


„  a  _  -16. 45s(s+0. 02996)  (s+1.097) 

11  (s+0. 015161j0. 02343) (s+O.8012±j6. 592) 


,  a  -2. 925s(s+0. 03459) (s+0. 6861) 

12  (s+0. 01516±j0. 02343) (s+0. 8012+ j6.592) 


a  _  -7. 084s (s+0. 408314. 916) 

*22  (s+0. 03448) (s+2.171) (s+0. 49961 j3. 129) 


,  a  -23.3s(s+0.3774lj3.848) _ 

21  (s+0. 03448)  (s+2.171) (s+0. 49961 j3. 129) 


(2-11) 

(2-12) 

(2-13) 

(2-14) 

(2-15) 

(2-16) 

(2-17) 

(2-18) 


(2-19) 


II. 3  Summary 

The  aircraft  equations  of  motion  in  state-variable 
format  are  presented  in  this  chapter.  The  longitudinal  and 
lateral  axes  are  coupled  within  the  control  matrix  allowing 
control  of  individual  elevators  and  flaperons.  Also  pre¬ 
sented  are  the  basic  plant  transfer  functions  of  interest 
computed  from  the  coupled  aircraft  equations  of  motion  for 
each  of  the  four  flight  conditions. 

The  next  chapter  presents  the  reconfiguration  sys¬ 
tem  model,  and  quantitative  feedback  theory  as  applied  to 
the  reconfiguration. 


III.  Reconf igura t ion  Theory 

111.1  Introduction 

The  concepts  of  QFT  can  be  used  to  design  fixed 
compensators  to  achieve  tolerance  specifications  (2;  3) . 

The  specifications  are  to  be  satisfied  despite  highly 
uncertain  parameters  which  can  arise  from  aircraft  control 
surface  failures  (partial  or  total) ,  and  from  widely  vary¬ 
ing  flight  conditions.  The  specif ications ’must  also  be 
satisfied  for  the  normal  system  (no  failures) .  The  advan¬ 
tage  of  a  fixed  compensator  design  is  that  fast  identifica¬ 
tion  of  failures  is  not  necessary.  A  more  optimum  compen¬ 
sator  can  be  switched  in  later,  if  required,  after  the 
failure  state  has  been  explicitly  identified.  This  chapter 
presents  the  MIMO  system  model,  and  develops  an  equivalent 
plant 'matrix  for  which  QFT  is  applicable. 

111. 2  Reconfiguration  Theory  and 

System  Model 

The  structure  of  the  system  model  is  seen  in 
Figure  3-1,  and  represents  a  four  input-two  output  system. 

A  two-elevator,  two-flaperon  system  is  represented  with 
each  individual  surface  separately  controllable.  The 
elevators  and  flaperons  are  capable  of  moving  individually, 
and  symmetrically  or  differentially  in  pairs.  The  right 
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.  3-1.  Plant  and  Compensation 
Model  Signal  Flow  Diagram 


Fig.  4-5.  Equivalent  SISO  Systems 


at  a)  *  0.1  rad/sec  are  seen  in  Figures  4-6  and  4-7,  respec¬ 
tively.  The  q  channel  is  represented  by  Q^,  and  the  p 
channel  by  Q22.  Each  template  represents  the  variable 
uncertainty  at  a  particular  frequency  for  the  six  CSC  modes 
and  four  flight  conditions.  Thus,  each  template  has  24 
points;  the  modes  are  numbered  1  through  6  as  indicated  on 
the  template  figures.  Each  CSC  mode  is  repeated  four  times 
at  a  particular  frequency  corresponding  to  the  four  flight 
conditions.  A  nominal  condition  is  chosen  and  indicated 
by  circling  that  condition  on  each  template.  The  nominal 
condition  for  is  flight  condition  2,  CSC  mode  4  (Mach  = 
0.6;  6^a  and  <S2a  failed),  while  flight  condition  4,  CSC 
mode  4  (Mach  =  1.6)  is  selected  as  the  nominal  condition 


TABLE  4-1 

FIGURE  OF  MERIT  TOLERANCES  FOR  PITCH  RATE 


Settling  Time 
tg  (sec) 

Rise  Time 
tr  (sec) 

Peak  Value 

M  (deg/sec) 

1.5  -  3.0 

0.33  -  1.63 

<  11 

FIGURE 

TABLE  4-2 

OF  MERIT  TOLERANCES  FOR 

ROLL  RATE 

Settling  Time 
ts  (sec) 

Rise  Time 
tr  (sec) 

Peak  Value 
Mp  (deg/sec) 

3.0  -  6.0 

0.65  -  3.23 

<  55 

Transfer  function  models  for  roll  rate  are  also  listed  in 
Appendix  D.  Once  acceptable  tolerances  are  established, 
bounds  on  the  loop  transmissions  are  formed. 

IV. 4  Loop  Bound  Design,  L2q 

The  decision  must  be  made  on  which  loop  transmis¬ 
sion  (L)  is  to  be  designed  first.  Since  the  equivalent 
system  is  a  two-by-two  structure  with  two  inputs  to  con¬ 
sider,  both  the  command  and  the  BNIC  elements  must  be  ana¬ 
lyzed  for  each  L  to  define  the  loop  bounds.  The  equivalent 
SISO  single-loop  systems  are  seen  in  Figure  4-5. 

As  to  which  L  to  design  first,  both  loops  require 
the  same  analysis,  so  neither  is  less  complicated  than 
the  other.  For  this  problem,  the  L  for  the  p  channel  is 
designed  first  because  of  the  smaller  uncertainty  associ¬ 
ated  with  this  variable.  This  is  seen  by  comparing  the 
plant  templates  for  p  and  q.  Templates  of  Q. .  and 


Fig.  4-4.  Roll  Rate  Tolerances 
in  the  Frequency  Domain 
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Fig.  4-2.  Pitch  Rate  Tolerances  in  the  Frequency  Domain 

merit  are  given  in  Reference  10.  Tolerances  for  these 
figures  of  merit  are  seen  in  Table  4-1.  The  p  tolerance 
b21  is  modeled  for  a  total  roll  angle  of  approximately 
five  degrees.  This  specification  is  established  with  severe 
failure  situations  in  mind,  such  as  one  elevator  and  one 
flapelron  failed  on  the  same  side.  Five  degrees  of  roll 
angle  seems  quite  reasonable  for  situations,  such  as  this. 
Tolerances  are  also  established  for  the  p  command  responses. 

Roll  rate  tolerances  b22,  a22 '  an<^  b12  are  also 
established  as  is  done  for  the  q  tolerances.  It  is  desired 
that  the  p  commanded  be  50  deg/sec  with  small  q  resulting. 
The  q  tolerance  b^2  is  modeled  for  a  total  pitch  angle 
of  approximately  three  degrees.  Figures  4-3  and  4-4 
illustrate  the  time  and  frequency  domain  tolerances  for  p, 
and  Table  4-2  lists  the  figure  of  merit  tolerances  for  p. 
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Fig.  4-1.  Pitch  Rate  Tolerances  in  the  Time  Domain 

The  p  response  is  a  result  of  loop  interaction,  and  is 
denoted  by  b2^.  The  time  response  tolerances  transformed 
into  the  frequency  domain  are  seen  in  Figure  4-2.  As  dis¬ 
cussed  earlier,  all  design  work  is  accomplished  in  the  fre¬ 
quency  domain.  Transfer  function  models  for  the  response 
tolerances  are  listed  in  Appendix  D.  Figures  of  merit  are 
established  for  the  response  tolerances.. 

Figures  of  merit  for  the  outputs  of  interest  are 
established  primarily  on  what  the  author  feels  are  reason¬ 
able  from  his  personal  experience  gained  from  flying  fighter 
aircraft.  The  main  concern  is  to  demonstrate  the  efficacy 
of  QFT  when  applied  to  a  reconfiguration  problem.  Settling 

time  t  ,  rise  time  t_,  and  peak  value  are  chosen  as  the 
s  r  c  p 

design  figures  of  merit.  Definitions  for  the  figures  of 
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The  second  condition. 


realizable  for  that  specific  Pg. 
also  derived  in  Reference  8,  requires  that 


Q 

Q 


11Q22 

12Q21 


<  1 


as  s  -*■  °° 


(4-1) 


This 'must  be  so  for  all  Pg  £  P,  and  all  u  >  some  where 
|1+L|  <1.  This  second  condition  is  an  inherent  constraint. 
Both  of  the  above  conditions  for  this  thesis  are  found  to 
be  satisfied  for  every  P  ,  so  design  can  proceed. 

IV. 3  Specifications 

In  general,  a  control  problem  exists  because  a  set 
of  performance  criteria  are  specified  (10).  Acceptable 
tolerances  for  the  closed-loop  responses  are  established, 
and  the  goal  is  to  find  compensation  which  guarantees  that 
these  tolerances  are  satisfied  over  the  range  of  plant 
uncertainty.  In  QFT  it  is  guaranteed  that  the  amount  of 
feedback  designed  into  the  system  is  such  that  the  desired 
responses  are  obtained  despite  highly  uncertain  plant 
parameters  (2) .  Compensation  is  synthesized  so  that  the 
specifications  are  satisfied  for  all  design  flight  condi¬ 
tions  and  CSC  modes. 

The  tolerances  for  the  q  and  p  responses  of  Figure 
3-1,  when  q  is  commanded,  are  illustrated  in  Figure  4-1. 

It  is  desired  that  the  q  commanded  be  ten  degrees  per  second 
(deg/sec)  with  small  p  resulting.  The  upper  bound  for 


q  response  is  denoted  by  b^  and  the  lower  bound  by  a^. 


and  CSC  mode.  For  this  study  the  equivalent  plant  set  has 
as  members  the  matrices  for  each  flight  condition  and  CSC 
mode,  giving  a  total  of  24  matrices  (4  flight  conditions 
times  6  CSC  modes) .  The  elements  of  each  Pe  are  listed 
in  Appendix  C.  After  developing  the  P  ,  £e_1  is  computed 
using  the  computer  program  Total,  obtaining  P  ^  =  (1/Q.  .). 

6  1 J 

The  become  the  plant  elements  for  the  equivalent  SISO 
systems  (2;  3) .  Appendix  C  also  lists  the  individual 
for  each  flight  condition  and  CSC  mode.  Notice  in  Appen¬ 
dix  C  that  some  of  the  Q's  have  unstable  poles.  This  does 
not  present  an  obstacle  for  this  technique.  This  is 
addressed  later  in  this  chapter.  Zeros  in  the  right-half- 
plane,  however,  represent  a  non-minimum  phase  system,  and, 
as  a  result,  the  design  is  more  difficult.  Dr.  Horowitz 
has  developed  QFT  to  design  for  non-minimum  phase  plants 
(9) .  Since  this  thesis  does  not  involve  non-minimum  phase 
plants,  they  are  not  addressed.  Before  proceeding  with  the 
design,  certain  conditions  must  first  be  satisfied  for  the 
MIMO  synthesis  technique  to  be  applicable. 

Two  conditions  must  be  satisfied  before  proceeding 
(8).  First,  Po  1  must  exist  in  order  to  obtain  the 
elements  for  the  single- loop  systems.  This  is  the  "con¬ 
trollability"  condition  and  must  be  true  for  each  flight 
condition  and  CSC  mode.  If  the  determinant  of  P£  equals 
zero  it  means  that  the  plant  outputs  are  not  independently 


IV.  1 


Introduction 


This  chapter  presents  the  steps  used  in  the  design 
of  the  compensators  g^,  g2  in  Figure  3-1.  The  procedure  is 
started  by  inverting  the  elements  of  the  equivalent  plant 
transfer  function  matrices  for  the  various  CSC  modes  to 
yield  the  plant  equations  for  the  single-loop  systems 
(see  Figure  4-5,  Section  IV. 4).  These  plant  equations 
are  used  in  the  detailed  design  execution.  The  p  channel 
compensator  is  designed  first  due  to  a  smaller  amount  of 
uncertainty  associated  with  this  variable.  Using  data 
derived  in  the  design  of  this  loop,  the  q  channel  com¬ 
pensator  is  then  derived  using  the  "improved  method"  (8) . 
Simple,  first-order  prefilters  are  designed  for  both 
channels.  The  synthesis  begins  with  the  plant  equations. 

IV. 2  Plant  Equations 

The  design  process  begins  by  obtaining  the  basic 
plant  transfer  function  equations  from  Eq.  (2-3)  for  each 
design  flight  condition.  These  equations  are  listed  in 
Chapter  II.  The  equivalent  plant  transfer  function 
matrices  (F^)  are  developed  as  described  in  Chapter  III  from 
these  basic  plant  equations.  Equivalent  plant  transfer 
function  matrices  are  developed  for  each  flight  condition 


I I I. 4  Summar 


Presented  in  this  chapter  is  the  multiple  input- 

multiple  output  system  model  and  reconfiguration  theory. 

An  equivalent  square  plant  matrix  is  also  developed  from 

a  non-square  plant  matrix  so  that  quantitative  feedback 
0 

theory  QFT  can  be  applied.  In  the  next  chapter,  the 
compensator  design  procedures  and  the  results  are  presented. 


Substitutions  are  made  into  Eq.  (3-7)  for  the 


various  control  surface  configurations  producing  closed- 
loop  transfer  functions  for  each  command  and  cross  coupling 
response.  For  example,  consider  the  case  in  which  the 
right  elevator  has  failed  (S^3  =  0)  ;  this  implies  that  P^a 
and  P2i3  are  ec3ual  to  zero-  With  these  substitutions,  the 
equivalent  plant  matrix  becomes 


P 

— e 


Pllb+vlP12 

-p2ib+yAP22 


Apl2+ypn 

p22+yp2ib 


(3-8) 


This  new  equivalent  plant  matrix  is  used  to  deter¬ 
mine  the  four  t^j  transfer  functions.  Each  control  surface 
failure  situation  can  be  treated  in  this  manner  to  obtain 
the  resulting  plant  transfer  function.  The  various  failures 
give  different  equivalent  plant  matrices  for  each  flight 
condition.  Quantitative  feedback  theory  is  applied  to  the 
entire  set.  The  number  of  elements  in  the  set  is  the 
product  of  the  number  of  CSC  modes  considered  by  the  number 
of  flight  conditions.  Appendix  C  lists  the  entire  Pfi  set 
used  in  this  thesis.  QFT  is  well  suited  for  this  type  of 
problem  since  the  surface  failures  merely  enlarge  the  plant 
set.  The  extra  "cost  of  feedback"  can  readily  be  seen  by 
the  designer. 


flgl  yf2g2 

gl  yg2 

^1?!  ”^^2g2 

rl 

gl  "yg2 

*i 

yflgl  f2g2 

r2 

ygi  g2 

_y2_ 

y^lgl  ~^2g2 

_ygi  "g2 

From  Eqs.  (3-1)  to  (3-5)  ,  and  by  letting  P„  = 

P.  .  +  P.  .  ,  and  AP.  •  =  P-  •  -  P .  .  ,  an  equivalent  plant 

1 3  13  13  x3  13 

matrix  Pg/  emerges  as  follows: 


pU+ypl2 

AP12+ViAPii 

P11 

e 

P12 

e 

AP21+yAP22 

P22+yP21 

*21e 

P22 

e 

Now  the  plant  matrix  is  square,  and  QFT  is  directly 
applicable.  This  two-by-two  MIMO  system  can  be  reduced 
to  four  SISO  single-loop  systems  for  the  design  synthesis 
(2;  3). 

After  completing  all  the  matrix  operations  the 
final  expression  for  T  becomes 

T  -  (ti:j} 

(?2P22  +1) -flglg2P12  P21  ^  ^f2g2P'1''  ^ 

e  e  e  e 

{flglP21  }  {f2g2P22  (glPll  +1)_f2glg2P12  P21  J 

= - ^ ^ ^ - - - -  (3-7) 

(glPH  +1)  (?2P22  +1*-flf2glg2P12  P21 
e  e  e  e 
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are  constants  which 


failure.  In  Figure  3-1,  and  V^i 
determine  the  division  of  control  effort  between  the 
flaperons  and  elevators.  For  this  preliminary  study, 
optimization  of  these  constants  is  not  addressed,  and  both 
u's  are  set  equal  to  0.25  for  the  majority  of  the  simula¬ 
tions.  Henceforth,  the  subscripts  are  dropped  on  these 
two  constants.  Preliminary  theory  has  been  developed  to 
enable  the  designer  to  specify  the  optimum  division  of 
control  effort.  Appendix  B  contains  this  work  along  with 
a  detailed  explanation  of  reconfiguration  theory  for 
larger,  more  complex  systems.  The  next  section  presents 
the  development  of  an  equivalent  plant  for  which  QFT  is 
applicable. 


III. 3  Equivalent  Plant  Development 

In  Figure  3-1,  Y  =  P6  =  P{GFr  -  GY} 


Thus, 

T  =  { t± j }  =  {I  +  PG }~1  PGF 
where  t^  =  J_I  +  PG)  1. 


P 


6 


a 

„  b 

_  a 

„  b 

11 

P11 

P12 

P12 

a 

_  b 

_  a 

„  b 

21 

_P21 

P22 

“P22 

[6xa,  6xbr  <S2a,  62b] 


(3-1) 


(3-2) 


(3-3) 


(3-4) 
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f 


respec- 


3  b 

and  left  elevators  are  designated  by  6^  and  6^ 
tively;  the  flaperons  are  designated  similarly  by  62a  and 
8 2^.  The  longitudinally  controlled  output  variable  is 
pitch  rate  (q) ,  while  the  laterally  controlled  output  vari¬ 
able  is  roll  rate  (p) .  As  stated  in  Chapter  II,  the  plant 
transfer  function  designator,  relates  the  ith  output 

to  the  jth  input  (control  surface) .  A  superscript  "a"  or 
"b"  denotes  the  right  or  left  control  surface.  For  example, 
P2^  denotes  the  effect  of  the  right  elevator  (6^  )  on 

roll  rate  (output  y2) .  For  a  normally  operating  system, 

_ 

P. .  =  P. .  .  Minus  signs  are  associated  with  P„.  and  P„„ 
lj  i]  ^  21  22 

so  that  the  two  pairs  of  surfaces  move  differentially  to 
generate  roll  rate.  In  the  no  fail  case,  surfaces  move  an 
^  «  equal  amount  in  either  direction  to  generate  the  commanded 

rates.  No  interaction  is  desired  between  the  longitudinal 
and  lateral  modes  under  normal  operating  conditions. 

Normally,  the  elevators  move  together  in  the  same 
direction  to  excite  only  the  longitudinal  mode,  and  the 
flaperons  move  opposite  each  other  to  excite  only  the 
lateral  mode.  Separate  control  of  the  surfaces  increases 
the  flexibility  of  a  system  by  exciting  both  aerodynamic 
modes  when  a  single  surface  is  actuated.  For  instance, 
using  the  structure  of  Figure  3-1,  the  surfaces  can  con¬ 
tinue  to  act  in  pairs  under  normal  operation,  but  act 
separately,  as  necessary,  to  offset  parameter  changes 
resulting  from  single  or  multiple  control  system  component 
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MAGNITUDE  ( DB  ) 

-40.00  -20.00  0.00  20.00  40.00 


Q22  (F  =  0.1) 


-360.00  -300.00 


■240.00  -180.00 

DEGREES 


(FAILURES) 

1  -  NONE 

2  -  DELTA  1A 

3  -  DELTA  2A 
4--  DELTA  1A»2A 

5  -  OELTA  l A»2B 

6  -  OELTA  2A»2B 

-1  20 . 00 


Fig.  4-7.  Template  of  Qj2 
at  w  =  0.1  Radians /Second 


A  '»  .%  V  /. 


for  Q22«  Dr’  Horowitz  recommends  using  the  same  nominal 
condition  for  each  loop  design.  This  eliminates  any  con¬ 
fusion  as  to  the  nominal  condition  when  solving  for  the 
compensator  explicitly.  Templates  are  made  for  frequencies 
approximately  one  octave  apart  until  the  locus  of  points 
is  essentially  a  vertical  line.  Appendix  E  contains 
templates  of  and  Q22  at  several  frequencies.  An 
analysis  of  the  p  channel  is  accomplished  next. 

The  BNIC  system  whose  output  is  y2^  is  analyzed 
first  since  only  the  input  due  to  loop  interaction  is 
considered.  A  response  due  to  r.^  does  not  exist  for  this 
single-loop  system  since  the  prefilter  f21  is  chosen  to  be 
zero.  Consider  the  expression  for  y2^ 


d21Q21 


21  l+g2Q22 


(4-2) 


where 


l21  ~  tll'/Q21^2^'  and  L2  ~  g2°22 


(4-3a  ,b) 


Rewriting  Eq.  (4-2) 
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(tll/Q21)Q22 

1+L„ 


-  b21 


(4-4) 


with  b21  being  the  specification  bound  for  p  due  to  the  q 
command.  To  ensure  total  compliance  for  all  acceptable 
til'  tbe  uPPer  bound  of  bll  for  the  q  response  becomes  bll* 
Upon  rearranging  terms,  Eq.  (4-4)  becomes 


In  order  to  determine  the  most  demanding  bound,  use  the 
minimum  value  of  Q2^  ^or  flight  conditions  and  CSC 

modes  at  the  frequency  in  question.  Also,  use  the  nominal 


value  for  Q22 '  Q22 Q"  This  procedure  results  in  overdesign 
unless  | |  >>  1.  If  | Lj |  >>  1  is  not  a  valid  assumption, 
then  Q21  and  Q22  in  Eq.  (4-5)  should  correspond  to  the  same 
flight  condition  and  CSC  mode.  This  takes  into  account 


the  correlation  between  Q21  and  Q22»  The  bound  on  L2q  is 
then  formed  from  the  most  demanding  of  all  the  bounds  from 
Eq.  (4-5)  at  each  phase  angle.  For  this  particular  problem 
the  demands  on  L2q  are  high  in  the  low  frequency  region 
making  |l2|  >>  1  a  valid  assumption.  Equation  (4-6)  is 
used  to  locate  bounds  on  L2q  up  to  to  =  2  r/s.  At  higher 
frequencies  another  technique,  described  later,  is  used. 
Equation  (4-5)  now  becomes 


b21°21  . 
_ mm 

bll°22o 


(4-6) 


The  resulting  bound  is  the  bound  for  the  nominal  L2,  I<2o. 
The  terms  on  the  right  hand  side  of  Eq.  (4-6)  are  con¬ 
stants  obtained  from  the  Bode  plots  of  their  transfer 
functions.  In  the  region  where  |1+L2q|  I  | L2o I r  Eq*  (4-6) 


becomes 


L2ol  i 
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b21°21 


min 


(4-7) 


and  the  demand  on  L2q  is  obtained  explicitly.  Equation 
(4-7)  results  in  a  contour  on  the  modified  Nichol's  Chart 
(MNCJ ,  also  referred  to  as  the  inverted  Nichol's  Chart,  and 
is  compared  with  an  analysis  of  the  y22  loop  of  Figure  4-5 
to  produce  the  most  demanding  bound  on  L2q.  A  detailed 
explanation  of  the  MNC  is  presented  in  Reference  2.  Using 
Eq.  (4-7)  ,  at  03  =  0.1  r/s  L2q  is  calculated  as  follows: 


|l2  |  >  0.172+27- (-26.  9-22. 5)  -  76  decibels  (dB) 


Parallel  vertical  lines  are  used  to  indicate  the  log  magni¬ 
tude  (LM)  operation.  The  bound  on  L2q  at  w  =  0.1  r/s 
coincides  with  the  horizontal  line  for  L2q  because  the 
(1+L2q)  1  contour  at  -76  dB  is  also  a  horizontal  line. 

The  y22  loop  in  Figure  4-5  has  two  inputs  to  con¬ 
sider,  r2  and  d22,  and  therefore  two  transfer  function 
expressions  must  be  analyzed.  The  expression  for  t22  is 


t 


22 


22 


(4-8) 


where 

=  £22q2Q22  f22L2 

22  1+q2<“!22  1+L2 


(4-9) 


Using  the  value  of  L2q  found  in  the  analysis  of  y2i'  x^ 
is  determined  from  the  expression 


b12°22o 

L2oQ21 


(4-11) 


Next,  the  template  of  q22  for  the  frequency  in  question  is 

taken  and  placed  on  the  Nichol's  Chart  to  determine 

J T  22  /t22  I*  Eq.  (4-12),  below,  is  satisfied,  then 

max  min 

the  bound  found  in  the  first  analysis  satisfies  both  loops 
and  is  used  as  the  bound  on  L2q  for  that  one  frequency. 


(b22“a22) 


+  x. 


If  Eq.  (4-12)  is  not  satisfied,  then  the  template  of  Q22 
is  translated  vertically  until  it  is  satisfied.  The  con¬ 
tour  through  the  Q22o  point  becomes  the  new  bound  on  L2q. 
For  this  thesis  it  is  found  that  the  demands  from  y21  are 
the  most  demanding  and,  as  a  result,-  determine  the  bounds 
on  L2q  for  uk  such  that  |1+L2q|  ~  Jl2q|.  fre" 

quencies  this  approximation  is  not  valid,  and  another 
method  is  used  to  locate  the  bounds  on  L2q. 

The  lower  specification  bound  a22  is  less  than 
-12  dB  for  u)  >  2  r/s.  The  LMf22  is  also  relatively  small 
in  the  region  since  it  is  sloping  down  at  approximately 
the  same  rate  as  a5_.  In  Eq.  (4-8)  the  small  magnitude  of 


f~~  makes  x_0  negligible  in  comparison  to  t , 


This 


leaves  x.  and  x,  to  determine  the  bounds  on  L~  for 
d22  d21  20 
w  >  2  r/s.  It  is  found  that  the  -4  dB  contour  satisfies 

the  constraints  imposed  by  both  x ,  and  x  ,  in  this 

a22  a21 

region.  This  greatly  simplifies  the  task  of  determining 
bounds.  The  template  of  Q22  is  translated  down  and  hori¬ 
zontally  as  necessary  on  the  Nichol's  Chart  so  that  the 

template  remains  outside,  but  as  close  as  possible  to  the 
-4  dB  contour.  The  trace  of  the  nominal  point  during  this 
procedure  becomes  the  bound  on  L2q.  procedure  is 

repeated  for  several  frequencies  until  the  universal  high 
frequency  bound  is  formed  (2) .  Figure  4-8  shows  L2q  and 
its  bounds  plotted  on  the  Nichol's  Chart,  and  the  transfer 
function  expression  for  L2q  is 


J2o 


5.95(10°)  (s+0. 3)  (s+12 . 9)(s+120) 
s  (s+0 . 2)  (s+3.46)  (s+48)  (s2+900s+562 ,500) 


(4-13) 


Inspection  of  Figure  4-8  shows  that  L2q  is  close  to  its 
bounds  for  the  lower  frequencies,  and  avoids  the  -4  dB 
contours  at  high  frequency.  If  L^0  were  optimally  designed, 
it  would  lie  exactly  on  its  bound  at  each  frequency,  and 
the  loop  transmission  bandwidth  would  be  minimum.  The  band¬ 
width  of  the  loop  transmission  is  an  important  considera¬ 
tion,  and  if  too  large,  the  system  can  be  susceptible  to 
noise.  Using  Total,  the  bandwidth  of  L2q  is  found  to  be 
34  r/s.  No  attempt  to  analyze  the  system  for  susceptibility 


is  derived. 


After  L^0  is  determined,  the  prefilter 

Using  the  methods  of  Reference  2,  f22  is  derived  to 
be 


_ 2.78 _ 

<s2+2.668s+2.78) 


(4-17) 


Due  to  a  limitation  on  the  number  of  equations  that  the 
simulation  program  can  accommodate,  f22  is  simplified  to 


f 


2 

22  (s+2) 


(4-18) 


Equation  (4-18)  is  a  good  approximation  of  Eq.  (4-17)  in 
the  important  low  frequency  region,  and  satisfies  the  pre¬ 
filter  design  criteria  of  Reference  2.  The  compensator  g2 
is  not  guaranteed  to  be  without  some  overdesign,  but  a 
method  to  reduce  some  of  the  overdesign  has  been  developed. 
This  method  is  known  as  the  "improved  method,"  and  is 
applied  in  the  design  of  compensator  g^  in  Section  IV. 5  (8). 


IV.  5 

into  t^ 


Loop  Bound  Design, 

The  exact  equations  for  and  t2^  are  substituted 
and  t12  to  yield  the  following  equations 


'11 


fllLl? 

1+L, 


(4-19) 


where 

qiQi  i  (l+L- ) 

L!e  =  qlQlle  =  U-Y12+V 


(4-20a) 


Y 


(4- 20b) 


v 

Q 


llw22 

12Q21 


It  is  found  that  in  this  problem  y  <<  1  for  all  cases. 

Simplifying  Eq.  (4-20a)  to  1^  =  L1#  then  Eq.  (4-19) 

e 

,  _  £11L1 
11  L+L1 


becomes 

(4-21) 


The  y12  equation  becomes 


f  22L2 


'12  Q12{  (1+Lx)  (l+L2)-y} 


(4-22) 


For  such  that  1+L2  ~  L2,  Eq.  (4-22)  reduces  to 
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"12  Q12(l+L1) 


(4-23) 


Equations  (4-21)  and  (4-23)  are  used  to  derive  the  bounds 
on  L^q.  The  procedure  is  similar  to  that  used  for  L2q, 
except  hi  consists  only  of  the  term  due  to  the  command 
input.  The  loop  interaction  term  disappears  since  this  is 
a  BNIC  problem.  In  the  low  frequency  region  where  1 1+L^J  I 
|LjJ/  the  t12  expression  creates  the  greatest  demands  on 
Lio,  and  the  following  inequality  is  used  to  derive  the 
bounds  in  this  region. 
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t12  =  Q.-U+L.  )  -  b12 


(4-24) 


b12Q12 


i+Ll|  - 


(4-25) 


1+L1 | ~ | Lx 


,Llo'  -  bl2Ql2  . 

min 


(4-26) 


As  an  example,  at  u  =  0.1  r/s,  using  Eq.  (4-26)  L1q  is 


calculated  to  be 


L.  >  0  -(-45.5  +  1.34)  :  44  dB 
lo  — 


This  becomes  the  bound  for  at  u  =  0.1  r/s  since  the 
contour  of  (1+L1q)-1  is  essentially  a  horizontal  line  on 
the  MNC  for  this  value  of  L1q.  Bounds  in  the  medium  and 
high  frequency  regions  are  obtained  as  is  done  for  the  L2q 
bounds  in  the  same  regions.  The  loop  transmission  L^q  is 
shaped  next. 


The  nominal  q.^  selected  for  L^Q  is 


qllo 


_ -2.235s(s+.0107) (s+,4845) 

(s-1.167) (s+2.028) (s2+.012944s+. 006131) 


(4-27) 


Notice  that  Eq.  (4-27)  has  an  unstable  pole  at  s  =  1.167. 
To  avoid  trying  to  achieve  exact  cancellation,  L^Q  must  be 
shaped  with  this  unstable  pole  as  a  factor.  One  technique 
is  to  shape  to  include  the  pole  and  zero  combination, 

-(s-l  iIt|  *  Tllis  expression  has  a  magnitude  of  one,  and 
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a  phase  angle  that  decreases  with  increasing  frequency. 
Each  bound  on  L1q  as  derived  above  is  shifted  by  the 

arg  ,  and  L^"  is  shaped  using  the  new  bounds. 

The  desired  is  then  obtained  by  shifting  back  in 

the  direction  of  the  original  bounds  by  the  same  number 
of  degrees.  The  relationship  between  L^q  and  Q '  is 


,  .  (s+1.167) 

lo  "  ^lo  (s-1.167) 


(4-28) 


For  ease  of  loop  shaping,  this  method  is  found  to  be  pre¬ 
ferable  to  including  only  the  unstable  pole.  The  Nichol's 
Chart  plots  of  L1q  and  are  seen  in  Figure  4-10,  and  the 

Bode  plot  of  L1q  is  seen  in  Figure  4-11.  Figure  4-10  shows 
that  there  is  overdesign  in  L1q.  The  bandwidth  of  L1q  is 
found  to  be  24  r/s.  Dividing  L1q  by  Q11q,  g-^  is  found  to 
be 


*1 


-2.39(106)  (s+.  3)  (s+1.1)  (s-H.167)  (s+2.028)  (s+20)  (s+ 
s2  (s+. 0107)  (s+,4845)  (s+. 5)2 (s+35)  (s+ 


430) (s2+. 01 2944s+. 00613) 
87)  (s2+1272s+l ,123,600) 


Notice  that  (s-1.167)  does  not  appear  as  a  factor  in  the 
denominator  of  g^,  and  exact  cancellation  is  not  a  problem. 

A  simple  design  is  found  to  suffice  for  f^.  The 
expression  for  f^  is 


fn  =  2/  (s+2) 


(4-30) 


LOOP  TRANSMISSION  -  LIO 


Fig.  4-11.  Bode  Plot  of 


IV. 6  Stability  Design  for  Three 

Simultaneous  Failures 

In  this  section  the  characteristic  polynomials 
are  derived  for  two  configurations  with  three  of  the  four 
control  surfaces  failed.  For  this  study,  it  is  desired  to 
achieye  only  a  stable  system.  No  attempt  is  made  to  meet 
performance  specifications.  Bounds  for  stability  are  formed 
using  templates  that  include  all  six  CSC  modes  previously 
analyzed  plus  the  configuration  of  interest  for  which  only 
one  surface  is  controllable.  As  before,  a  nominal  condi¬ 
tion  is  chosen  and  the  stability  bounds  are  drawn  on  the 
Nichol's  Chart.  The  3  dB  contour  is  used  for  gain  and 
phase  margins.  The  nominal  loop  transmission  is  then 
shaped  as  is  done  for  L2q  in  Chapter  IV.  By  examining  the 
zeros  of  the  characteristic  polynomial,  which  is  a  function 


TABLE  5-2 


PITCH  RATE  RESPONSE  FIGURES  OF  MERIT 


2.16 

3.12 

2.20 

3.24 

3.84 

2.16 


compliance  with  the  performance  tolerances  listed  in 
Table  4-1,  repeated  here  as  Table  5-1.  Settling  time, 
rise  time,  and  peak  value  are  calculated  for  each  response 
and  listed  in  Table  5-2. 

TABLE  5-1 

FIGURE  OF  MERIT  TOLERANCES 
FOR  PITCH  RATE 

Settling  Time  Rise  Time  Peak  Value 

t  (sec)  t  (sec)  M  (deg/sec) 

S  P 

1.5  -  3.0  0.33  -  1.63  <  11 


From  Table  5-2  it  is  seen  that  design  tolerances 

for  t  are  not  satisfied  in  all  instances.  The  t  calcu- 
s  s 

lated  for  modes  2,  4,  and  5  of  flight  conditions  1  and  2 

is  found  to  exceed  the  maximum  allowed  t  of  3  seconds. 

s 

An  explanation  for  the  increased  settling  times  at  flight 
conditions  1  and  2  may  lie  in  the  relative  locations  of  the 
poles  and  zeros  of  the  basic  plant  equations  P.^3  and 
Sec.  II. 2.  The  poles  of  these  basic  plant  equations  for  q 
at  flight  conditions  3  and  4  tend  to  be  located  farther 
to  the  left  in  relation  to  their  zeros  than  are  corres¬ 


ponding  poles  of  the  same  plant  equations  for  flight  condi¬ 
tions  1  and  2.  Therefore,  system  closed-loop  roots  are 
probably  more  negative  at  flight  conditions  3  and  4  and  the 
resulting  settling  times  will  be  greater  at  flight  condi¬ 
tions  1  and  2.  The  t  is  easily  reduced  by  decreasing  the 


PITCH  RATE  (DEG/SEC) 

.00  4.00  8.00 


TIME  (SEC) 


Fig.  5-2  (c).  Pitch  Rate  Responses 
Flight  Condition  Three 


O 

O 


Fig.  5-2  (d) .  Pitch  Rate  Responses 
Flight  Condition  Four 


Fig.  5-2 (b) .  Pitch  Rate  Responses 
Flight  Condition  Two 


is  necessary  due  to  a  limitation  on  the  number  of  equations 
that  the  simulation  program  can  accommodate.  Figure  5-1 
shows  a  Bode  plot  comparison  of  Eqs.  (4-28)  and  (5-1). 
Reducing  the  order  to  g^  by  one-half  is  certain  to  make 
the  implementation  of  the  controller  more  practical. 


Fig.  5-1.  Bode  Plot  Comparison  of  and 
the  Approximation  for 

The  q  responses  for  the  four  flight  conditions 
and  six  CSC  modes  are  seen  in  Figures  5-2a  to  5-2d.  Eac 
figure  represents  the  responses  for  one  flight  condition 
and  the  six  CSC  modes.  The  responses  are  examined  for 


are  explained  for  compliance  with  the  performance  toler¬ 
ances  detailed  in  Chapter  IV.  Surface  deflections  are 
examined  for  saturation,  and  a  more  efficient  division  of 
control  effort  is  investigated  to  reduce  or  eliminate 
saturation.  Control  surface  initial  rates  are  found  to 
saturate  for  more  cases  when  the  input  is  a  pulse,  but 
representative  examples  with  ramped  inputs  are  also  pre¬ 
sented  showing  rates  that  are  not  saturated.  Mr.  Yin-Kuei 
Liao  of  the  Chung- Shan  Institute  of  Science  and  Technology, 
Republic  of  China,  developed  the  computer  program  used  to 
simulate  the  system  of  Figure  2-1.  It  is  recommended  that 
follow-on  thesis  efforts  in  QFT  concentrate  on  developing  a 
computer-aided  design  (CAD)  package  adaptable  to  higher- 
order  systems.  A  listing  of  the  computer  programs  used  in 
the  simulation  is  presented  in  Appendix  F.  Compliance  with 
the  specifications  for  q  are  investigated  first. 

V.2.1*  Pitch  Rate  Command  Response 

Before  proceeding  with  the  q  responses,  it  is 
noted  that  the  compensator  g^,  Eq.  (4-28)  ,  is  approximated 
as  follows: 

_  -6.6  (105)  (s+4)  (s+100)  (s2+0.013s+0.006)  ,r  1V 

g  =  - - -  (5-1) 

i  s(s+0.0223)  (s+25)  (s+630) 

Thus,  the  order  of  g^  is  reduced  from  tenth  to  fifth-order. 
Equation  (5-1)  is  obtained  from  a  straight-line  approxi¬ 
mation  of  the  Bode  plot  for  Eq.  (4-28).  The  approximation 


V.  Results 


V. 1  Introduction 

The  results  of  the  computer  simulation  are  pre- 
sented  in  this  chapter.  The  digital  simulation  program 
integrates  the  system  equations  in  state-variable  form  to 
yield  the  desired  outputs.  The  q  and  p  command  and  cross¬ 
coupling  responses  are  examined  for  compliance  with  per¬ 
formance  tolerances.  Control  surfaces  are  examined  for 
proper  sense  of  movement  and  saturation.  The  insights  that 
QFT  affords  a  designer  are  investigated,  and  the  results 
of  control  surface  "hardover"  simulations  are  illustrated. 

A  brief  description  of  the  simulation  program  follows. 

V. 2  Simulation 

The  system  represented  in  Figure  2-1  is  simulated 
with  the  plant,  prefilter,  and  compensator  equations  in 
state-variable  form.  A  digital  integration  routine  is 
applied  to  the  system  over  a  desired  time  interval  to  yield 
a  nonlinear,  real-time  simulation  which  includes  control 
surface  position  saturation.  The  outputs  of  the  simulation 
are  p,  q,  and  the  four  control  surface  positions  (right/ 
left  elevators  and  right/left  flaperons) .  Outputs  are 
obtained  for  each  design  flight  condition  and  control 
surface  configuration  (CSC)  mode,  and  the  q  and  p  responses 


transmissions  are  then  shaped  on  the  Nichol's  Chart.  Simple 
expressions  for  the  prefilters  are  derived  completing  the 
designs.  Lastly,  stability  only  designs  for  configurations 
with  three  simultaneous  surface  failures  are  synthesized, 
and  further  illustrates  the  potential  of  quantitative  feed¬ 
back  'theory.  Chapter  V  presents  the  system  simulations. 


(4-39) 


v-.v 


n+r  x  (s+1.56)  (s+32. 5ij46. 9)  (s+424. 6ij567) 
'  2o  '  "  s(s+22.9) (s+34.6) (s+429. 2± j 572 . 2) 

and 


ll+Llo 


) 

e 


(s+1.996)  (s+1388±j!838)  (s-t-1426±jl914) 
s  (s+1408±  jl877)^ 


(4-40) 


The  characteristic  polynomial  is  formed  by  the  product  of 
Eqs.  (4-39)  and  (4-40).  Again,  the  system  is  stable  as  all 
the  zeros  of  Eq.  (4-37)  lie  in  the  left-half-plane. 

These  designs  for  stability  further  illustrate  the 
power  of  QFT.  Loss  of  aircraft  control  is  the  probable 
result  if  time  is  required  to  identify  and  switch  in  com¬ 
pensation  for  the  two  failure  cases  examined  in  this  sec¬ 
tion.  No  identification  and  switching  in  of  compensation 
is  required  for  the  fixed  compensators  designed  using  QFT. 
This  buys  the  aircrew  valuable  time;  for  example,  time  to 
maneuver  the  aircraft  into  more  favorable  ejection  param¬ 
eters*. 


IV. 7  Summary 

This  chapter  presents  a  complete  description  of 
the  procedure  used  in  the  compensator  designs.  First,  the 
equivalent  plant  expressions  are  derived,  and  the  equiva¬ 
lent  single-loop  plant  equations  (Q^j)  are  obtained. 
Specifications  are  derived  for  the  command  and  cross¬ 
coupling  responses.  Based  on  these  specifications,  bounds 
on  the  loop  transmissions  are  obtained,  and  the  loop 
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60.00  -40.00  -20.00  0.00  20.00  40,00 


and 


(4-36) 


n  .  (s+5.7)  (s+28 . 7±j  37 . 3)  (s+278. 5tj 371 . 1) 

U loj  s  (s+22. 7)  (s+34.6)  (s+281 . 4±  j 375 . 2) 


The  characteristic  polynomial  is  the  product  of  Eqs.  (4-35) 
and  (4-36),  and  as  can  be  seen,  all  zeros  are  in  the  left- 
half-plane,  indicating  a  stable  system.  Figures  4-12  and 
4-13  show  L2o  and  ,  respectively,  and  the  bounds  on 


The  same  procedure  as  above  is  used  to  design  for 
stability  when  only  one  flaperon  is  controllable.  In  this 

case  only  P22a  an<*  P12a  ex;“-st  =  P2ia  =  ^l^  =  P21b 

=  P22  ~  P12  =  0) •  When  these  substitutions  are  made 
into  Eq.  (3-7)  ,  the  characteristic  polynomial  is 


(1  +  yP12  g^  +  P22  g2) 


(4-37) 


Equation  (4-37)  may  be  factored  as  follows 


<i+p22ag2) 


(l+upl2agi) 

1+P22S 


(1+L2o  ).  (l+Llo  ) 
e  e 


(4-38) 


where. 


22 


and  L 


lo_ 


tlP12agl 

1+P22a92 


The  factors  in  Eq.  (4-38)  are  derived  in  a  manner  similar 
to  that  used  to  derive  the  factors  in  Eq.  (4-34) .  The 
design  process  yields 
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yglP12 

yglP223 


g2P12a 
g2P22 — 


l  +  yP12  9i  +  p 


22 


(4-32) 


The  Characteristic  polynomial  for  this  configuration  is 

1  +  ypi2  gl  +  P22  g2  (4-33) 

Equation  (4-33)  is  factored  as  follows: 


a  (i+y pi2  gx) 

(1+P22  ’2'  ,  -  <1+L2o  »  <1+L 


1+p22a9; 


lo 


(4-34) 


Where, 


t  _  D  a  _  .  _  _  yP12apgl 

L2o  P22  og2 '  and  Llo  “  ,  p  a 
e  e  1+P22  og2 


Prefilters  f^  and  f2  in  Eq.  (3-7)  are  made  equal  to  unity 

for  the  convenience  of  this  analysis. 

Loop  transmissions  are  shaped  separately  on  the 

Nichol's  Chart  for  each  factor  in  Eq.  4-34.  The  L2o  loop 

e 

is  shaped  first  so  that  the  factor  (1+L_  )  can  be  used  in 

°e 

forming  the  bounds  for  the  loop.  From  Total,  the 

6 

factors  of  Eq.  (4-34)  are 


<1+L2o  > 


e 


(s+4.9)  (s+18. 4)  (s+299. 7fj399. 7) 
s  (s+22. 9)  (s+300± j 400) 


(4-35) 


47 


of  and  L 2Q /  the  stability  of  the  system  is  determined. 
Configurations  for  which  only  one  elevator,  and  only  one 
flaperon  are  controllable  are  examined  separately. 

The  first  configuration  analyzed  is  the  system  for 
which  only  one  elevator  surface  remains  controllable.  The 
remaining  elevator  and  the  two  flaperons  are  failed  to  their 
neutral  reference  positions.  The  design  begins  by  deter- 
mining  £e.  Only  and  ^22  ex;*-st  f°r  this  configuration 

(pnb  =  P21b  -  P22a  =  P12a  =  P22b  =  P12b  =  0).  These  sub- 
sti tut ions  are  made  into  Eq.  (3-6)  to  yield 
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— e 


P 
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a 

11 

uPn 

a 

21 

ppjl 

(4-31) 


As  before,  in  order  to  determine  the  equivalent 

single-loop  plant  equations  Q. .,  P  -1  must  be  computed. 

1  j  e 

From  Eq.  (4-31)  |P  |  =0,  and  P  does  not  exist.  In 
this  case  the  system  effectively  becomes  a  SISO  system 
with  only  one  output  independently  achievable.  The  Pj^g 
(i=j)  plant  equations  are  now  used  in  the  detailed  design 
procedure. 

When  the  elements  of  Eq.  (4-31)  are  substituted 
into  the  system  transfer  function  matrix  Eq.  (3-7)  ,  the 
system  matrix  for  one  controllable  elevator  becomes 
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time  constant  of  f 


Scheduling  of  f^  may  need  to  be 


11* 

employed  if  the  t  is  too  small  for  some  flight  conditions. 
Close  examination  of  the  data  in  Table  5-2  also  reveals 
that  the  difference  between  the  maximum  and  minimum  t 

s 

is  2.14  sec.  This  exceeds  the  t  tolerance  of  1.5  sec  by 

s 

0.64 'sec.  Quantitative  feedback  theory  guarantees  that 
the  original  performance  tolerances  will  be  satisfied 
despite  the  uncertainty  (2)  .  An  exact  reason  v?hy  this  has 
occurred  is  not  known,  but  one  possible  explanation  is  that 
the  approximation,  for  g^  may  not  be  adequate  to  satisfy 
the  tolerances  in  all  cases.  Despite  this  fact,  q  is  found 
to  be  controllable  for  each  flight  condition  and  CSC  mode. 

A  roll  angle  tolerance  when  q  is  commanded  must  also  be 
satisfied. 

As  stated  earlier  in  Section  IV. 3,  it  is  desired 
that  the  roll  angle  be  less  than  five  degrees  when  q  is 
commanded.  The  roll  angle  is  estimated  by  graphically 
calculating  the  area  under  the  p  response  curves  resulting 
from  cross-coupling.  The  cross-coupling  responses  are 
seen  in  Figures  5- 3a  to  5- 3d.  Only  four  responses  are  seen 
per  figure  as  cross-coupling  does  not  occur  for  CSC  modes 
one  and  ?ix  due  to  the  symmetry  of  the  control  surfaces 
for  these  two  modes.  The  estimated  worse  case  total  roll 
angle  at  each  flight  condition  is  seen  in  Table  5-3. 
Inspection  of  this  data  reveals  that  the  roll  angle  specifi 
cation  is  satisfied  for  all  cases.  Control  surface  deflec¬ 
tions  when  q  is  commanded  are  examined  next. 


TIME  (SEC) 


Fig.  5-3 (a).  Roll  Rate  Cross-Coupling 
Responses,  Flight  Condition  One 


Fig.  5-3  (b).  Roll  Rate  Cross-Coupling 
Responses,  Flight  Condition  Two 
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Fig.  5-3 (c).  Roll  Rate  Cross-Coupling 
Responses,  Flight  Condition  Three 


Fig.  5-4  (d) .  Roll  Rate  Cross-Coupling 
Responses,  Flight  Condition  Four 
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Fig.  5-4 (a) .  Elevator  and  Flaperon  Deflections 
for  Pitch  Rate,  CSC  Mode  One  and 
Flight  Condition  One 
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Fig.  5-4  (b) .  Elevator  and  Flaperon  Deflections 
for  Pitch  Rate,  CSC  Mode  Two  and 
Flight  Condition  One 


Fig.  5-4 (c) .  Elevator  and  Flaperon  Deflections 
for  Pitch  Rate,  CSC  Mode  Three  and 
Flight  Condition  One 
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Fig.  5-4 (d) .  Elevator  and  Flaperon  Deflections 
for  Pitch  Rate,  CSC  Mode  Four  and 
Flight  Condition  One 


are  taken  from  Reference  6,  and  repeated  in  Table  5-4. 
Inspection  of  Figures  5-4a  to  5-4 f  show  that  surface 
position  saturation  does  not  occur  for  the  q  command  at 
flight  condition  one.  Initial  surface  rates  are  seen  to 
be  quite  high,  and  are  estimated  to  saturate  or  nearly 
saturate  for  the  pulse  input.  Saturation  is  further 
addressed  in  a  later  section  of  this  chapter.  In  summary, 
control  surface  movements  appear  to  be  in  the  proper  sense 
to  generate  positive  q  for  all  flight  conditions  and 
CSC  modes. 


TABLE  5-4 

SURFACE  POSITION  AND  RATE  LIMITS 

Position 

Rate  Limit 

Surface 

Limit  (deg) 

(deg/sec) 

Elevators 

±25 

60 

Flaperons 

±20 

52 

V.2.2  Roll  Rate  Command  Response 

The  p  responses  for  the  four  flight  conditions 
and  CSC  modes  are  seen  in  Figures  5-5a  to  5-5d.  Each 
response  is  examined  for  compliance  with  the  tolerances 
listed  in  Table  4-2,  repeated  here  as  Table  5-5. 
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Fig.  5-5  (c).  Roll  Rate  Responses 
Flight  Condition  Three 


Fig.  5-5 (d) .  Roll  Rate  Responses 
Flight  Condition  Four 


TABLE  5-5 


FIGURE  OF  MERIT  TOLERANCES 

FOR  ROLL  RATE 

Settling  Time 

Rise  Time 

Peak  Value 

tg  (sec) 

tr  (sec) 

Mp  (deg/sec) 

3.0  -  6.0 

0.65  -  3.23 

<  55 

The  figures  of  merit  calculated  for  each  p  response  are 
listed  in  Table  5-6.  As  is  necessary  for  g^,  the  com¬ 
pensator  g2,  Eq.  (4-16),  is  approximated  in  the  same  manner 
as  g^.  The  approximation  for  g2  is 

_  -3.68  (104)  (s+3)  (s+900)  (s2+2s+4)  /r  ^ 

g2  s(s+0.25)  (s+50)  (s+158)  (s+3750) 

The  original  ninth-order  equation  for  g2  is  reduced  to  an 
effective  fifth-order  equation.  Figure  5-6  shows  a  Bode 
plot  comparison  of  Eqs.  (4-16)  and  (5-2).  The  data  in 
Table  5-6  is  examined  for  compliance  with  performance 
tolerances. 

The  data  in  Table  5-6  indicates  that  performance 

tolerances  for  t  ,  t  ,  and  M  are  satisfied  for  all  but, 

s  r  p 

possibly,  one  case.  The  t  for  CSC  mode  6  and  flight 

S 

condition  one  is  greater  than  five  seconds  due  to  the 
oscillatory  nature  of  the  response.  Also,  most  of  the 
responses  have  settling  times  that  are  faster  than  desired. 
These  responses  may  be  slowed  by  using  a  larger  time  constant 


TABLE  5-6 


ROLL 

RATE  RESPONSE  FIGURES  OF  MERIT 

S 

Flight 

CSC 

Settling  Time  Rise  Time 

Peak  Value 

1 

Condition 

• 

Mode 

t  (sec)  t  (sec) 

b  i 

M  (deg/sec) 

1 

4.25  1.15 

50.7 

- 

2 

4.16  1.22 

51.1 

1 

3 

3.84  1.17 

51.0 

X 

4 

3.84  1.30 

51.9 

5 

4.00  1.42 

52.4 

6 

5.00+  1.14 

53.0 

1 

1.85  1.06 

50.1 

2 

4.16  1.21 

51.1 

3 

1.80  1.10 

50.0 

A 

4 

1.76  1.10 

50.0 

5 

1.70  1.02 

49.9 

6 

1.50  1.00 

50.1 

1 

1.92  1.10 

50.0 

— 

2 

1.88  1.09 

50.0 

'•  • 

o 

3 

1.86  1.10 

50.0 

X 

4 

1.76  1.16 

50.0 

5 

1.84  1.10 

50.0 

6 

1.76  1.03 

50.0 

• 

1 

1.92  1.00 

50.0 

- .  ■■■ 

2 

1.98  1.04 

50.1 

A 

3 

1.88  1.07 

50.0 

•* 

4 

4.48  0.95 

51.1 

■-Tv 

5 

1.94  1.05 

50.1 

.* 

6 

1.76  1.06 

50.1 

**  .  -I 

•  *•  *> 

-  .--I 
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Fig.  5-6.  Bode  Plot  Comparison  of  G2  and 
the  Approximation  for  G2 

in  the  prefilter  f22*  A  higher-order  prefilter  can  be 
used  to  shape  the  responses  to  obtain  more  desirable 
figures  or  merit.  This  is  not  attempted  in  this  study  due 
to  a  restriction  on  the  number  of  equations  that  the  simu¬ 
lation  program  can  accommodate.  The  q  resulting  from  the 
p  command  is  examined  for  compliance  with  performance 
specifications. 

The  pitch  angle  due  to  cross-coupling  is  estimated 


in  the  same  way  that  roll  angle  due  to  cross-coupling  is 
estimated.  It  is  desired  that  the  pitch  angle  remain 
less  than  three  degrees  during  the  p  command.  The  q 


cross-coupling  responses  are  shown  in  Figures  5-7a  to  5-7(3, 
and  Table  5-7  lists  the  worse  case  pitch  angle  estimated 
at  each  flight  condition. 


• 

TABLE  5-7 

ESTIMATED  PITCH  ANGLE  WITH 
ROLL  RATE  COMMANDED 

Flight 

CSC 

Pitch  Angle 

Condition 

Mode 

(deg) 

1 

5 

-3.10 

2 

4 

-0.70 

3 

4 

-1.00 

4 

4 

-1.25 

The  magnitude  of  the  pitch  angle  for  CSC  mode  5 
and  flight  condition  1  is  seen  to  slightly  exceed  3  degrees. 
If  the  accuracy  of  the  pitch  angle  estimations  is  con¬ 
servatively  estimated  at  5  percent,  the  performance  spe¬ 
cification  for  pitch  angle  is  satisfied.  Control  surface 
deflections  due  to  the  p  command  are  estimated  for  proper 
sense  of  movement  and  saturation. 
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Fig.  5-7 (a).  Pitch  Rate  Cross-Coupling 
Responses,  Flight  Condition  One 


Fig.  5-7  (b) .  Pitch  Rate  Cross-Coupling 
Responses,  Flight  Condition  Two 
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Fig.  5-7  (c) .  Pitch  Rate  Cross-Coupling 
Responses,  Flight  Condition  Three 
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Fig.  5-7 (d) .  Pitch  Rate  Cross-Coupling 
Responses,  Flight  Condition  Four 


Elevator  and  flaperon  deflections  for  the  p  command 
at  flight  condition  one  are  seen  in  Figures  5-8a  to  5-8f. 

As  is  seen  in  Figure  5-8a,  the  surface  pairs  are  moving 
differentially  in  the  proper  sense  to  generate  positive 
roll  rate.  The  flaperons  are  also  commanded  in  excess  of 
their  20  degree  limit,  while  the  elevators  are  well  below 
their  limit  of  25  degrees.  This  is  essentially  the  case 
for  the  remaining  figures  in  this  sequence.  It  is  obvious 
from  Figures  5-8  that  more  of  the  control  effort  needs  to 
be  shared  by  the  elevators  to  reduce  flaperon  movement. 

This  point  is  addressed  later  when  saturation  is  discussed. 
Initial  surface  rates  for  all  surfaces  are  near  or  at 
saturation  for  the  pulse  input;  these  rates  are  greatly 
reduced  when  the  input  is  ramped,  as  is  seen  later.  The 
surface  deflection  figures  for  p  at  flight  conditions  two, 
three,  and  four  are  displayed  in  Appendix  G. 

It  is  concluded  from  the  pitch  rate  and  roll  rate 
results  that  the  compensators  g^  and  g2  provide  robust 
control  for  all  flight  conditions  and  CSC  modes  for  which 
the  system  is  designed  despite  minor  qualifications  for  the 
pitch  rate  results.  It  is  also  reasonable  to  conclude  that 
the  system  provides  robust  control  for  all  points  encom¬ 
passed  by  the  design  flight  conditions.  Saturation  is 
the  next  topic  of  discussion. 


Fig.  5-8 (e) .  Elevator  and  Flaperon  Deflections 
for  Roll  Rate,  CSC  Mode  Five  and 
Flight  Condition  One 
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Fig.  5-8 (f ) .  Elevator  and  Flaperon  Deflections 
for  Roll  Rate,  CSC  Mode  Six  and 
Flight  Condition  One 


V. 3  Saturation 

Different  results  are  obtained  for  several  of  the 
cases  when  the  provision  for  position  saturation  is 
included  in  the  simulation  program.  Flight  condition  one, 
CSC  modes  one  and  two  are  taken  as  cases  in  point.  The 
surface  deflections  when  p  is  commanded  for  CSC  mode  one 
without  saturation  are  presented  in  Figure  5- 8a,  and  it  is 
evident  that  the  flaperons  show  deflections  in  excess  of 
their  limit.  The  same  simulation  is  run  with  position 
saturation  included,  and  the  results  are  seen  in  Figure 
5-9.  The  flaperons  are  seen  to  saturate,  and  the  elevator 
deflections  increase  to  generate  the  commanded  roll  rate. 
The  p  responses  with  and  without  saturation  are  seen  in 
Figure  5-10.  Note  the  more  oscillatory  nature  of  the 
response  with  saturation.  The  result  is  much  more  dramatic 
for  CSC  mode  two;  Figure  5-11  displays  the  responses  for 
this  CSC  mode  with  and  without  saturation.  Cross- 
coupling  is  found  to  increase  significantly  when  surfaces 
saturate,  as  is  seen  in  Figure  5-12a.  Surface  deflections 
corresponding  to  Figure  5-12a  are  seen  in  Figure  5-12b. 

It  is  obvious  from  these  illustrations  that  position 
saturation  reduces  the  stability  of  the  system.  Satura¬ 
tion  can  be  reduced  or  eliminated  by  better  distribution 
of  the  control  effort. 
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Fig.  5-9.  Elevator  and  Flaperon  Deflections 
With  Saturation,  CSC  Mode  One 
and  Flight  Condition  One 
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Fig.  5-10.  Roll  Rate  Response  With  and 
Without  Saturation,  CSC  Mode  One 
and  Flight  Condition  One 
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Fig.  5-11.  Roll  Rate  Response  With  and 
Without  Saturation,  CSC  Mode  Two 
and  Flight  Condition  One 
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Fig-  5-12 (a) .  Pitch  Rate  Cross-coupling  Response 
With  and  Without  Saturation,  CSC  Mode  Two 
and  Flight  Condition  One 


Fig.  5-12  (b) .  Elevator  and  Flaperon  Deflections 
for  CSC  Mode  Two  and  Flight  Condition  One 

The  division  of  control  effort  between  the  ele¬ 
vators  and  flaperons  may  be  varied  by  changing  the  con¬ 
stant's  y^2  and  y2^  in  Figure  2-1.  As  discussed  in  Sec¬ 
tion  III. 2,  these  constants  are  arbitrarily  set  at  0.25. 
Increasing  y12  causes  the  elevators  to  share  more  of  the 
effort  in  generating  p.  This  is  seen  in  Figure  5-13  where 
y^2  is  increased  to  0.5.  By  comparing  Figure  5-13  with 
Figure  5-9,  it  is  seen  that  the  flaperons  no  longer 
saturate,  as  more  of  the  command  is  being  generated  by  the 
elevators.  For  a  reconfiguration  strategy  to  be  effec¬ 
tive,  some  sort  of  scheme  to  optimize  the  division  of  con¬ 
trol  surface  effort  is  essential.  It  is  recommended  that 


Fig.  5-13.  Elevator  and  Flaperon  Deflections 
for  Roll  Rate  with  V12  ~  0*5,  CSC  Mode  One 
and  Flight  Condition  One 


follow-on  thesis  efforts  in  flight  control  system  recon¬ 
figuration  investigate  strategies  to  optimize  the  division 
of  control  effort.  Appendix  B  contains  theory  developed 
by  Dr.  Horowitz  to  optimize  the  division  of  control  effort. 
Initial  surface  rates  are  noted  earlier  to  be  high. 

So  far,  all  inputs  have  been  pulse  commands. 

Figures  5-14  and  5-15  illustrate  the  effect  on  surface 

rates  when  the  inputs  are  ramped.  In  Figure  5-14  the  input 

2 

is  a  clipped  ramp  with  a  slope  of  20  deg/sec  (0.5  sec  to 
peak  value) .  The  elevator  and  flaperon  initial  rates  are 


found  to  be  well  below  saturation  levels.  The  same  is 


Fig.  5-14.  Elevator  and  Flaperon  Deflections  for 
Pitch  Rate  with  a  Ramped  Input,  0.5  sec  to 
Peak  Value  (CSC  Mode  One  and 
Flight  Condition  One) 
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Fig.  5-15.  Elevator  and  Flaperon  Deflections  for 
Roll  Rate  with  a  Ramped  Input,  1  sec  to 
Peak  Value  (CSC  Mode  One  and 
Flight  Condition  One) 
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TABLE  A- 2 


AIRCRAFT  DATA  FOR  FLIGHT  CONDITION  TWO 
'  (Mach  =  0.6  and  30,000  feet) 


q 

2 

(dynamic  pressure--lbs/f t  )  =  158. 

91 

VT 

(trim 

velocity--f t /sec) 

=  596.91 

(trim 

angle-of- attack — 

deg)  =  4.705 

Longitudinal  Body  Axis  Primed  Dimensional  Derivatives 

Z  " 

=  -0.526422 

M  "  =  2.52708 

X  " 

23.0402 

a 

a 

a 

h  ' 

=  -0.066156 

M6  "  =  -5.86214 

X6  ' 

= 

3.17035 

e 

e 

e 

Z,  " 
6f 

=  -0.111711 

M,  "  =  -0.211773 
5f 

X*  " 
6f 

= 

-2.09855 

z  ' 
q 

=  0.997184 

M  "  =  -0.341902 

q 

X  " 

q 

= 

-48.8785 

Z  " 
u 

=  -0.000109 

Mu"  =  -0.000337 

Xu' 

= 

-0.005142 

V 

=  -0.004425 

MQ"  =  0.000313 

xe' 

= 

-32.0915 

Lateral-Directional 

Body  Axis  Primed 

Dimensional 

Derivative; 

V 

.=  2.29583 

L'  =  -19.2246 

p 

Y  " 

e 

= 

-0.154099 

N  " 

P 

=  -0.000888 

L  "  =  -0.893601 

P 

.  Y  " 
p 

= 

0.082387 

N  " 
r 

=  -0.278676 

Lr'  =  0.318845 

Y  ' 
r 

= 

-0.998322 

N<S 

=  -1.96651 

L6  "  =  3.92325 

Y  " 

= 

0.021165 

r 

r 

r 

N6 

a 

=  -0.268403 

L6  "  =  -17.4468 
a 

Y  ' 

a 

r= 

0.000357 

N, 

"  =  -1. 50547 

hx  '  =  -13.5832 

Y 

1  & 

= 

0.014398 

6dt 

6DT 

6dt 

N(5  " 

=  1.51008 

L6  '  =  0.414519 

Y  " 

c 

= 

0.007335 

c 

°c 

TABLE  A-l 


AIRCRAFT  DATA  FOR  FLIGHT  CONDITION  ONE 
(Mach  =  0.2  and  30  feet) 


q 

VT 

•  a 

T 

2 

(dynamic  pressure--lbs/f t  )  = 

(trim  velocity — ft/sec)  =  223 

(trim  angle-of-attack — deg)  = 

59.18 

.  26 

14.98 

Longitudinal  Body  Axis  Primed  Dimensional  Derivatives 

V 

=  -0.459802 

M^'  ■  0.542375 

X  "  = 

a 

12.5375 

V 

e 

=  -0.074045 

Mfi  "  =  -2.23634 
e 

X6  '  = 
e 

0.357830 

Z,  ' 
5f 

=  -0.061526 

Mr  -  =  0.241866 

st 

X.  "  = 
6f 

-1.97404 

Z  ' 

q 

=  0.995737 

M  "  =  -0.633651 

q 

X  "  = 

q 

-58.0974 

z  ' 

u 

=  -0.000786 

Mu"  =  -0.000651 

v  - 

0.002886 

V 

=  -0.037268 

Mq'  =  -0.000058 

V  - 

-31.1064 

Lateral-Directional 

Body  Axis  Primed  Dimensiona  1  Derivatives 

N3 

=  0.919806 

=  -11.0532 

Y  '  = 

8 

-0.143895 

N  " 

P 

-0.003864 

Ip'  =  -0.934188 

Y  "  = 

P 

0.261703 

N  " 
r 

=  -0.257624 

L  "*  =  1.03044 

Y  "  = 

r 

-0.994485 

N6  ' 
r 

=  -0.814818 

L6  "  =  1.42280 

r 

Y  '  = 

6 

r 

0.021662 

NS  ' 
a 

=  0.075009 

L6  .  =  -4.47725 

a  , 

V 

0. 010422 

\r 

"=  -0.603086 

L,  -4.28449 

V 

Y  "= 

DT 

0.022018 

V 

=  0.328110 

Lj  -  =  -0.296006 
c 

Y  '  = 

X 

0.006234 

Appendix  A:  Flight  Parameters  and 

Aerodynamic  Derivatives 

Flight  parameters  and  aerodynamic  derivatives  are 
presented  in  this  appendix  for  the  four  design  flight  con¬ 
ditions.  These  data  are  taken  directly  from  Reference  6, 
and  repeated  here  for  the  convenience  of  the  reader.  An 
extensive  aerodynamic  package  was  developed  by  the  General 
Dynamics  Corporation  for  the  AFTI/F-16.  Derivatives  are 
formulated  in  the  body  axis  system,  and  have  the  units  of 
radians,  radians/sec,  or  ft/sec;  i.e.  ,  XQ"  (radians),  X  " 
(radians/sec),  Xu"  (ft/sec),  etc. 
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VI. 3  Recommendations 

The  following  recommendations  follow  directly  from 
the  preceding  conclusions: 

1.  Research  using  quantitative  feedback  theory 
for  the  design  of  multivariable  flight  control  laws  should 
continue,  particularly,  for  reconf igurable  flight  control 
systems.  Future  investigations  should  be  for  systems  of 
third  or  higher-order.  Optimizing  the  division  of  control 
effort  and  digital  control  are  fundamental  prerequisites 
for  a  practical  reconf igurable  flight  control  system,  and 
should  be  included  in  future  research  efforts.  Also,  a 
more  accurate  aircraft  model  should  be  developed  to  include 
coupling  in  the  plant  matrix. 

2.  A  computer-aided  design  package  should  be 
developed  for  quantitative  feedback  theory.  This  package 
should  include,  but  not  be  limited  to  the  following  capa¬ 
bilities:  Bode  plot,  Nichol's  Chart,  plant  template  formu¬ 
lation,  loop  transmission  boundary  derivation,  loop  trans¬ 
mission  shaping,  curve  fitting,  and  system  simulation. 

In  addition  to  allowing  more  difficult  design  problems  to 
be  solved,  these  capabilities  would  greatly  speed  up  the 
design  process  and  result  in  more  precise  control  laws. 


systems  may  be  significant.  A  fixed-compensator  design 
has  the  potential  to  significantly  minimize  the  complexity 
of  fault  identification  schemes  that  may  be  required  for  a 
reconf igurable  flight  control  system.  The  next  section 
presents  some  of  the  conclusions  that  may  be  drawn  from 
this' thesis . 

VI. 2  Conclusions 

1.  Procedures  developed  from  quantitative  feed¬ 
back  theory  are  effective  in  the  design  of  multivariable 
control  laws  for  normal  or  impaired  flight  control  systems. 
The  technique  reveals  important  insights  to  the  designer 
about  compensation  requirements.  This  fixed-compensator 
design  approach  may  significantly  reduce  the  anticipated 
complexity  of  reconf igurable  flight  control  systems. 

2.  The  application  of  quantitative  feedback  theory 
in  the  design  of  multivariable  control  laws  is  a  reason¬ 
ably  'straightforward  process.  The  technique  is  relatively 
simple  once  the  basic  theory  is  understood,,  and  the 
designer  has  some  practice  in  its  use. 

3.  The  controllers  developed  in  this  thesis  pro¬ 
vide  robust  control  for  all  flight  control  surface  con¬ 
figurations  and  flight  conditions  for  which  the  design  is 
synthesized. 

4.  An  interactive  computer-aided  design  package 
should  be  developed  for  quantitative  feedback  theory  so 
that  design  problems  of  larger  scope  can  be  attempted. 


VI.  Conclusions  and  Recommendations 

VI. 1  Discussion 

The  primary  purpose  of  this  thesis  is  to  demon¬ 
strate  the  effectiveness  of  quantitative  feedback  theory 
as  applied  to  flight  control  system  reconfiguration.  Com¬ 
pensators  are  synthesized  to  be  robust  for  an  unimpaired 
flight  control  system  as  well  as  for  several  different 
cases  in  which  the  flight  control  system  is  severely 
impaired.  The  compensators  are  also  to  be  robust  for  a 
large  portion  of  the  aircraft  flight  envelope.  Quantita¬ 
tive  feedback  theory  is  applied  to  the  linearized,  small 
perturbation  equations  of  motion  for  a  longitudinally 
unstable  aircraft  (AFTI/F-16).  The  application  of  quanti¬ 
tative  feedback  theory  is  found  to  be  relatively  straight- 
forward  for  this  difficult  multiple  input-multiple  output 
problem.  The  transparency  of  the  technique  provides  impor¬ 
tant  insights  into  the  compensator  requirements  from  the 
initial  stages  of  the  design  process.  The  results  of  this 
thesis  demonstrate  that  quantitative  feedback  theory  is 
well  suited  for  flight  control  system  reconfiguration  design. 
The  fixed-compensator  designs  are  seen  to  provide  satis¬ 
factory  control  at  all  design  points.  The  importance  of 
fixed-compensator  designs  in  reconf igurable  flight  control 


compensator.  Control  surface  "hardover"  is  also  investi 
gated,  and  simulation  results  for  this  condition  further 
illustrate  the  utility  of  quantitative  feedback  theory. 
Recommendations  and  conclusions  are  presented  in  the 
following  chapter. 


TIME  ( SEC  ) 


Fig.  5-19  (b).  Roll  Rate  Response  from  Flaperon 
Hardover,  Flight  Condition  Two 

V.  6  Summary 

This  chapter  presents  a  description  of  the  simula¬ 
tion  program  and  the  results  for  the  pitch  rate  and  roll 
rate  commands.  From  an  examination  of  the  command  and 
cross-coupling  responses  it  is  determined  that  the  com¬ 
pensators  g^  and  g2  provide  reasonably  robust  control  over 
the  design  flight  conditions  and  control  surface  configura¬ 
tion  modes.  Saturation  is  examined,  and  the  necessity  to 
optimize  the  division  of  control  effort  is  evident.  Quanti¬ 
tative  feedback  theory  is  seen  to  provide  the  designer  with 
valuable  insights  prior  to  the  actual  synthesis  of  the 


tihe  (sec) 


Fig.  5-18  (b) .  Roll  Rate  Response  from  Elevator 
Hardover,  Flight  Condition  Two 


Fig.  5-19  (a).  Pitch  Rate  Response  from  Flaperon 
Hardover,  Flight  Condition  Two 


example,  let  6^  fail  to  the  hardover  position,  the  result 

output  vector  (P..a5.a  •  ^l^l*  ^  is  treate<*  as  a 

max  max 

disturbance  to  be  attenuated.  The  equivalent  plant  trans¬ 
fer  function  matrix  P  for  this  condition  then  becomes 

— e 

identical  to  P  for  6.a  =  0  (CSC  mode  2)  since  6.a  is  no 
— e  1  i 

longdr  controllable.  Mr.  Yin-Kuei  Liao  has  simulated 
separate  elevator  and  flaperon  hardovers  using  the  com¬ 
pensators  derived  in  Chapter  IV.  Simulation  results  are 
seen  in  Figures  5-18  to  5-19.  The  system  is  seen  to  effec 
tively  attenuate  the  hardover  surface's  outputs.  The 
results  of  this  section  further  demonstrate  the  utility  of 
QFT  for  the  challenge  of  FCS  reconfiguration.  See  Appen¬ 
dix  B  for  Dr.  Horowitz's  explanation  of  the  hardover  sur¬ 
face  problem. 


Fig.  5-18  (a).  Pitch  Rate  Response  from  Elevator 
Hardover,  Flight  Condition  Two 


Figure  5-16.  The  additional  area  of  plant  uncertainty 
is  enclosed  by  the  dashed  line.  Flight  conditions  one 
through  four  are  identified,  respectively,  by  the  letters 
a,  b,  c,  and  d.  The  difference  between  the  bound  on  L 
at  a)  =1.0  r/s  with  and  without  the  plant  uncertainty  due 
to  both  elevators  failed  is  seen  in  Figure  5-17.  The 
bounds  in  Figure  5-17  represent  the  demands  on  L  due  to 
the  command  response  t^.  Bound  A  is  developed  from  the  Q-^ 
template  with  the  additional  uncertainty  of  both  elevators 
failed,  while  bound  B  is  developed  from  the  same  template 
without  the  extra  uncertainty.  Clearly,  it  is  seen  from 
Figure  5-17  that  the  larger  the  amount  of  plant  uncertainty 
the  higher  is  the  bound  on  L,  and  higher  bounds  mean  a 
larger  loop  transmission  bandwidth.  The  L  bandwidth  is 
increased  by  more  than  one  and  one-half  octaves  from  bound 
B  to  bound  A. 

This  example  further  illustrates  the  utility  of 
QFT  as  a  design  technique  for  FCS  reconfiguration.  This 
method  provides  the  designer  with  valuable  insights  very 
early  in  the  design  process  as  to  the  "cost"  of  feedback. 

V.5  Control  Surface  "Hardover" 

A  condition  known  as  "hardover"  occurs  when  a  sur¬ 
face  without  command  suddenly  deflects  to  its  maximum 
position  and  remains  there.  The  output  from  this  control 
surface  may  be  treated  as  an  external  disturbance  vector 
which  is  to  be  attenuated  by  the  feedback  system.  As  an 


2 

true  when  the  p  input  is  ramped  with  a  slope  of  50  deg/sec 
(1  sec  to  peak  value) ;  initial  rates  are  also  well  below 
saturation  levels  as  seen  in  Figure  5-15.  Figures  5-14 
and  5-15  may  be  compared  with  Figures  5-4a  and  5-8a  where 
the  inputs  are  pulsed  to  see  the  marked  reduction  in  the 
surface  rates. 

V. 4  QFT  Insights 

Quantitative  feedback  theory  provides  the  designer 
with  insight  at  the  initial  design  stages  into  the  nature 
of  the  loop  transmission.  By  examining  the  extent  of  plant 
uncertainty,  and  the  narrowness  of  the  performance  specifi¬ 
cations,  the  designer  gains  valuable  insight  about  the 
bounds  on  L  and,  hence,  on  L  itself.  For  instance,  the 
larger  the  area  of  plant  uncertainty  and/or  the  tighter 
the  performance  specifications,  the  higher  are  the  bounds 
on  L.  Thus,  the  designer  has  a  feel  early  in  the  design 
for  the  compensation  required.  If  he  feels  that  the  demands 
on  L  are  too  high  he  may  choose  to  reduce  the  amount  of 
plant  uncertainty  in  the  design  or  relax  the  performance 
tolerances.  For  this  thesis,  it  was  initially  planned 
that  the  compensator  designs  would  be  robust  when  both 
elevators  are  failed.  This  was  decided  against  after 
examining  the  template  of  and  seeing  the  extra  "cost" 
in  terms  of  loop  transmission  bandwidth  for  including  this 
configuration.  The  additional  plant  uncertainty  in  the 
pitch  channel  when  both  elevators  are  failed  is  seen  in 


TABLE  A- 4 


AIRCRAFT  DATA  FOR  FLIGHT  CONDITION  FOUR 
(Mach  =  1.6  and  30,000  feet) 

q  (dynamic  pressure — lbs/ft2)  =  1129.31 
VT  (trim  velocity — ft/sec)  =  1591.75 
oirp  (trim  angle-of -attack — deg)  =  1.68 


Longitudinal  Body  Axis  Primed  Dimensional  Derivatives 
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This  appendix  contains  quantitative  inherent 
reconfiguration  theory  developed  by  Dr.  Isaac  Horowitz. 
Reconfiguration  synthesis  theory  is  developed  for  2  input 
2  output,  5  input-3  output,  and  10  input-6  output  systems 
utilizing  control  of  individual  flight  control  surfaces. 
Theory  to  optimize  the  division  of  control  surface  effort 
under  no  failure  and  failure  conditions  is  also  presented 


A  QUANTITATIVE  INHERENT  RECONFIGURATION  THEORY 


FOR  A  CLASS  OF  SYSTEMS 

Isaac  Horowitz* 

Summary 

In  aircraft  flight  control,  most  control  surfaces  are  in 
pairs  (elevators,  ailerons,  canards  etc),  with  each  pair  normally 
controlled  as  a  single  unit.  If  a  surface  fails,  the  usual 
approach  is  to  attempt  explicit  identification  and  switch-in  of 
compensation  prepared  for  that  contingency.  In  this  paper  each 
surface  is  separately  controlled,  permitting  "inherent  recon¬ 
figuration",  wherein  the  design  is  apriori  made  such  that  despite 
one  or  several  simultaneous  surface  failures,  the  system  still 
satisfies  the  original  performance  tolerance  (of  course  over  a 
smaller  dynamic  range) ,  with  the  same  original  fixed  compensa¬ 
tion. 


Inherent  reconfiguration  is  a  natural  extension  of  "Quan¬ 
titative  Feedback  Theory". (QFT) ,  wherein  the  system  design  is 
tuned  to  the  plant  uncertainty  set  F  =  {P},'  and  to  the  accept¬ 
able  system  output  set  A.  In  QFT  one  apriori  designs  so  that 
the  system  output  is  in  A  for  all  P  in  P.  Surface  failures 
simply  enlarge  the  set  P.  The  transparency  of  QFT  enables  the 
designer  to  readily  see  the  extra  "cost  of  feedback"  for  this 
enlargement  of  P  by  inclusion  of  surface  failures. 
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A  QUANTITATIVE  INHERENT  RECONFIGURATION  THEORY 
FOR  A  CLASS  OF  SYSTEMS 

Isaac  Horowitz 

1.  INTRODUCTION 

Consider  a  system  with  uncertain  constrained  multiple 
input  -  multiple  output  (MIMO)  part  denoted  as  the  plant  P  with 
m  control  inputs  and  n  outputs  to  be  controlled.  If  m>n  then 
n-n  degrees  of  freedom  are  available  for  optimization  purposes. 
In  an  "output-controllable"  plant  (Horowitz  and  Shakei  1975) , 
the  n  desired  outputs  uniquely  determine  the  remaining  n  free 
inputs.  If  n>m,  only  m  outputs  are  independently  achievable. 
Hence,  an  effective  n*n  plant  is  normally  assumed  in  uncertain 
MIMO  systems,  if  the  problem  is  to  achieve  desired  input-output 
relations  despite  uncertainty  (Horowitz  1979,  1983) . 

There  is  an  important  plant  class  wherein  one  or  more  of 
the  n  plant  inputs  is  actually  a  pair  e.g.  elevator,  canard  and 
aileron  pairs  in  an  aircraft.  The  elevator  pair  is  usually 
counted  as  one  control  input,  because  its  two  parts  are  tied 
or  driven  together  to  achieve  longitudinal  mode  control,  and 
therefore  do  not  affect  the  lateral  variables  such  as  roll  and 
yaw.  Similarly  the  two  ailerons  are  usually  tied  or  driven 
differentially  for  lateral  control,  and  therefore  this  pair  aces 
not  affect  longitudinal  variables  such  as  pitch  and  normal 
acceleration.  However,  the  individual  elements  of  each  of  these 
pairs  do  affect  both  the  longitudinal  and  lateral  variables  and 
this  fact  is  indeed  used  for  "reconfiguration"  purposes  in  case 
of  significant  control  surface  damage.  Such  reconfiguration 


requires  explicit  identification  of  the  failure  and  the  switch- 
in  of  standby  compensation  elements  designed  for  this  specific 


purpose.  The  challenge  (Rubertus  1983,  Chandler  1934,  Chandler 
and  Potts  1983)  treated  in  this  paper 

is  to  devise  a  design  theory  in  which  no  explicit  identifica¬ 
tion  is  required.  Rather,  even  under  severe  control  surface 
failures,  the  tolerances  are  to  be  satisfied  automatically,  with 
no  need  for  switching  in  of  new  compensations.  The  same  compen¬ 
sation  functions  are  to  be  adequate  under  normal  (no  failure) 
and  under  the  maximum  possible  control  surface  failures.  This 
is  denoted  here  as  "inherent  reconfiguration”.  The  question 
of  how  optimum  such  a  single  design  is  for  the  variety  of  failure 
conditions,  is  discussed  in  Section  5.  Three  configurations  of 
increasing  complexity  are  treated  in  Sections  2-4. 

1 . 1  Synthesis  objectives 

A  linear  time-invariant  plant  model  is  assumed.  (See 
Section  5  for  discussion  of  nonlinear  plants) .  There  are 
specified  tolerances  on. the  closed-loop  system  responses  to 
commands  r^  in  Figure  lr  as  follows.  Let  y  =  Tr,  T  =  [t^^] , 
and  each  t^ has  its  set  of  acceptable  functions  T^_.  either  in 
the  time  or  frequency  domain.  It  is  particularly  convenient 
to  assume  "basically  noninteracting"  (BNIA)  tolerances  of  the 
form  It.j (ju) |  <_  b^j(w)  for  i^j  (see  Section  5  for  discussion). 
There  is  also  a  set  of  possible  external  disturbance  vectors 
V  -  {cf}  acting  on  the  plant  whose  effect  on  the  system  output 
is  given  by  v^  =  T^d,  =  J'  tolerances  of  the  form 

1 1  j  |  <_  b^  (w)  ,  giving  sets  of  acceptable  disturbance  response 

‘  i  j  1  i  j 

functions  T  ,  .  These  tolerances  on  the  elements  of  T  and  T  . 


are  to  be  achieved  over  the  entire  set  of  plant  variation 
(different  flight  conditions  in  flight  control)  and  over  the 
maximum  possible  simultaneous  control  input  failures,  by  means 
of  a  fixed  set  of  compensation  functions. 


2.  A  TWO-PAIR  STRUCTURE 
2 . 1  Introduction 

The  simplest  plant  for  illustrating  the  theory  is  shown 

(in  darker  lines),  together  with  its  compensations,  in  Figure  1, 

ci  b 

For  example,  the  (6^,6“)  pair  could  be  the  right  and  left 
3  Id 

elevators  (S2,<52)  the  right  and  left  ailerons,  y^  the  pitch 

angle,  y2  the  roll  angle.  The  effect  of  each  individual  j 

surface  on  each  i  output  is  given  explicitly  by  P^j'P^j  with 

pa .  =  p*?.  under  normal  (no  failure)  conditions.  If  6~,6, 
d  Id 

(or  6^,6^)  have  equal  values,  i.e.  the  right  and  left  ailerons 
(or  elevators)  are  deflected  equally  in  the  same  direction. 


then  y2(roll)  =  r(6ap2i  -  S^p^)  =  That  is  why  minus  signs 


,b_b 


are  associated  with  P^2'^22  F^9ure  i •  and  plus  signs  for  all 
a  b  ;  • 

Pli'Pli  ' 

If  the  command  inputs  (r^,r2)  and  the  plant  outputs 
(yi,y2) *  hut  no  internal  plant  variables  can  be  measured 
(Horowitz  1982a),  then  the  most  general  feedback  structure  has 
16  free  functions:  8 f ^ ^  from  the  two  r ^  to  the  four  6^  and 
8itk  j  to  the  four  6^  from  the  two  y^  in  Figure  1.  But  it  is 
assumed  here  that  in  the  normal  (no  failure)  case,  zero 
interaction  between  the  lateral  and  longitudinal  outputs  is 


desired  i.e.  t^  =  0  for  i  /  j,  also  that  pa^  =  p  —  . 


This 


accounts  for  the  various  symmetries  in  the  f  ^  ^ of  Figure  1. 
Thus,  if  r^  -  0,  r2  7s  0  the  desired  y^  =  0 ,  y2  ^  0  are  achievabl 
by  5a  =  -6°,  62  =  ~Z>2'  so  opposite  signals  via  (m^2,-m^2), 
(m22,-m22)  '  ^22 '""^22^  are  nee^e^-  However,  the  same 

(m11,m11)  (m21,m21)  and  f^  give  6a  =  6°  if  y1  jt  0,  y2  =  0  are 

wanted  in  response  to  r^  /  0,  r2  =  0.  Hence,  only  8  free  com¬ 
pensating  functions  are  left.  Note  that  Figure  1  is  only  one 
of  an  infinitude  of  possible  canonic  structures  (Horowitz  1963) 
exploiting  the  available  freedom.  Simple  equivalence  relations 
permit  any  one  to  be  replaced  by  another  structure  with  no  need 
for  redesign. 

2 . 2  Development  of  Synthesis  Theory 

In  Figure  1,  y  =  Tr  =  P  =  P(Fr-My) ,  so 

T  =  [ t±  j  ]  =  (I+PM)  "*1PF  with  (1) 
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Quantitative  synthesis  theories  for  satisfying  each  of  the 


(3c) 


n  t^j  tolerances  of  (1)  despite  the  uncertainty,  have  been 


developed  (Horowitz  1979,  1982b)  and  applied  (Horowitz  and 
Loecher  1981,  Horowitz  et  al  1981,  19£l,  Betzold  et  al  1984)  to 
significant  examples.  These  techniques  have  the  merit  of 
reducing  the  n*n  MI MO  problem  into  a  number  of  single  input  - 
single  output  (SISO)  single-loop  problems,  whose  solutions  are 
guaranteed  to  solve  the  original  MIMO  problem.  It  would  be 
very  desirable  to  apply  this  theory  to  the  present  problem 
class  (see  also  Section  2.2.2).  For  this  purpose  it  is 


necessary  to  identify  a  n*n  equivalent  plant  P  and  compensator 


M  , F  matrices  (n=2  in  this  section) ,  such  that 
e  e  . 


T  =  (I+PM)~1PF  Of  £q{l)  =  (I+P  M  )  XP  F 

6  6  6  6 


(4) 


Also  for  the  BN I A  tolerances  of  Section  1.1,  and  F^  should 

e  e 


be  diagonal.  Examination  of  (3a-c)  reveals  this  is  achieved 
by  setting 
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The  available  quantitative  design  theory  is  directly 

aoolicable  to  the  new  equivalent  P  ,M _, F  system,  including 
-  6  6  6 

surface  failures,  as  is  next  seen.  Note  in  (6a)  that 

p.  .  (i^i)  =0  in  the  no  failure  case.  See  Section  2.2.2  for 

additional  strong  justification  of  relations  (5a, b). 

2.2.1  Failure  cases 

Various  kinds  of  failure  are  possible.  Total  failure  of 
is  defined  as  =  0  for  all  i,  which  may  involve  new 

values  of  the  other  j »  if  the  aerodynamic  coefficients  in¬ 
volved  in  these  p^j  are  affected  by  the  failure.  The  design 

theory  easily  accommodates  this  with  the  new  value  P^je(i/j) 

b  b 

no  longer  zero,  but  p12e  =  -Pi2pll'  p21e  =  ~p21  and 

Pile  "  Pll+“21p12'  P22e!='P22+l,12Pb2r  Partial  failures  like- 

wise  result  in  new  values  of  the  P^je»  without  changing  the 

form  of  the  MIMO  problem.  The  quantitative  feedback  theory 

(QFT)  of  (Horowitz  1979,  1982b)  was  devised  precisely  for  such 

situations  of  the  p. .  assuming  different  values  due  to  un- 

6 

certainty.  The  various  failure  cases  simply  give  additional 

elements  in  the  set  P  =  { [ p . .  ] }  of  possible  plant  matrices. 

6  1 3  6 

Suppose  a  fixed  compensation  design  is  being  made  for  m^ 
representative  flight  conditions  involving  different  Mach  and 
altitude  values  and  m2  different  failure  cases.  Then  at  most 


n^m2  (instead  of  m^)  sets  pf  matrix  element  values  constitute 
the  plant  uncertainty  set  P  ,  to  which  the  design  theory  is 
applied  in  precisely  the  same  manner  as  if  there  were  only  m^ 
sets.  it  is  even  conceivable  that  the  extent  of  the  new  un¬ 
certainty  (measured  by  the  ranges  of  the  individual  frequency 
res'ponses  {p.  .  (jo)}  at  any  o>)  is  not  much  enlarged  by  the 

i  j  e 

various  failure  cases. 

"Hardover"  is  a  type  of -failure  in  which  one  or  more  control 
surfaces  are  frozen  at  some  fixed  values  and  can  no  longer  be 
changed.  If  6^  is  so  frozen  at  (say)  its  maximum  value  , 

then  let  its  resulting  vector  of  output  components 
{pll5lx  '  p2151x) ’  =  de  be  treated  as  an  external  disturbance 
vector,  which  is  to  be  suitably  attenuated  by  the  feedback 
system,  in  order  to  satisfy  the  system  performance  tolerances. 

The  effective  plant  upon  which  d  acts  is  P  of  (6a)  in  which 

S  G 

S3  3 

p 2^  =  =  0/  because  6^  is  no  longer  controllable.  The 

synthesis  theory  can  handle  simultaneous  disturbance  attenuation 
and ’command  input  response  under  plant  uncertainty. 

•  *  4 

2.2.2.  Additional  justification  of  constraints  (5a, b) 

One  might  argue  that  (5,ab)  may  not  be  optimum  relations 
inasmuch  as  they  are  apparently  required  only  so  that  the 
specific  synthesis  techniques  of  Horowitz  (1979,  1982b)  may  be 
directly  used.  It  is  next  shown  that  they  are  in  fact 
necessary  in  order  that  failures  may  lead  to  small  and  con¬ 
trollable  channel  interaction  despite  plant  uncertainty.  Thus, 

consider  the  system  transfer  matrix  T  under  6^  failure,  fcr  which 

a  a  .  .  . 

P11  =  P?1  =  9lving 


*  b  lb  b 

(f llPll+f21?12) (1+m22?22)  |  P11P22 {m12f 22~m22f 125 _f12Pll 

+  f21P12rn12^21_  '  _ 

p12P21(mllf2l“fllm21)_f21P21  [  (f12P21+f22P22)  (1+m2lP12) 

»  ,  b 

1  +  f22P22mllPll 


b  b 

(X+m2^p^2)  (l+rn2.2P21+m22P22^  +  mllPll  ^+m22P22^  '  T^e  dominant  terms 
in  the  numerator  of  any  t  involve  m  p  p  products,  so  small 
tuv  (u^v)  relative  to  t  require  such  products  disappear  in  the 
numerators  of  t  independently  of  the  actual  p  values.  This  is 
achieved  in  (7)  only  by  (£a,b) ,  because  the  other  alternative  of 

setting  for  example  PiiP22^m  12^22-m22^12^  =  ^12pll  imPract^ca^- 
when  there  is  non-negligible  plant  uncertainty.  It  is  readily 
found  that  (5a,b)  suffice  for  this  purpose  for  6^  ,  6°  ,  <S2  •  & 2 
failures  one  at  a  time,  or  any  two  at  a  time. 

It  is  arguable  that  the  same  result  is  achieved  in  (7)  by 
setting  m . .  =  f. .  =  0  for  all  i/ j ,  which  is  a  special  case  of 

1  J  1  J 

(5a,b) .  However  it  would  then  be  impossible  to  achieve  t22  under 
simultaneous  failure  of  $2  an<3  ^2  '  or  to  ac^ieve  t^l  ^  both 
6^  and  f2  failed.  Also’'  62  ,  62  would  not  contribute  to  y^  when 
r^O  and  6^  ,  6^  would  not  contribute  to  y2  when  r^O.  Thus 
(5a, b)  are  inherently  necessary  in  any  synthesis  technique  which 
achieves  the  objectives  in  this  paper. 


2.2.3  Division  of  control  effort 

(a)  No  failure  case .  In  order  to  find  the  6’s  needed  to  achieve 
the  desired  outputs  note  that  6  =  Fr-My  =  (F-MT)r.  Take  any  r^O 
the  other  r^=0  so  Fr  =  f1  the  ith  column  of  F,  and  Tr  =  t1  the  ith 
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column  of  T.  In  the  no  failure  case  if  i=l ,  and  if  =  p  — 

giving  t2^  =  0,  the  result  is  <5^  =  <5^  =  anc^ 

62  =  62  =  f21-m21tll'  SO  that  due  to  (5a)  '  =  fll/,f21  = 

l/u21(s).  Similarly  for  r2^0,  r1=0  and  p^V  =  ^ ,  in  the  no 

failure  case  <5^  =  -6^  =  f]_2~mi2t22'  ^2  =  ~^2  =  ^22-m22t2':>  ' 
so  '(5b)  gives  6^/ <52  =  ^12^22  =  ^12^'  T^e  designer  maY  choose 
v21,lJ12  to  ac^aeve  the  optimum  division  of  control  effort  between 
<5^  and  6^.  Note  that  in  the  traditional  older  design  procedure 
where  the  lateral  and  longitudinal  modes  are  decoupled,  there  is 


no  such  division  of  control  effort  —  only  the  6^  are  used  to 
achieve  y^. 

The  following  appears  a  reasonable  way  to  select  the  u  — • 
Under  command  t^  =  p^6^  +  P^2“2"  Let  ^2  =  ^l^l^ll^^  ' 

<{>21  a  constant,  so  that  t^  =  Pil^l  d+^2i^  ,  order  that 

ci  Id 

51  '  5i  work  together  to  achieve  t^ ;  <}>21  can  possibly  be  chosen 

cl  cl 

on  the  basis  of  the  relative  saturation  levels  of  6-^  ,  62-  This 

gives  (see  2.2.3)  u2^  =  ^21P11/P12“  Similarly,  consideration  of 

r2  command  leads  to  <5^  =  52(*>12p22//p21  (r50'ri=0)  =  yi262  '  0f 

* 

course,  the  p^^  values  change  for  various  flight  conditions,  but 
the  ratios  p^/p^#  P22//p21  c^ani?e  muc^  less.  A  compromise 
fixed  Pj^/P^j  ratio  can  be  chosen.  Also,  if  scheduling  is  used, 
then  the  in  j  values  may  also  be  scheduled. 

(b)  Control  effort  division  under  failure  The  above  equation 
<5  =  (f1-Mt1)r^  must  be  properly  used  under  failure,  because  it 
refers  to  the  commanded  5  values .  Hence,  if  for  example  6^  fails, 
the  above  gives  the  commanded  but  not  actual  value.  Simple 
equations  are  obtained  for  the  6's,  but  the  difficulty  is  in  the 


t.  .  values.  If  cf  fails,  then  t,..,  t_.  are  no  lonaer  zero  so 
1  j  1  1221 

for  r~=l/  r^=0 


f 22“ (m21t12+m22t22)  (f 22_m22t22) “y21mllt12 


(f22_n22t22)  "  y21mllt12 


y2(f22  m22t22*  “  mllt12 


It  is  readily  found  (Horowitz  1982b) ,  that 

t .  .p.  . 

4.  -  Jrjog 

t’-(lr])  p.  (1+t. )  ' 


f . . -m . . t . . 
n  li  li 


n,utii<1+VT) 

rriran 


where  T  =  P12e!?21e/E’llep22e  '  la  "  muuquu  ' 

• 


The  above  are  actually  general  results,  giving 


.b  -  l 

°1  L 


*a  .b 
52  62 


lp^  +  Ul2(1+£ 

•>  •  1+VY 


rY)„ 


which  is  **h^2  an  tbe  no  case,  because  then  P12e  =  °* 

6^  fails,  p^  =  p^  =  0  so  (10)  gives,  recalling  (6a), 

”^12 11 2 1^1 2^  ^ll+y21^12^  ;''n  tbe  control-frequency  range  where 
|  ( jtu)  |  >>l-y.  For  r^=l,  tbe  result  is  in  general 


62+62  "  U21 


1  +  y12£2P21e 

P22e(1  +  VY> 


which  is  1/y,.  in  the  no  fail  case  where  p0.  =  0,  and 
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the  inputs  to  the  systems  of  Figures  1-3.  Suppose  that  the 
desired  response  function  set  for  mode  x  is  x .  H  is  designed 
so  that  its  outputs  are  T.  c.  in  response  to  command  input 
number  x  of  value  c  .  Then  the  r.  of  the  systems  in  this  paper 

»  X  X 

have  the  values  r’.  =  T.  c.  .  H  has  as  amny  input  points  as  the 

X  X  X  X  X 

0 

modes  F  being  considered,  and  new  ones  can  be  added  without 
disturbing  the  designs  of  F,  M  of  this  paper.  The  latter  point 
is  the  principal  advantage  of  this  approach,  but  the  optimum 
under  BNIA  design  may  not  be  optimum  for  each  or  any  of  the 
various  modes.  Clearly  compromises  between  these  two  approaches 
are  feasible  with  the  design  theory  presented  here. 

It  is  also  worth  noting  that  the  well-known  "scheduling" 
technique  is  obviously  easily  accommodated  into  the  synthesis 
procedure.  An  advantage  of  quantitative  synthesis  is  that  the 
designer  clearly  sees  the  relations  between  "cost  of  feedback" 
(Horowitz  et  al  1982)  and  the  range  of  uncertainty  and  the 

tolerances,  so  the  trade-off  between  scheduling  complexity  and 

* 

cost  of  feedback  reduction  is  highly  transparent. 


condition  for  obtaining  the  equivalent  linear  set  is  to  solve 
the  nonlinear  plant  equations  backwards  -  given  the  plant  output 
to  find  the  plant  inputs.  This  means  the  parameters  equivalent 


to  the  of  this  paper  must  first  be  chosen  i.e.  one  must 
first  decide  how  to  use  the  freedom  due  to  the  plant  having 
more  inputs  than  outputs. 


BNIA  design 

Optimization  with  respect  to  the  has  been  previously 

discussed.  The  assumption  of  desired  "basically  noninter¬ 
acting  closed-loop  tolerances"  is  related  to  the  optimization 
problem,  as  follows.  In  advanced  fighter  aircraft  there  are 
many  different  operating  modes.  In  mode  i,  a  specific  com¬ 
bination  of  outputs  may  be  wanted  in  response  to  command 
(a  specific  set  of  t^  for  j  =  1  to  n)  .  In  mode  m  an  entirely 
different  set  of  tjm  in  response  to  rm,  may  be  specified  with 
its  set  of  tolerances.  It  is  conceivable  that  the  optimum  use 
of  the  available  freedom  (more  plant  inputs  than  outputs) ,  may 

f 

be  significantly  different  for  the  different  modes.  If  this 

A 

is  sufficiently  important,  different  designs,  with  conceivably 
different  compensations  inside  the  feedback  loops,  may  be  made, 
and  the  new  systems  switched  in  accordance  with  the  mode  desired 
The  design  theory  presented  here  is  clearly  readily  applied 
to  this  alternative,  as  the  design  technique  is  not  exclusively 
for  BNIA  systems  (Horowitz  and  Loecher  1931) . 

This  paper  offers  an  alternative  approach  to  the  many-mode 
design  problem.  A  prefilter  H  is  added  whose  output  provides 


based  upon  control  effort  division,  have  been  discussed  in  the 
paper.  Another  way  of  using  the  available  design  freedom,  is 
to  minimize  the  Pg  uncertainty.  Thus  in  (6a),  Plle  =  Pn+u21P12 
with  correlations  often  between  tne  parameters  of  p^  and  P-^* 

One  could  consider  choosing  so  as  to  minimize  the  uncertainty 

range  of  p^  <3ue  to  the  various  flight  conditions  and  failures. 
Other  optimality  criterion  could  also  be  considered,  e.g.  drag 
minimization. 

Sections  2-4  presented  the  synthesis  theory  by  application 
to  structures  whose  complexity  ranged  from  (4  plant  inputs  -  2 
outputs)  to  (10  plant  inputs  -  6  outputs)  .  In  all  of  these, 
two  "groups"  (longitudinal  and  lateral  for  aircraft)  were 
assumed.  But  it  is  obvious  that  there  can  be  as  many  groups 
as  outputs.  Recall  also  (Section  2.1)  that  different  compensa¬ 
tion  structures  may  be  used,  in  particular  G  =  H  in  the  forward 
path,  and  a  new  F'  prefilter  in  place  of  F  with  GF'  =  F  and 
the  outputs  of  -P  and  F'  being  the  inputs  into  G. 

.  Nonlinear  uncertain  plants 

Linear  time- invariant  plant  models  were  used  in  this  paper, 
permitting  use  of  transform  equations.  The  synthesis  theory  is 
readily  extended  to  uncertain  nonlinear  MIMO  systems  by  means 
of  the  "linear  plant  set  equivalence  technique"  (Horowitz  1976, 
1981a, b,c,  1932c,  Horowitz  et  al  1980a, b)  •  In  this  technique  the 
uncertain  nonlinear  plant  set  is  replaced  by  linear  time  in¬ 
variant  equivalent  sets.  The  solutions  (compensations)  for 
the  latter  problem  are  guaranteed  to  work  for  the  original 
nonlinear  problem  for  the  set  of  command  and  disturbance  inputs 
for  which  the  equivalence  has  been  developed.  A  necessary 


5.  DISCUSSION 


This  paper  has  presented  a  unified  synthesis  procedure  for 
flight  control  design  which  apriori  takes  into  account  in  a 
precise,  quantitative  manner  many  possible  combinations  of 
control  surface  failures,  the  various  flight  conditions  and 
external  disturbances  (gusts  etc.).  This  can  all  be  done  by 
means  of  a  fixed  compensat ion  design ,  providing  the  equivalent 
plant  matrix  Pe  satisfies  the  usual  necessary  conditions 
(Horowitz  1982b) .  The  minimum-phase  condition  (det  Pq  has 
no  right  half-plane  zeros)  guarantees  that  any  arbitrarily 
narrow  tolerances  can  be  satisfied.  If  Pg  is  nonminimum-phase 
it  is  necessary  that  the  performance  tolerances  are  compatible 
with  the  constrained  loop  bandwidths  achievable. 

It  is  conceivable  that  a  single  fixed  compensation  design 
to  handle  such  a  large  variety  of  conditions,  may  require  im- 
practically  large  loop  bandwidths.  If  so,  one  could  consider 
scheduling  for  different  classes  of  mach  number  and  altitude, 
in  order  to  reduce  the  range  of  uncertainty.  Scheduling  is 
quite  reliable,  whereas  'failure  identification  may  sometimes 
be  fairly  unreliable. 

Optimization 

Since  each  control  surface  is  assumed  separately  control- 
able,  the  number  of  plant  control  inputs  exceeds  the  outputs 
to  be  controlled  in  the  no  failure  condition.  The  exploita¬ 
tion  of  this  freedom  constitutes  an  optimization  problem,  which 
must  be  solved  before  the  design  for  uncertainty  and  failure 
can  proceed.  The  difficulties  involved  when  optimization  is 


m33t31  P31e6/P33e 


«  =  f11-n»11t11- 


From  the  matrix  eciuation 


<5  =  (F-MT)r,  for  r^  =  1  and  other  r^  =  0  ,  there  result  five 


scalar  equations  <5^  =  ^i2m22t21+'Ji3m33t31^  '  yii  =  *  for 

i  =  1  to  5.  The  m22t21  m33fc31  are  rePlaced  bY  P21e4//p22e  ' 


of  above,  giving  four  equations  for  j  =  2  to  5 : 


Xji  -  («/«,) 


J.UjlP22eP33e~  (yj2P21eP33e+tJ2  3P2  2eP31e) 


p22ep3  3e_(li12p21ep33e+y13p22ep31e) 


(20) 


The  Puve  are  of  course  functions  of  the  p^j  as  per  (18,19),  so 

the  above  equation  relates  the  p^j  to  the  p^  and  Xj^.  Similarly 

four  equations  relationg  X.j  (j=l,3,4,5)  to  the  p^j  and  p^^  are 

obtained  by  setting  only  =  the  others  zero,  etc.  By  this 

means  (r^  =  1^  all  others  2ero  for  i  =  1  to  6) ,  a  total  of  24 

equations  relating  the  to  the  24  p^j  are  obtained. 

4 . 2  Number  of  tolerable  failures 

At  least  six  independent  control  inputs  are  needed  to 

achieve  linear  control  of  the  six  outputs  in  Figure  3,  so  any 

four- of  the  ten  control  surfaces  may  fail  simultaneously.  The 

* 

design  can  be  apriori  made  so  that  the  desired  performance 

tolerances  on  all  36  t. .  of  T  are  satisfied  for  any  such  four 

13 

simultaneous  failures.  Five  failures  will  permit  only  twenty- 
five  t^j  tolerances  satisfied,  and  again  this  can  be  incorporated 
into  the  design,  etc.  -  all  this  to  be  achieved  by  a  single 


fixed  set  of  compensators. 


As  in  Sections  2,  3  interaction  between  the  two  sets  (1-3  and 
4-6)  is  zero  in  the  no  fail  case.  But  as  in  Sections  2,  3,  the 
paired  6^^  of  set  1  (61»62  here)  participate  in  achieving  set  2 
outputs  in  response  to  set  2(r2~rg)  commands,  and  vice  versa. 

The  design  procedure  is  the  same  as  in  Sections  2,  3.  Once  the 
y.T  are  chosen,  the  P  set  is  available  by  listing  all  flicht 
conditons  and  surface  failures  for  which  the  design  is  to  apply, 
and  proceeding  precisely  as  in  any  quantitative  MIMO  synthesis 
problem. 

4.1  Choice  of  y^j 

Essentially  the  choice  of  is  an  optimization  problem  — 

how  to  use  the  available  freedom  (more  plant  inputs  than  outputs) 

for  each  r.  input.  There  are  30  each  m.  f.  .  to  be  chosen: 

Setting  m^/m^  *  f^/fjj  of  (20b)  gives  24  relations,  the  actual 

y^j  values  gives  24  more  and  finally  quantitative  design  leads 

to  six  each  of  diagonal  ,  f^.  However,  quantitative  design 

can  proceed  only  after  the  y. .  of  (20b)  have  been  chosen,  per- 

«  — 

mitting  P  of  (18,19)  to  be  evaluated  for  the  various  flight 

C  • 

conditions  and  surface  failure  cases.  Suppose  the  X^j  = 

(6^/<5j)r  (all  ri  =  0  for  i^j)  have  been  chosen,  if  optimum 

division  of  control  effort  under  no  fail  is  the  optimality 
criterion.  Then,  as  in  Section  3.1,  the  X^  must  be  related  to 
the  y^j  and  a  procedure  similar  to  that  in  Section  3.1  can  be 
followed:  e.g.  for  r^  =  1  (all  other  r^  -  0) ,  only  t^  ,  t21  , 

t^  exist.  From  (I+PeMe)T  =  PgFe  and  neglecting  t^  or  t2^  cf 
t^  (justifiable  in  BNIA  systems)  and  |  ^  f  >  > 1  (highly  justifiable 
in  control  frequency  range)  ,  one  gets  ^22^21  ~  P21e<?/,p22e  ' 
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The  =  0  for  i  =  4-6  in  order  that  6^  =  0  for  r4  *  r3  or  rg 

commands  as  it  is  assumed  that  t^  =  0  is  desired.  Similarly  • 

cigj  =  0  for  j  =  1-3  so  6g  =  0  for  or  commands,  as 

t, .  =  0  is  assumed  desired. 

6d 

P  is  derived  and  the  (j .  .  defined  in  the  same  manner  as  in 
e  13 

Sections  2.2,  3  for  the  reasons  given  in  2.2,  2.2.2  giving 

T  =  (I+P  M  ) _1P  F  ,  and 
e  q  g  c 

First  three 


columns  of  P  = 
e 


P11  P12  P13  P14  P15 

P21  P22  P23  P24  P25 

P31  P32  P33  P34  P35 


AP4 1  Ap42  ° 
Ap51  Ap52  0 


Ap44  Ap4 5 
Ap54  AP55 


Ap61  Ap62  0  AP64  Ap65 


y  12  y  13 
1  v-n 


y31  y32  1 

y41  y42  y43 
y 51  y52  y53 


Last  three 

columns  of  Pq  =  Tp^  Ap12  Ap14  Ap15  0 

Ap21  Ap22  AP24  AP25  0 

Ap31  A]p32  Ap34  Ap3 5  0 


^42  H4  4  ^45  4  6 

P52  P54  P55  P56 

P62  P64  P65  P66 


y 14  y 15  y 16 
y 24  y 25  y 26 
1  y4 5  y 46 

Uc«  1  Pec 


y 64  y65  1 


-  u  •  • 

Api j  =  Pi j  "  pxj  (^3,6),  u.j  =  ^ 

33 

Pij  =  Pij  +  Pij  (^3,6). 

Me  =  diag.mii  ,  Fg  =  diag.f.. 


f  .  . 


(20a-e) 


VJ/.V.V. 


situations  can  be  included  in  the  design. 

4.  TEN  INPUT  -  SIX  OUTPUT  PLANT 

The  final  structure  treated  here  (Figure  3)  has  one  single 

and  two  pairs  of  control  surfaces  for  one  set  of  variables 

(e.g.  longitudinal)  and  a  similar  number  for  a  second  set  (e.g. 

lateral).  The  resulting  P,F,M  matrices  are 

Fa  ba  b  a  ba  b„~-l 

11  P11  P12  P12  P13  P14  P14  P15  P15  0 


a 

b 

’21 

P21 

P22 

a 

b 

a 

’31 

P31 

P32 

a 

_b 

a 

41 

_P41 

P42 

a 

b 

a 

51 

-p51 

P52 

a 

b 

a 

61 

•p61 

P62 

b 

*12 

P13 

p?4 

p14 

a 

p15 

p15 

b 

*22 

P23 

P24 

P?4 

a 

p25 

4 

b 

*32 

P33 

P34 

P34 

P35 

p35 

b 

'42 

0 

„a 

P44 

-p44 

P?5 

-p45 

& 

0 

P?4 

-pS4 

P55 

■52 

0 

P64 

~p64 

p65 

"P65 

It  is  assumed  6^  affects  only  the  second  set  of  outputs  and  6^ 
only  the  first  set,  hence  p; c  =  0  for  i  =  1-3,  =  0  for 

j 

F 


done  as  follows.  As  noted  in  2.2.3,  in  general  in  Figures  1,2 


6  =  (F-MT)r,  giving  for  r^  command  in  Figure  2: 


61  =  fll"(mlltll+m12t21) ' 

62  =  f2l“  (m21fcll+ra22t21)  =  y21 (f  ll-"1!!*!^  “  m22fc21 

6 3  =  f3l"  (m31tll+m32t21)  =  y31  (f  ll_mllfcll)  ~v'32m22t12  ' 


Equations  (9)  are  used  to  eliminate  ,  f^-m^t^. 

The  result,  after  writing  the  P. .  in  terms  of  p  and  y  ,  of 

*  13 e  ruv  Huv 

(13) ,  is 

y21  ~^21 

X.  (i-nPvi+Pni  )  +  X.p^+X.L.P,,  -  -  =  - 

1  21c22  C21  2r23  3  21r23  m 22  m22 

y31  _X3l 

xlp22X31  “X2P22+X3  (X31P23+P21)  “  = 

X1  =  1“p12u21  '  X2  =  y31~y21U32  '  x3  =  y32_y12y: 


xlp22X31 


.  P21  _ 

_X21 

m22 

m22 

)  "  — = 
m22 

~X  31 

m22 

X1  1-y12y; 


32  ^12^31 


Two  more  equations  are  obtained  for  r2^0,  r^=0,  setting  &±/&2  ~ 
X12  '  ^3^2  =  X32  9ivin9  a  total  of  four  equations  in  the  four 
y.  .  unknowns.  Unfortunately  the  as  yet  unknown  irs,  appears  in 
the  "first  two  and  m^  ir;  the  last  two,  but  the  terms  containing 
these  m^  appear  small  compared  to  the  others,  so  the  result 
should  be  relatively  insensitive  to  the  values. 

3 . 2  Number  of  tolerable  failures 

It  is  readily  seen  that  there  exists  a  linear  range  of 
control  of  all  three  outputs  in  Figure  2,  even  if  any  two 
(of  the  five)  surfaces  totally  fail  simultaneously.  But 
failure  of  any  three  leaves  only  two  outputs  or  any  two 
independent  combinations,  independently  controllable,  while 
failure  of  four  leaves  only  one  output  controllable.  All  these 


inherently  necessary  for  it  to  be  possible  to  achieve  under 
failure  longitudinal  output  with  small  lateral  effect,  and  vice 


versa . 


3.1  Division  of  control  effort  and  choice  of 


One  is  tempted  to  follow  the  procedure  in  Section  2.2.3,  if 
optimum  control  effort  division  is  desired  i.e.  for  r^l,  r2  = 

r3  =  °  in  tlx=  +  P12,S2  +  p1363  '  let  52/5l  =  <{,21?11/?12  ' 

6^/6^  =  <f3]_P1i/pl3  with  ^21' ^31  constants*  in  or{3er  that  the 
three  <?s  work  together  to  achieve  t^.  However,  the  resulting 
t21  =  Pn^i  ^1+<^21+(^31^  t^!en  not  independently  realizable. 


so  the  above  cannot  be  followed  for  r^,r 2  commands.  It  can  be 


applied  to  the  r^  command  because  t^  =  1 23  =  fc33  ~ 

P31lSl  +  p3353*  Let  6ly/|53  =  *13P33^P31*  Als°  61  =  f13~m13t33 

P13  ^f33-m33t33^  '  63  =  f33~m33t33  so  51^53  =  y13  which  9lv®s 

=  1i,i3P33//P3i  once  the  optimum  choice  of  6^/6^  has  been  made  . 

But  for  x±’r2  commands  the  optimal  division  of  the  control  effort 
is  more  complicated  as  it  is  subject  to  the  constaints  of  t21,t^2 

9 

satisfying  their  BNIA  tolerances  of  the  form  1 123_  ( j  (w) )  /t^  ( ju)  | 

£  e21(w),  [t12,t22(  -  e12^w^  with  e2l'e12  Jcnown*  Thus  if 

tll>/t21  =  Iplil5i//Ip2i6i  =  1^t21'  with  e21  very  this  9ives 

the  constraint  p21  +  P22*21  +  p23^31  ~  0  where  *n  =  (|5j/<51)r  ^0; 

for  r2  command  the  resulting  constraint  due  to  t^2  is 

12  +  P13X32  _  w'  "i2 


Pll'l2 


+  Pi  o  +  Pi  0,  -  (6^/6-,) 


1'  ~2'r2/0  * 


Suppose  the  A^.  values  have  been  "optimally"  chosen,  it  it 


then  necessary  to  realte  them  to  the  of  (14-c)  in  order  for 


P  of  (13)  to  be  known  and  the  design  to  proceed.  This  can  be 


a  b  a 

P11  Pll  P12  p13 

P^3 

“11 

a12  a13 

p  = 

a  b  a 

P21  P21  P22  P23 

P23 

;  m,f  = 

°11 

a12  ~a13 

dL  _3  n  3 

P31  _P31  0  P33 

^3 

a21 

a22  0 

a31  a32  a33 


a31  a32  ~a33 


(12a, b) 


a  =  m  for  M,  a  =  f  for  F. 

Note  that  p32  =  0  as  the  rudder  has  no  effect  on  the  longitudinal 

output;  (*23  =  0  as  62  =  0  is  desired  for  r^  command  giving  y3 

ci  b 

output  with  desired  y.^  =  =  0;  as  before  p^j  =  p°j  is  assumed. 

As  in  Section  2.2,  one  demands  T  =  (I+PeMe)  "'’PeFe  = 

(I+PM)-1PF  with  3x3  P  and  diagonal  M  ,F  giving  for  the  analogs 

S  G  ® 

of  (5,6) 


P  = 
e 


APij 


P11  P12  P13 

1 

P21  P22  P23 

y21 

Ap3l  0  Ap3^ 

»31 

-Apn  Ap13- 

~U13~ 

AP21  AP23 

1 

— 

_ P31  P33 _ 

a  b 

Pij  "  Pij' 

pij  =  f 

‘12 


for  its  first  2  columns 


for  its  third  column 


a  ,  b 

.  .  +  p .  .  ,  p 

13  *13 


m .  .  f  .  . 

.  .  =  -^1  =  41 
13  Bjj  rrr 


(13) 


(14a-e) 


Me  =  diag.  mjLi  ,  Fg  =  diag.  f£i 

Koto  that  P31e»  P32g  '  p13e  '  p23e  =  0  an  the  n°  case> 

The  results  in  Sections  2.2.1,  2.2.2  apply  as  again  (14c)  is 
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be  incorporated  in  the  design,  to  guarantee  the  desired 
performance  under  these  circumstances. 

2.2.5  Design  procedure 


The  plant  of  Figure  1  has  four  control  inputs  (the  four 
6s)  and  only  two  outputs,  so  there  exist  two  degrees  of  freedom 
(Pi2/t2i)  available  for  optimization.  It  was  shown  in  Section 
2.2.3  how  these  may  be  chosen  for  optimum  division  of  the  control 
effort.  However  other  optimal  criteria  may  be  perferable  (see 
Section  5).  In  any  case,  it  is  assumed  that  U12,li21  have  been 
chosen,  which  determines  Pg.  The  various  modes  and  conditions 
for  which  the  system  is  to  operate  are  listed  by  the  designer, 
the  flight  conditions,  the  failure  cases,  the  disturbances 
(gusts  etc)  etc.  These  give  the  sets  PQ  ,  V  of  Sections  1.1, 

2.2.1  and  the  performance  tolerances  T,  T^of  Section  1.1.  The 
design  technique  of  (Horowitz  1979,  1982b)  is  then  applied, 
resulting  in  required  m^(s),  fj^(s),  i=l/2.  The  values  of  the 
mij'^ij  ^or  are  avaiiabie  from  (5a, b)  and  the  design  is 

complete. 

*  * 

3.  A  FIVE  INPUT  -  THREE  OUTPUT  PLANT 

The  slight  increase  in  plant  complexity  in  Figure  2  results 

in  a  more  difficult  problem  of  selection,  if  control  effort 

division  is  the  optimality  criterion.  The  added  single  control 

input  and  the  Pair  (5®, 5®)  are  assumed  all  of  one  type  say 

lateral  e.g.  $2  is  a  rudder,  (6^,6^)  is  an  aileron  pair  while 
Si  b 

( 5  ^ ,  <5  ^ )  is  a  longitudinal  type,  e.g.  elevator  pair.  If  so  the 
P,M,F  matrices  are 
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«-•  *-• 


.  -  .  -v  -  .  ... . 


.  .  V.V, 


~  P22/M21  (P22+IJ12^21^  t‘lie  contro*  frequency  range  v/here 
1 12  ( juj)  [  >>1~y .  Equations  (10,11)  can  be  used  to  obtain  the 
control  effort  ratios  for  other  failure  cases. 


Thus,  if  1^2' ^21  provide  the  optimum  division  of  control 
efforts  between  6^  and  <S2  i-n  the  no  fail  case,  then  in  view  of 
the 'above,  the  division  under  failure  may  no  longer  be  optimum. 
But  in  failure  cases,  the  dominating  concern  is  acceptable 
performance,  giving  more  time  for  reliable  identification  and 
switching  in  of  compensations  optimum  for  the  failure  situation, 
should  this  be  considered  justifiable. 

2.2.4  Performance  under  failures 

The  quantitative  MIMO  synthesis  techniques  (Horowitz  1979, 
1932b)  guarantee  the  achievement,  theoretically,  of  arbitrarily 
narrow  performance  tolerances  for  a  large  class  of  uncertain 
MIMO  plants.  The  design  may  not  be  practical  if  the  sensors  are 
not  sufficiently  accurate  or  have  unduly  large  noise  levels. 

The  principal  condition  for  this  theoretical  ability  is  that 
the  MIMO  plant  is  minimqm-phase  over  its  range  of  uncertainty 
i.e.  det.  P£  has  no  right  half-plane  zeros  and  that  Pe  is 
output  controllable  (Horowitz  and  Shaked  1975,  Horowitz  1982b) 
meaning  that  despite  the  failures,  the  remaining  control  sur¬ 
faces  still  have  the  physical  ability  to  achieve  two  indepen¬ 
dent  outputs  (of  course  for  a  smaller  output  range) .  This  is 
so  in  the  present  problem  for  any  single,  or  any  two  simultan¬ 
eous  control  surface  failures.  If  there  are  three  simultaneous 
failures,  then  only  one  output  is  achievable  and  the  second 
output  is  not  independently  achievable.  Such  cases  too  can 


Appendix  C:  Equivalent  Plant  Transfer  Function  Matrices 


The  equivalent  plant  transfer  function  matrices 
are  listed  in  this  appendix  for  each  flight  condition  and 
control  surface  configuration  mode.  These  matrices  are 
derived  from  the  Pe  matrix  below  using  the  basic  plant 
transfer  function  equations  developed  in  Chapter  II. 

P  = 

— e 


P11  +  yP12 


Api2  + 

AP21  +  yAP22  P22  +  yP21 


where 


P 


ij 


P.  .a  +  P.  .b  and 
il  il 


The  reciprocal  plant  expressions  (1/Q^j)  for  the  SISO 

systems  of  Figure  4-5  are  obtained  by  inverting  the 

matrices.  These  transfer  functions  are  expressed  in  terms 

of  the  Q .  .  equations  in  this  appendix.  P  1  =  (1/Q-  •). 
i  j  e  1 3 
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TABLE  C-l 

EQUIVALENT  PLANT  MATRICES  FOR  FLIGHT  CONDITION  ONE 


P  * 
lie 


P  * 
F12e 


-2.176  (0)  (-.0175) (-.4617) 
(.363)  (-1.3)  (-. 0768+ j. 207) 


P  * 
21e 


P  * 
*22e 


-5.549(0)  (-.2237±il.029) 
(-.104)  (-.684) (-.274±jl. 91) 


Mode 

2: 


-1.058(0)  (-.01677)  (-.4669) 
(.363) (-1.3) (-.0768±j.207) 


.2795(0)  (-.01822)  (-.4568) 
(.363)  (-1.3)  (-.0768±j.207) 


2.142(0) (-. 3017+ jl. 562) 


(-.104) (-.684) (-.274±jl.91) 


-5.041(0)  (-.2153+1.9547) 


(-.104)  (-.684) (-.274+jl.91) 


2.206(0) (-.01787K-.4592) 


(.363) (-1.3) (-.0768+ j. 207) 


-.121(0)  (-.0654)  (-.259) 


(.363)  (-1.3)  (-.0768±j.207) 


,56(0)  (-.21±j.853) 


-3.3(0)  (-.2363111.13) 


TABLE  C-l — Continued 


Plle 

P12e 

\  -  * 

m 

-1.09(0)  (-.0175)  (-.462) 

.159(0)  (-.003)  (-.587) 

(.363)  (-1.3)  (-. 0768+ j. 207) 

(.363)  (-1.3)  (-.0768± j.207) 

Mode 

£: 

P21e 

P 

22e 

”  ' 

2.7(0) (-.282±jl.45) 

-2.78(0)  (-.224+jl.03) 

(-.104)  (-.684) (-.274+jl. 91) 

(-.104) (-.684) (-.274± jl. 91) 

Plle 

P12e 

-1.09(0)  (-.0175)  (-.462) 

.4(0)  (-.027)  (-.402) 

(.363)  (-1.3)  (-.07681  j. 207) 

1.363)  (-1.3)  (-. 0768+ j.  207) 

Mode 

5: 

t 

P21e 

P22e 

1.582(0)  (-.336+11.74) 

-2.78(0) (-.224+jl.03 

(-.104) (-.684)  (-.274+jl. 91) 

1-.104) (-.684) (-.274+jl. 91) 

• 

Plle 

P12e 

h  * 

-2.236(0)  (-.0182)  (-.457) 

(.363)  (-1.3)  (-.0768+ j. 207) 

Mode 

6: 

P21e 

P22e 

-1.07(0) (-.302±il.56) 

(-.104) (-.684) (-.274± jl.  91) 

**  7*,' 

TABLE  C-2 


EQUIVALENT  PLANT  MATRICES  FOR  FLIGHT  CONDITION  TWO 


P  * 

^lle 

P  * 

x2e 

-5.92(0)  (-.0099)  (-.562) 

• 

(1.167)  (-2.02)  (-. 0065+ j. 078) 

Mode 

1: 

P  * 

*21e 

P  * 

*22e 

-20.8(0)  (-.223+H.697) 

(.078)  (-.827)  (-.211±jl. 95) 


Plle 

-2.98(0)  (-.0098)  (-.574) 

P12e 

.73(0)  (-.01)  (-.55) 

Mode 

2: 

(1.167) (-2.02) (-.0065±j.078) 

(1.167) (-2.02) (-.0065±j.078) 

P21e 

P22e 

6.79(0) (-.244±i2.101) 

-19.14 (0) (-.221+11.66) 

(-.078) (-.827) (-.211±jl.95) 

(.078)  (-.827)  (-.211± jl. 95) 

“lie 


-5.89(0)  (-.01)  (-.56)  .106(0)  (-.0067)  (-1.86) 

(1.167) (-2.02) (-.0065± j .078)  (1.167)  (-2.02)  (-. 0065+ j. 078) 

Mode 

3: 

P21e  P22e 


2.18(0) (-.219ijl.61) 


-12.12(0) (-.22±11.7) 


TABLE  C-2 — Continued 


Plle 

-2.96(0)  (-.01)  (-.56) 

P12e 

.84(0)  (-.01)  (-.72) 

Mode 

4s 

(1.167)  (-2.02) (-.0065+ j. 078) 

P21e 

8.97(0) (-.24±il.99) 

(1.167)  (-2.02)  (- . 00651 j. 078) 

P 

*22e 

-10.4(0)  (-.22±il.7) 

(-.078)  (-.827) (-.211±jl.95) 

(.078) (-.827) (-.2111 jl. 95) 

Plle 

P12e 

-2.96(0)  (-.01)  (-.56) 

.63(0)  (-.013)  (-.326) 

Mode 

5: 

(1.167)  (-2.02) (-.0065+ j. 078) 

(1.167)  (-2.02)  (-.00651;). 078) 

P21e 

P22e 

4.6(0)  (-.256±i2.3) 

-10.4(0)  (-.22H1.7) 

(-.078) (-.827) (-.211±jl.95) 

(.078) (-.827) (-.211ijl.95) 

-2.93(0)  (-.01)  (-.55) 

(1.167)  (-2.02) (-.0065±j.078) 

Mode 

6: 

P21e  P22e 


-1.7(0)  (-.244±j2.1) 


(.078)  (-.827)  (-.211±jl. 95) 


TABLE  C-3 


EQUIVALENT  PLANT  MATRICES  FOR  FLIGHT  CONDITION  THREE 


P  * 
lie 


Mode 

1: 


-25.7(0)  (-.013)  (-1.52) 
(.965)  (-3.22) (-.0076±j.538) 


-63.7(0)  (-.358±i3.07) 
(-.0272  (-2.7)  (-.39±j2.96) 


Mode 

2: 


lie 

-13.7(0)  (-.0126)  (-1.526) 
(.965) (-3.22) (-.0076±j.538) 


21e 


25.4(0)  (-.375±i3.58) 
(-.0272)  (-2.7) (-.39±j2.96) 


12e 

3.75(0)  (-.013)  (-1.52) 

(.965)  (-3.22) (-.0076+ j. 538) 


22e 


-57.4(0)  (-.354±i3.01 
(-.0272  (-2.7) (-.39±j2.96) 


P  P 

lie  12e 


-24.9(0)  (-.0126M-1.514)  3.24 (0) (-.013) (-1.65) 

(.965) (-3.22) (-.0076+ j. 538)  (.965)  (-3.22)  (-.0076±j.538) 

Mode 

3: 

P  p 

21e  22e 


6.38(0)  (-.354±j2.93) 


-22.9(0)  (-.352+j2.84) 


TABLE  0  3 — Continued 


Plle 

-12.8(0)  (-.0126)  (-1.52) 

P12e 

6.24(0)  (-.0125)  (-1.58) 

Mode 

4s 

(.965) (-3.22) (-.0076±j.538) 

P21e 

31.7(0)  (-.37+13.46) 

(.965)  (-3.22) (-.0076+ j. 538) 

P 

22e 

-31.8(0)  (-.358±j3.07) 

(-.0272)  (-2.7) (-.39±j2.96) 

(-.0272) (-2.7) (-.39±j2.96) 

Plle 

P12e 

-12.8(0)  (-.0126)  (-1.52) 

-.229(0)  (-.012)  (-3.44) 

Mode 

5: 

(.965) (-3.22) (-. 0076+ j. 538) 

(.965)  (-3.22)  (-. 0076+ j. 538) 

P21e 

P22e 

*  4 


19(0)  (-.382+j3.77) 


31.8(0)  (-.358±j3.07) 


TABLE  C-4 


EQUIVALENT  PLANT  MATRICES  FOR  FLIGHT  CONDITION  FOUR 


.035)  (-2.17)  <-5±j3.13) 


♦Equivalent  Plant 


TABLE  C-4 — Continued 


P 


lie 


P 


12e 


Mode 

4: 


-17.2(0)  (-.03)  (-1.08)  7.04(0)  (-.03)  (-.927) _ 

(-.0152+ j. 0234)  (-8.01±j6.6)  (-.0152± j.0234)  (-.801±j6.6) 


P 


21e 


P 


22e 


25.1(0)  (-.38±j3.93) 
(-.035)  (-2.17)  (— 5± j 3 . 13) 


-12.9(0)  (-.394±j4.47) 
(-.035)  (-2.17)  (-.5±j3.13) 


12e 


Mode 

5: 


-17.2(0)  (-.03)  (-1.08)  1.19(0)  (-.026)  (-2.1) 

(-.0152± j.0234)  (-8.01±j6.6)  (-.0152± j.0234)  (-.801±j6.6) 


21.5(0) (-.37±i3.75)  -12.9(0) (-.39±j4.45) 

(-.035)  (-2.17) (-5±j3.13)  (-.035)  (-2.17) (-.5±j3.13) 


Plle  P12e 


-32.9(0)  (-.03)  (-1.1)  __ 

(-.0152± j.0234)  (-8.01+ j6. 6) 

Mode 

6; 

P21e  P22e 


-11.7(0)  (-,38±j3.85) 
(-.035)  (-2.17)  (— . 5± j3. 13) 
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TABLE  C-5 — Continued 


Lower  Bound — a 


22 

£  _  _ 6.937 _ 

R  (s  +  0.75)  (s  +  2.5)  (s  +  3.7) 

Figures  of  Merit: 

Peak  Time,  T  =  Very  Large 

ir 

Settling  Time,  Tg  =  5.99  sec 
Peak  Value,  =  50.0 

Loop  interaction  Bound — b^2 

3  _  0.04s  (s  +  1)  (s  +  10) _ 

R  (s  +  0.75) (s  +  5) (s2  +  2s  +  2) 

Figures  of  Merit: 

Peak  Time,  Mp  =  0.992  sec 

Settling  Time,  T  =4.30  sec 

s 

Peak  Value,  M^  =  1.4  deg /sec 


Loop  Interaction  Bound — b2- 

£  _  0.335s  (s  +  1) (s  +  : 

R  (s  +  0.75) (s  +  5) {s2  +  2s  +  2) 

Figures  of  Merit: 

Peak  Time,  =  0.992  sec 

Settling  Time,  T  =4.30  sec 

s 

Peak  Value,  M^  =  2.344  deg/sec 

Roll  Rate  Models 

Upper  Bound — b22 

p  _  2.136  (s+1 . 3) _ 

R  (s2  +  2.667s  +  2.776) 

Figures  of  Merit: 

Peak  Time,  T^  =  1.538  sec 

Settling  Time,  T  =3.01  sec 

s 

Peak  Value,  =  55.0  deg /sec 
P 
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Pitch  Rate  Models 


Upper  Bound — 

5  _  4.167  (s  +  2.667) 

R  (s2  +  5.334s  +  11.113) 

Figures  of  Merit: 


Peak  Time,  T  =0.785  sec 
P 

Settling  Time,  T  =  1.5  sec 

s 

Peak  Value,  M^  =  11.0  deg /sec 
Lower  Bound — a^ 


S  = 
R 


63.75 

(s  +  1.5)  (s  +  4.25)  (s  +  10) 


Figures  of  Merit: 


Peak  Time,  Tp  =  Very  Large 

Settling  Time,  T  =  3.01  sec 

s 

Peak  Value,  M  =10  deg /sec 
P 


Appendix  D:  Pitch  Rate  and  Roll  Rate 


Response  Models 


This  appendix  contains  the  transfer  function 
models  for  pitch  rate  and  roll  rate.  These  models  are 
synthesized  with  desired  specifications  in  mind,  and  the 
models  are  used  to  develop  bounds  on  the  loop  transmission 
Also  listed  in  this  appendix  are  the  figures  of  merit  for 
the  response  models.  See  Chapter  IV  for  the  time  and 
frequency  domain  plots  of  these  models. 
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TABLE  C-8 — Continued 


TABLE  C-8 — Continued 


TABLE  C-7 — Continued 


TABLE  C-7 — Continued 


EQUIVALENT  PLANT  EQUATIONS  FOR  THE  SINGLE  INPUT-SINGLE  OUTPUT  SYSTEMS 
FOR  FLIGHT  CONDITION  THREE  (P”1  =  [11Q. ■]) 


♦Equivalent  Plant. 


TABLE  C-5 — Continued 


*■.  n 


^7 


O 

o 

o 


Ctl 


'  •360.  CO 


! 


!  F  -  8  .  G  1 


( FAILURES  ) 

1  -  SCNC 

2  -  D£LM  !fl 

3  -  OEL'n  38 

1  -  Btl'H  11.28 
S  -  OtLTI  lfi.28 
8  -  GEt’1  21*28 


-3CC.CC  -2*0. GG  -180. CO  -120. CO 

DEGREES 


-SO. CO 


Fig.  E-3  (a)  . 


Template  of  at  to 


8  Rad/Sec 


O 

O 


O 

O 


OO 
3-  cs» 


o 

o 


G ZZ  IF  =  8.0) 


°  I 

7-1 - t - 

-360. CO  -300. GO 


TTioTcO  -180.  CO 

DEGREES 


(FqiLURES) 

t  -  *iopie 

2  -  0EL »ft  »ft 

3  -  DLL  *ft  2ft 

4  •  OUffl  lft»2ft 

5  -  Otirfl 

6  -  OCLTft  2ft»28 

-120.00  -60.00 


Fig.  E-3(b).  Template  of  at  to  =  8  Rad/Sec 


Appendix  F:  Simulation  Programs 

The  computer  programs  listed  in  this  appendix  are 
the  programs  used  to  conduct  the  simulation  of  the  system 
represented  by  Figure  2-1.  The  simulation  is  structured 
around  a  main  program  and  three  subroutines.  The  main 
program  is  called  Zsimu;  Zsimu  calls  subroutines  Derv, 
Zdata,  and  Degol.  Poles  and  zeros  of  the  basic  plant  equa¬ 
tions  are  listed  in  Zdata. 

Data  for  the  four  basic  plant  equations  (p^j)  at 
each  flight  condition  are  listed  in  Zdata.  These  data  are 
passed  as  poles  and  zeros  to  Zsimu.  Derivatives  of  the 
output  variables  are  passed  from  Derv  to  Zsimu. 

Subroutine  Derv  uses  the  state-variable  format  to 
write,  basic  plant,  prefilter,  and  compensator  equations. 
Desired  outputs  are  derived  from  the  state  equations  based 
upon  the  structure  of  the  system  being  simulated.  For  this 
simulation,  a  maximum  of  30  state-variables  can  be  defined 
in  Derv.  This  limitation  is  due  to  the  subroutine  Degol. 
Position  saturation  for  the  elevators  and  flaperons  is 
also  included  in  this  subroutine.  Derivatives  of  the  out¬ 
put  variables  are  passed  to  the  main  program  for  integra- 


The  main  program  Zsimu  takes  the  outputs  of  Derv 


and  integrates  these  using  the  subroutine  Degol.  Degol  is 
an  integration  routine  for  linear,  ordinary  differential 
equations  similar  to  the  program  ODE.  The  "case"  block  of 
Zsimu  allows  inputs  (IR) ,  flight  conditions  (IC) ,  or  CSC 
modes  (IM)  to  be  selected  as  desired.  The  "integration" 
block  allows  the  integration  time  and  the  number  of  points 
in  the  int  gration  to  be  specified.  For  this  simulation  the 
outputs  of  Zsimu  are  pitch  rate,  roll  rate,  and  the  four 
surface  deflections.  These  outputs  are  identified  by  the 
variables  YYl  to  YY6,  respectively.  Integration  time  is 
denoted  by  the  variable  XX.  If  other  outputs  are  desired, 
they  are  defined  in  Derv.  A  step  or  ramp  command  input  may 
be  specified  in  the  "Degol"  block. 


L6 


ooooo  o  n  o  o  o  ooo  nnooooooo 


Zsimu 


PROGRAM  Z  S IMU ( I NP  UT , OUTP  UT , DAT ALP , D AT AOUT , PLOT , TAPE  9 , 
+TAPE6  OUTPUT ,  T  AP  E  7  =D  AT  ALP  ,  TAPE8=DATA0UT) 


UPDATE %  AUG.  15,  1984 

ORIGINAL%  AUG. 12,  1984 

***************  DIMENSION  *************** 


—  DEGOL  - 

REAL  Y (30) 

EXTERNAL  DERV 

DIMENSION  XX (500) ,YY(500)  ,YYl (500) ,YY2  (500) ,YY3  (500) , 
DIMENSION  YY4  (500) ,YY5  (500)  ,YY6  (500) 

—  PLANT  - 

DIMENSION  VP11K  (4) ,VPllZl (4)  ,VP11Z2  (4) ,VP11Z3  (4) 

DIMENSION  VP12K (4)  ,VP12Z1  (4) ,VP12Z2(4)  .VP12Z3 (4) 

COMPLEX  VP11P1 (4)  ,VP11P2 (4) ,VP11P3  (4) ,VP11P4  (4) 

DIMENSION  VP21K  (4) ,VP21Z1 (4) 

COMPLEX  VP21Z2 (4) ,VP21Z3 (4) 

*  DIMENSION  VP22K  (4 )  ,VP21() 

COMPLEX  VP22Z2 (4) ,VP22Z3 (4) 

DIMENSION  VP22P1  (4) ,VP22P2 (4) 

COMPLEX  VP27P3 (4) ,VP22P4 (4) 


***************  COMMON  *************** 


DERV  - 

COMMON /P11Z/  P11K,P11Z1 ,P11Z2 , P11Z 3 
COMMON/P11P/  P11P1 ,P11P2 ,P11P3 ,P11P4 
C 

COMMON /P 1 2 Z /  P12K,P12Z1,P12Z2,P12Z3 
C 

COMMON /P21Z/  P21K,P21Z1,P21Z2,P21Z3 


on  o o n o o  o  non  non  no  on  non  on 


COMMON  /P22Z/  P22K,P22Z1 ,P22Z2 ,P22Z3 
COMMON  /P22P /  P22P1 ,P22P2 ,P22P3 ,P22P4 


COMMON  /F/  FlK,FlP,F2K,F2P 

COMMON  /G1Z /GlK ,G1Z1 , G1Z2 , G1Z  3 ,G1Z4 
COMMON  /GlP/GlPl ,G1P2 ,G1P3 ,GlP4 

*  COMMON  /G2Z /G2K ,G2Z1 ,G2Z2,G2Z3,Z2Z4 
COMMON  /G2P/G2P1 ,22G2P3 ,G1P4 

COMMON  /R/  R1,R2 

COMMON  /U/  U12 , U21 

COMMON  /Cl  01,02 

COMMON  /D/  DlA,DlB,D2A,D2B 

COMMON  /IS/  IS1A, ISlB, IS2A, IS2B 

-  MONITOR  - 

COMMON  /M/  F1,F2,G1,G2,E1,E2,D11,D12,D21,D22 
COMMON  /M/  Cll ,C12 ,C21 , C22 
-  ZDATA  - 

COMMON  /VP11Z/VP11K,VP11Z1,VP11Z2 ,VP11Z3 
COMMON  /VP11P/VP11P1 ,VP11P2 ,VPllP3 ,VPllP4 

COMMON  /VP12Z /VP12K,VP12Z1 ,VP12Z2 ,VP12Z3 

COMMON  /VP21Z/VP21K,VP21Z1,VP21Z2,VP21Z3 

COMMON  /VP22Z /VP22K,VP22Z1 ,VP22Z2 , VP22Z 3 
COMMON  /BP22P/VP22P1 , VP22P2 ,VP22P3 ,VP22P4 

-  CASE  - 

DIMENSION  IRV (30) ,ICV (30) ,IMV (30)  • 

DATA  NC/1/ 

DATA  IRV/1, 1,1, 1,1, 1,2, 23*0/ 

DATA  ICV/2 ,2,4,2,4,2,2,23*0/ 

DATA  IMV/4 ,2,3,4,5,6,7,23*0/ 

***************  MAIN  PROGRAM  *************** 

-  P  (S)  - 

CALL  ZDATA 

-  F1(S)  - 
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F1K  =  2. 

FlP  =  2. 

c 

c— 

-  F2  (S)  - 

.  -  . 

F2K  =  2. 

c 

F2P  =  2. 

c 

■ 

c— 

-  G1(S)  - 

„  -  -  » 

G1K  =  -660000. 

c 

* 

G1Z1  =  4. 

G1Z2  =  100. 

G1Z  3  =  .013 

G1Z4  =  .006 

“* 

c 

G1Z3  =  .013  +  .006 

c 

G1Z4  =  .013*. 006 

c 

GlPl  =  .0223 

G1P2  =  25. 

c 

G1P3  =  630. 

m  i 

c 

G1P4  =  630. 

G1P3  =60  +  630. 

G1P4  =  630. *630. 

'  • 

c 

c— 

-  G2  (S)  - 

G2K  =  -36800. 

- 

G2Z1  =  3. 

**  "■  m  ’ 

G2Z2  =  900. 

G2Z3  =  2. 

G2Z4  =  4. 

c 

G2Z 3  ■  1.  +  10. 

c 

G2Z4  =  10. 

"  “"1 

c 

• 

G2P1  =  .25 

G2P2  =  50. 

J 

c 

G2P3  =  158. 

c 

G2P4  =  3750. 

G2P3  =  158.  +  3750. 

G2P4  =  158. *3750. 

—  " 

C 

C 

c -  INPUT  R  (I ) 


R1  -  10. 

R2  =  50. 

C 

C -  SCALE  U(I , J) 

U12  =  .25 
U21  =  .25 
C 
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*************** 


0* **************  DO  CASES 


C 


DO  10  III  =  1,NC 


C 


C 

C 

C 

C 


C 

c 


c 

c 


IR  =  IRV  (III) 

IC  =  ICV(III) 

IM  =  IMV(III) 

R1  =  0. 

R2  =  0. 

FIR. EQ.  1)  Rl  =  10. 

IF  (IR  .EQ.  2)  R2  *  50. 


IS1A  =  1 

IS1B  =  1 

IS2A  =  1 

IS2B  =  1 

IF ( IM  .EQ. 

2) 

IS1A 

IF (IM  .EQ. 

3) 

IS2A 

IF (IM  .EQ. 

4) 

IS1A 

IF (IM  .EQ. 

4) 

IS2A 

IF (IM  .EQ. 

5) 

IS1A 

IF  (IM  .EQ. 

5) 

IS2B 

IF  (IM  .EQ. 

6) 

IS2A 

IF  (IM  .EQ. 

6) 

IS2B 

IF  (IM  .EQ. 

7) 

IS1A 

IF (IM  .EQ. 

7) 

IS1B 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


Pll(S)  - 

PllK  =  VP11K (IC) *. 5 
P11Z1  =-VPllZl(IC) 

P11Z2  =  VP11Z2 (IC) 

P11Z3  =  VP11Z3  (IC) 

P11Z2  =  -VP11Z2 (IC)  -  VP11Z 3 (IC) 
P11Z3  =  VP11Z2 (IC)  1Z3  (IC) 


C 

C  P11P1  =  -2 .  *REAL  (VPUPl  (IC)  ) 

C  P11P2  =  CABS  (VPUPl  (IC)  )  *2 

P11P1  =  -REAL  (VPUPl  (IC) +VP11P2  (IC)  ) 
P11P2  =  REAL(VP11P1  (IC)  *VP11P2  (IC)  ) 
P11P3  =  -2.*REAL  (VP11P3  (IC) ) 

P11P4  =  CABS (VP11P3 (IC) ) *2 

C -  P12  (S)  - 

P12K  =  VP12K (IC) *. 5 
P12Z1  =  -VP12Z1 (IC) 

C  P12Z2  =  VP12Z2 (IC) 

C  P12Z3  =  VP12Z3 (IC) 

P12Z2  =  -VP12Z2  (IC)  -  VP12Z 3 (IC) 
P12Z3  =  VP12Z2 (IC) *VP12Z3 (IC) 


C 

0 -  P21(S)  - 

P21K  =VP21K(IC)  *.5 


P21Z1  =  -VP21Z1  (IC) 

P21Z2  =  -2 . *REAL (VP21Z2 (IC) ) 

P21Z3  =  CABS (VP21Z2 (IC) ) *2 
C 

c -  P22  (S)  - 

P22K  =VP22K (IC) *. 5 
P22Z1  =  -VP22Z1 (IC) 

P22Z2  =  -2. *REAL (VP22Z2 (IC) ) 

P22Z3  =  CABS (VP22Z2 (IC) ) *2 
C 

C  *  P22P1  =  VP22P1 (IC) 

C  P22P2  =  VP22P2 (IC) 

P22P1  =  -VP22P1  (IC)  -  VP22P2 (IC) 

P22P2  =  VP22P1 (IC) *VP22P2 (IC) 

P22P3  =  -2 . *REAL (VP22P3 (IC) ) 

P22P4  =  CABS (VP22P3 (IC) ) *2 
C 
C 

C* **************  INTEGRATION  *************** 
C 

NE  =  28 

AE  =  l.E-5 

RE  =  l.E-6 

MAX  =  1000 
H  =  l.E-4 
N  =  60000 
NP  =  300 
TO  =  0. 

Tl  =  H 
C 

C -  INITIAL  CONDITION  - 

Cl  =  0. 

C2  =  0. 

*D1A  =  0. 

DIB  =  0. 

D2A  =  0. 

2  =  0.  • 

C 

DO  30  I  =  1 ,NE 


30  Y  (I)  =  0. 

C 

C -  ITERATION 


II  =  0 
IIXY  =  1 
KKXY  =  N/NP 

IF (KKXY  .LT.  1)  KKXY  =  1 
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NL 

C -  YY (II)  - 

IF (I  .EQ.  IIXY)  THEN 
II  =11+1 

C  WRITE (8,122) II , Gl ,G2 ,Dll ,D12 ,D21 ,D22 

C  WRITE (8,122) II, Cl, C2,C11,C12,C21,C22 


PRINT  122, II, TO, Cl, C2 , DlA ,D1B ,D2A,D2B 
XX  (II)  =  TO 
YYl  (II)  =  Cl 
YY2 (II)  =  C2 
.  YY  3  ( 1 1 )  =  DlA 

YY4  (II)  =  DIB 
YY5  (II)  =  D2A 
YY6  (II)  =  D2B 
IIXY  =  IIXY  +  KKXY 
END  IF 
C 


C -  CALL  DEGOL  - 

C  PRINT  122, I, TO, Cl, C2 ,D1A,D1B,D2A,D2B 


IK  =  1 

CALL  DEGOL (DERV,NE , Y,T0 ,Tl ,RE, AE , IK, MAX) 
IF (IK  .NE.  2)  THEN 
PRINT  * , IK, TO ,T1 
STOP 
END  IF 
T1  =T1  +  H 
C 

C  IF  (T1  .GT.  5.)  R1  3  0. 

C  IF  (Tl  .GT.  5.)  R2  =  0. 

C  IF (IR  .EQ.  1  .AND.  Tl  .LE.  .5)  Rl  =  20.*T1 

IF (IR  .EQ.  1  .AND.  Tl  .GT.  .5)  Rl  =  10. 

IF (IR  .EQ.  1  .AND.  Tl  .GE.  5.)  Rl  =  0. 

C 

.  IF (IR  .EQ.  2  .AND.  Tl  .LE.  1.)  R2  =  50.*Tl 
IF (IR  .EQ.  2  .AND.  Tl  .GT.  1.)  R2  =  50. 

IF (IR  .EQ.  2  .AND.  Tl  .GE.  5.)  R2  =  0. 

C 

40  CONTINUE 
C 

C******** *******  PLOT  *************** 

NY  =  II 

PRINT  * ,NXY 
C 

WRITE (8, 121) NXY 
C  WRITE (7, 121) NXY 

DO  50  I  =  1 ,NXY 

WRITE (8,122) I ,XX (I ) ,YYl (I) ,YY2 (I)  ,YY3  (I) 
YY4 (I) , YY5 (I) , YY6 (I) 

WRITE  (7,122)  I, XX  (I)  ,YY1(I)  ,YY2  (I)  ,YY3(I) 
YY4 (I) ,YY5 (I ) ,YY6 (I) 

50  CONTINUE 
C 

CALL  LPL (NXY ,  XX,  YYl,IM) 


n  o 


c -  end  cases  loop  - 

10  CONTINUE 
C 

C* **************  print  out  *************** 

c 

CALL  LPRINT(l) 

CALL  LPLl (NXY,XX,YY2 ,1) 

CALL  LPRINT (2) 


c***************  format  *************** 

121  FORMAT (IX, I 5) 

122  FORMAT (1X,I5,7 (1X,E9. 3)) 

C 


RETURN 


n  n  n  n  n  o  o 


SUBROUTINE  DEGOL 


Degol 


SUBROUTINE  DEGOL (F , NEQN , Y , T , TOUT , RELERR, ABSERR , I FLAG , 
ITMAX) 

.  INTEGER  K 
REAL  HOLD 

LOGICAL  START, CRASH, STIFF 
DIMENSION  Y (NEQN) ,PSI (12) 

DIMENSION  YY (30) ,WT(30) ,PHI (30,16) ,P(30) ,YP (30) , 

YPOUT (30) 

EXTERNAL  F 

DATA  FOURU/ . 44E-15/ 

MAXN  UM  =M AX  0  (500, ITMAX) 

IF  (NEQN  .LT.  1  .OR.  NEQN  .GT.  30)  GO  TO  10 
IF (T  .EQ.  TOUT)  GO  TO  10 

IF (RELERR  .LT.  0.0  .OR.  ABSERR  .LT.  0.0)  GO  TO  10 
EPS  =  AMAXl (RELERR, ABSERR) 

IF (EPS  .LE.  0.0)  GO  TO  10 
IF (IFLAG  .EQ.  0)  GO  TO  10 
ISN  =  ISIGN(1, IFLAG) 

IFLAG  =  IABS (IFLAG) 

IF (IFLAG  .EQ.  1)  GO  TO  20 

IF (IFLAG  .GE.  2  .AND.  IFLAG  .LE.  5)  GO  TO  20 
10  IFLAG  =  6 
RETURN 

20  DEL  =  TOUT  -  T 

ABSDEL  =  ABS  (DEL) 

.TEND  =  T  +  10 . 0*DEL 
IF (ISN  .LT.  0)  TEND  =  TOUT 
NOSTEP  =  0 
KLE4  =  0 
STIFF  =  .FALSE. 

RELEPS  =  RELERR/EPS 

ABSEPS  =  ABSERR/EPS 

IF (IFLAG  .EQ.  1)  GO  TO  30 

IF(ISNOLD  .LT.  0)  GO  TO  30 

IF (DELSGN*DEL  .GT.  0.0)  GO  TO  50 

ON  START  AND  RESTART  ALSO  SET  WORK  VARIABLES  X  AND  YY (*) 
STORE  THE  DIRECTION  OF  INTEGRATION  AND  INITIALIZE  THE 
STEP  SIZE 

30  START  =  .TRUE. 

X  =  T 

DO  40  L  ■  1 ,NEQN 


DELSGN  =  SIGN (1 . 0 ,DEL) 

H  =  SIGN (AMAX1 (ABS (TOUT-X) ,FOURU*ABS (X) ) , TOUT-X) 

IF  ALREADY  PAST  OUTPUT  POINT,  INTERPOLATE  AND  RETURN 

50  IF (ABS  (X-T)  .LT.  ABSDEL)  GO  TO  60 

CALL  INTRP (X , YY , TOUT , Y , YPOUT , NEQN , KOLD , PHI , PS I ) 

I FLAG  =  2 
T  =  TOUT 
*  TOLD  =  T 
ISNOLD  =  ISN 
RETURN 

IF  CANNOT  GO  PAST  OUTPUT  POINT  AND  SUFFICIENTLY  CLOSE 

EXTRAPOLATE  AND  RETURN 

60  IF (ISN  .GT.  0  .OR.  ABS (TOUT-X)  .GE.  FOURU*ABS (X) ) 

GO  TO  80 
H  =  TOUT  -  X 
CALL  F  (X , YY , YP) 

DO  70  L  =  1 , NEQN 

70  Y (L)  =  YY (L)  +  H*YP (L) 

I FLAG  =  2 
T  =  TOUT 
TOLD  =  T 
ISNOLD  =  ISN 
RETURN 

TEST  FOR  TOO  MUCH  WORK 

80  IF (NOSTEP  .LT.  MAXNUM)  GO  TO  100 
I FLAG  =  ISN*4 
‘IF (STIFF)  I FLAG  -  ISN*5 
DO  90  L  =  1 ,NEQN 

90  Y  (L)  =  YY  (L) 

T  =  X 
TOLD  =  T 
ISNOLD  =  1 
RETURN 

LIMIT  STEP  SIZE,  SETWEIGHT  VECTOR  AND  A  STEP 

100  H  =  SIGN (AMIN1 (ABS (H) , ABS (TEND-X) ) ,H) 

DO  110  L  =  1 ,NEQN 

110  WT (L)  =  RELEPS  *ABS (YY (L) )  +  ABSEPS 
CALL  STEP (X,YY,F, NEQN, H, EPS, WT, START, 

+  HOLD, K, KOLD, CRASH, PHI, P,YP, PSI) 

TEST  FOR  TOLERANCE  S  TOO  SMALL 
IF (  .NOT.  CRASH)  GO  TO  130 
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I FLAG  =  ISN*3 
RELERR  =  EPS*RELEPS 
ABSERR  =  EPS*ABSEPS 
DO  120  L  =  1 /NEQN 
120  Y (L)  =  YY  (L) 

T  =  X 
TOLD  =  T 
ISNOLD  =  1 
RETURN 

AUGMENT  COUNTER  ON  WORK  AND  TEST  FOR  STIFFNESS 

130  NOSTEP  =  NOSTEP  +  1 
KLE4  =  KLE4  +  1 
IF (KOLD  .GT.  4)  KLE4  =  0 
IF (KLE4 .  GE.  50)  STIFF  =  .TRUE. 

GO  TO  50 
END 

*********************************************** 

SUBROUTINE  STEP 

*********************************************** 

SUBROUTINE  STEP  (X ,  Y ,  F , NEQN , H , EPS , WT , START , 

+  HOLD, K, KOLD, CRASH, PHI ,  P ,  YP , PS I) 

LOGICAL  START , CRASH , PHASEl , NORND 

DIMENSION  Y (NEQN) ,WT (NEQN) ,PHI (NEQN, 16) ,P  (NEQN) , 

+  YP (NEQN) ,PSI (12) 

DIMENSION  ALPHA  (12),BETA  (12)  ,SIG(13)  ,W(12)  ,V(12)  ,G(13) 
+  GSTR(13)  , TWO  (13) 

EXTERNAL  F 

DATA  TWOU/.22E-15/ 

DATA  FOURU/.44E-15/ 

DATA  TWO/2 . 0,4. 0,8. 0,16. 0,32. 0,64. 0,128. 0,256. 0,512.0 
+  1024.0,2048.0,4096.0,8192.0/ 

DATA  GSTR/0. 500, 0.0833, 0.0417, 0.0264, 0.0188, 0.0143, 

+  0.0114,0.00936,0.00789,0.00679,0.00592,0.00524, 

+  0.00468/ 

DATA  G (1) ,G(2) /l. 0 ,0 . 5/ ,SIG (1) /1.0/ 

CRASH  =  .TRUE. 

IF (ABS (H)  .GE.  FOURU*ABS (X) )  GO  TO  5 
H  =  SIGN (FOURU* ABS (X) ,H) 

RETURN 

5  P5EPS  =  0. 5*EPS 

IF  ERROR  TOLERANCE  IS  TOO  SMALL,  INCREASE  IT  TO  AN 
ACCEPTABLE  VALUE 

ROUND  *  0.0 
DO  10  L  =  1 , NEQN 

10  ROUND  =  ROUND  +  (Y  (L) /WT  (L)  )  **2 

ROUND  =  TWOU*SQRT (ROUND) 
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IF (P5EPS  .GE.  ROUND)  GO  TO  15 
EPS  =  2.0*ROUND* (1.0  +  FOURU) 

RETURN 

15  CRASH  =  .FALSE. 

IF (.NOT. START)  GO  TO  99 

INITIALIZE.  COMPUTE  APPROPRIATE  STEP  SIZE  FOR  FIRST 
C  STEP 

CALL  F (X, Y ,YP) 

SUM  =0.0 
*  DO  20  L  =  1 ,NEQN 
PHI  (L,  1)  =  YP(L) 

PHI (L,2)  =0.0 

20  SUM  =  SUM  +  (YP (L) /WT (L) ) **2 
SUM  =  SQRT (SUM) 

ABSH  =  ABS (H) 

IF (EPS. LT.  16. 0*SUM*H*H)  ABSH  =  0 . 25*SQRT (EPS/SUM) 

H  =  SIGN (AMAX1 (ABSH /FOURU*ABS (X) ) ,H) 

HOLD  =0.0 
K  =  1 
KOLD  =  0 
START  =  .FALSE. 

PHASE1  =  .TRUE. 

NORND  =  .TRUE. 

IF (P5SEPS 

.GT.  100. 0*ROUND)  GO  TO  99 
NORND  =  .FALSE. 

DO  25L  =  1 ,NEQN 
25  PHI (L , 15)  =  0.0 

99  IFAIL  =  0 

***  END  BLOCK  0  *** 

***  BEGIN  BLOCK  1  *** 

COMPUTE  COEFFICIENTS  OF  FORMULAS  FOR  THIS  STEP.  AVOID 
COMPUTING  THOSE  QUANTITIES  NOT  CHANGED  WHEN  STEP  SIZE  IS 

NOT  CHANGED. 

** * 


100  KP1  =  K+l 
KP2  =  K+2 
KM1  =  K-l 
KM2  =  K-2 

NS  IS  THE  NUMBER  OF  STEPS  TAKEN  WITH  SIZE  H,  INCLUDING 
THE  CURRENT  ONE.  WHEN  K.LT.NS,  NO  COEFFICIENTS  CHANGE 

IF (N  .NE.  HOLD)  NS  =  0 
NS  =  MIN0 (NS+l,KOLD+l) 

NSP1  =  NS+1 

IF (K  .LT.  NS)  GO  TO  199 


C  COMPUTE  THOSE  COMPONENTS  OF  ALPHA (*) , BETA (*) PSI (*) , 

C  SIG (*)  WHICH  ARE  CHANGED 
BETA (NS)  =1.0 
REALNS  =  NS 

ALPHA (NS)  =  1. 0 /REALNS 
TEMPI  *  H*REALNS 
SIG(NSPl)  =1.0 
IF (K  .LT.  NSP1)  GO  TO  110 
DO  105  I  =  NSP1 ,K 
IM1  =  1-1 

TEMP2  =  PSI(IMl) 

PSI(IMl)  =  TEMPI 

BETA (I)  =  BETA(IMl) *PSI (IMl) /TEMP2 
TEMPI  =  TEMP 2  +  H 
ALPHA  (I)  =  H /TEMPI 
REALI  =  I 

105  SIG (1+1)  =  REALI * ALP HA (I ) *SIG ( I ) 

110  PSI (K)  =  TEMPI 
C 

C  COMPUTE  COEFFICIENTS  G(*) 

C 

C  INITIALIZE  V(*)  AND  W(*),  G(2)  IS  SET  DATA  STATEMENT 

C  IF (NS  .GT.  1)  GO  TO  120 

DO  115  IQ  =  1,K 

TEMP 3  =  IQ* (IQ+1) 

V ( IQ)  =  1 . 0 /TEMP 3 
115  W (IQ)  =  V (IQ) 

GO  TO  140 
C 

C  IF  ORDER  WAS  RAISED,  UPDATE  DIAGONAL  PART  OF  V  ( *) 

C 

120  IF (K  .LE.  KOLD)  GO  TO  130 
. TEMP4  =  K*KPl 
V (K)  =  1. 0/TEMP4 
NSM2  =  NS- 2 

IF (NSM2  .LT.  1)  GO  TO  130 
DO  125  J  =  1 ,NSM2 
I  =  K-J 

125  V(I)  =  V(I)  -  ALPHA  (J+1)*V(I+1) 

C 

C 

C  UPDATE  V(*)  AND  SET  W(*) 

C 

130  LIMITl  =  KPl  -  NS 
TEMPS  =  ALPHA (NS) 

DO  135  IQ  =  1, LIMITl 

V (IQ)  =  V (IQ)  -  TEMP5*V (IQ+1) 

135  W (IQ)  =  V (IQ) 

G(NSP1)  =  W(l) 


non  nooooonoooo  ,  no 


COMPUTE  THE  G(*)  IN  THE  WORK  VECTOR  W(*) 


145 

150 

199' 


NSP2  =  NS  +2 

IF (KPl  .LT.  NSP2)  GO  TO  199 
DO  150  I  =  NSP2 ,KP1 
LIMIT2  =  KP2  -I 
TEMP  6  =  ALPHA (1-1) 

DO  145  IQ  =  1,LIMIT2 

W  (IQ)  =  W  (IQ)  -  TEMP6*W  (IQ+1) 
G  (I )  =  W(l) 

CONTINUE 

***  END  BLOCK  1  *** 


***  BEGIN  BLOCK  2  *** 

PREDICT  A  SOLUTION  P(*),  EVALUATE  DERIVATIVES  USING 
PREDICTED  SOLUTION,  ESTIMATE  LOCAL  ERROR  AT  ORDER  K 
AND  ERRORS  AT  ORDERS  K,  K-l ,  K-I  AS  IF  CONSTANT  STEP 
SIZE  WERE  USED. 


CHANGE  PHI  TO  PHI  STAR 

IF  (K  .LT.  NSP1 )  GO  TO  215 
DO  210  I  =  NSPl ,K 
TEMPI  «  BETA (I) 

DO  205  L  =  1 , NEQN 

205  PHI  (  L, I  )  =  TEMP1*PHI (L,I) 

210  CONTINUE 

PREDICT  SOLUTION  AND  DIFFERENCES 

215  DO  220  L  »  1 ,NEQN 

PHI (L,KP2)  =  PHI (L, KPl) 

PHI (L, KPl)  =0.0 
220  P (L)  =  0.0 

DO  230  J  =  1 ,K 
I  =  KPl  -  J 
IP1  =  1+1 
TEMP2  =  G (I) 

DO  225  L  =  1 ,NEQN 

P (L)  =  P (L)  +  TEMP2*PHI (L,I) 
225  PHI (L , I)  =  PHI (L,I)  +  PHI(L,IPl) 

230  CONTINUE 

IF (NORND)  GO  TO  240 
DO  235  L  =  1,NEQN 

TAU  =  H*P (L)  -  PHI (L, 15) 

P (L)  =  Y (L)  +  TAU 

235  PHI (L, 16)  =  (P (L)  -  Y (L) )  -  TAU 

GO  TO  250 

240  DO  245  L  =  1,NEQN 
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245  P (L)  -  Y (L)  +  H*P (L) 

250  XOLD  =  X 
X  =  X  +  H 
ABSH  =  ABS(H) 

CALL  F(X,P,YP) 

ESTIMATE  ERRORS  AT  ORDERS  K,  K-l,  K-2 

ERKM2  =  0.0 
ERKMl  =0.0 
*  ERK  =0.0 
DO  265  L  =  1,NEQN 
TEMP 3  =  1.0/WT(L) 

TEMP4  =  YP  (L)  -  PHI (L,  1) 

IF (KM2) 265,260,255 

255  ERKM2  =  ERKM2  +  ( (PHI (L ,KM1) +TEMP4 ) *TEMP3) **2 
260  ERKMl  =  ERKMl  +  (  (PHI (L , K) +TEMP4 ) *TEMP3) * *2 
265  ERK  =  ERK  +  (TEMP4 *TEMP3) **2 
IF (KM2) 280,275,270 

270  ERKM2  =  ABSH*SIG (KMl ) *GSTR (KM2 ) *SQRT (ERKM2 ) 
275  ERKMl  =  ABSH*SIG (K) *GSTR (KMl) *SQRT (ERKMl) 

280  TEMP5  =  ABSH*SQRT (ERK) 

ERR  =  TEMP5* (G (K)  -  G(KPl)) 

ERK  =  TEMP5*SIG(KPl) *GSTR(K) 

KNEW  =  K 

TEST  IF  ORDER  SHOULD  BE  LOWERED 
IF (KM2) 299,290,285 

285  IF (AMAX1 (ERKMl, ERKM2)  .LE.  ERK)  KNEW  =  KMl 
GO  TO  299 

290  IF (ERKMl  .LE.  0.5*ERK)  KNEW  =  KMl 
TEST  IF  STEP  SUCCESSFUL 
299  IF (ERR  .LE.  EPS)  GO  TO  400 


*  END  BLOCK  2  *** 

***  BEGIN  BLOCK  3  *** 

THE  STEP  IS  UNSUCCESSFUL.  RESTORE  X,  PHI  (*,*),  PSI(*) 
IF  THIRD  CONSECUTIVE  FAILURE,  SET  ORDER  TO  ONE.  IF 
STEP  FAILS  MORE  THAN  THREE  TIMES,  CONSIDER  AN  OPTIMAL 
STEP  SIZE.  DOUBLE  ERROR  TOLERANCE  AND  RETURN  IF 
ESTIMATED  STEP  SIZE  IS  TOO  SMALL  FOR  MACHINE  PRECISION 
*** 


RESTORE  X,  PHI (* , *)  AND  PSI(*) 
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305 

310 


315  * 


PHASEl  =  .FALSE. 

X  =  XOLD 
DO  310  I  =  1,K 

TEMPI  =  1.0 /BETA (I) 

IP1  =  1+1 

DO  305  L  =  1 ,NEQN 

PHI  (L,  I )  =  TEMPI* {PHI (L, I)  -  PHI(L,IP1)) 
CONTINUE 

IF  (K  .LT.  2)  GO  TO  320 
DO  315  I  =  2  ,K 

PSI  (1-1)  =  PSI  (I)  -  H 


C  ON  THIRD  FAILURE,  SET  ORDER  TO  ONE.  THEREAFTER,  USE 
C  OPTIMAL  STEP  SIZE 

C  IFAIL  =  IFAIL  +  1 

TEMP2  =0.5 

IF  (IFAIL  -  3)  335,330,325 

325  IF  (P5EPS  .LT.  0.25*ERK)  TEMP2  =  SQRT (P5EPS/ERK) 

330  KNEW  =  1 
335  H  =  TEMP2*H 
K  =  KNEW 

IF (ABS (H)  .GE.  FOURU*ABS (X) )  GO  TO  340 
CRASH  =  .TRUE. 

H  =  SIGN (FOURU*ABS (X) ,H) 

EPS  =  EPS  +  EPS 
RETURN 

340  GO  TO  100 

***  END  BLOCK  3  *** 


***  BEGIN  BLOCK  4  *** 

THE  STEP  IS  SUCCESSFUL.  CORRECT  THE  PREDICTED  SOLUTION, 
EVALUATE  THE  DERIVATIVES  USING  THE  CORRECTED  SOLUTION 
AND  UPDATE  THE  DIFFERENCES.  DETERMINE  BEST  AND  SIZE  FOR 
NEXT  STEP. 


400  KOLD  =  K 
HOLD  =  H 


CORRECT  AND  EVALUATE 
TEMPI  =  H*G (KP1) 


IF 

(NORND) 

GO  TO  410 

DO 

405  L 

=  1 , NEQN 

RHO  = 

TEMPI*  (YP(L)  - 

PHI (L, 1) )  - 

PHI  (L 

Y(L)  - 

P(L)  +  RHO 

405 

PHI (L,15) 

=  (Y  (L)  -  P  (L)  )  -  RHO 

GO 

TO  420 

410 

DO 

415  L 

=  1 ,NEQN 

415 

Y(L)  = 

P(L)  +  TEMPI* 

(YP  (L)  -  PHI 

(L, 1) ) 

420 

CALL  F (X, 

Y ,  YP) 
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UPDATE  DIFFERENCES  FOR  NEXT  STEP 

DO  425  L  =  1 ,NEQN 

PHI (L ,KPl)  =  YP (L)  -  PHI (L,  1) 

425  PHI  (L ,KP2)  =  PHI(L,KP1)  -  PHI(L,KP2) 

DO  435  I  =  1 ,K 

DO  430  L  =  1 ,NEQN 

430  PHI (L ,  I )  =  PHI (L,  I)  +  PHI (L,KP1) 

435  CONTINUE 


ESTIMATE  ERROR  AT  ORDER  K+l  UNLESS 
IN  FIRST  PHASE  WHEN  ALWAYS  RAISE  ORDER, 

ALWAYS  RECIDED  TO  LOWER  ORDER, 

STEP  SIZE  NOT  CONSTANT  SO  ESTIMATE  UNRELIABLE 

ERKPl  =  0.0 

IF (KNEW  .EQ.  KM1  .OR.  K.EQ.  12)  PHASEl  =  .FALSE. 

IF  (PHASED  GO  TO  450 
IF (KNEW  .EQ.  KM1)  GO  TO  455 
IF (KP1  .GT.  NS)  GO  TO  460 
DO  440  L  =  1 , NEQN 

440  ERKPl  =  ERKPl  =  (PHI (L ,KP2) /WT (L) ) **2 

ERKPl  =  ABSH*GSTR(KP1) *SQRT (ERKPl) 

USING  ESTIMATED  ERROR  AT  ORDER  K+l,  DETERMINE 
APPROPRIATE  ORDER  FOR  NEXT  STEP 

IF (K  .GT.  1)  GO  TO  445 

IF (ERKPl  .GE.  0 . 5*ERK)  GO  TO  460 

GO  TO  450 

445  IF (ERKM1  .LE.  AMIN1 (ERK, ERKPl ) )  GO  TO  455 

IF (ERKPl  .GE.  ERK  .OR.  K  .EQ.  12)  GO  TO  460 

HERE  ERKPl  .LT.  ERK  .LT.  AMAXl (ERKM1 ,ERKM2 )  ELSE  ORDER 
WOULD  HAVE  BEEN  LOWERED  IN  BLOCK  2.  THUS  ORDER  IS  TO 
BE  RAISED 

RAISE  ORDER 

450  K  =  KPl 

ERK  =  ERKPl 
GO  TO  460 

LOWER  ORDER 
455  K  =  KMl 
ERK  =  ERKMl 

WITH  NEW  ORDER  DETERMINE  APPROPRIATE  STEP  SIZE  FOR 
NEXT  STEP 
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460  HNEW  =  H  +  H 

IF  (PHASED  GO  TO  465 

IF  (P5EPS  .GE.  ERK*TW0 (K+l) )  GO  TO  465 
HNEW  =  H 

IF  (P5EPS  .GE.  ERK)  GO  TO  465 
TEMP 2  =  K+I 

R=  (P5EPS/ERK) ** (1. 0/TEMP2) 

HNEW  =  ABSH*AMAXl (0.5 , AMAXl  (0 . 9 ,  R)  ) 

HNEW  =  SIGN (AMAXl (HNEW,FOURU*ABS (X) ) ,H) 

465  .  H  =  HNEW 
RETURN 

***  END  BLOCK  4  *** 

END 

c********************************************** 

C  SUBROUTINE  INTRP 

C********************************************** 

SUBROUTINE  INTRP (X ,  Y , XOUT , YOUT , YPOUT , NEQN , KOLD , PHI , PSI ) 
DIMENSION  Y (NEQN)  , YOUT (NEQN)  , YPOUT  (NEQN)  ,PHI  (NEQN, 16)  , 
PSI  (12) 

DIMENSION  G (13) ,W (13) ,RHO (13) 

DATA  G (1 ) /l. 0 / ,RHO (1 ) /I . 0 / 

HI  =  XOUT  -  X 
KI  =  KOLD  +  1 
KIP1  =  KI  +  1 

INITIALIZE  W(*)  FOR  COMPUTING  G(*) 

DO  5  I  =  1,KI 
TEMPI  =  I 

5  W(I)  =  1 . 0 /TEMPI 

TERM  =0.0 

COMPUTE  G (*) 

DO  15  J  =  2 ,KI 
JM1  =  J  -  1 
PSI JM1  -  PSI ( JM1) 

GAMMA  =  (HI  +  TERM) /PSIJM1 
ETA  =  HI/PSIJM1 
LIMITl  =  KIP1  -  J 
DO  10  I  =  1, LIMITl 

10  W(I)  =  GAMMA *W (I)  -  ETA*W (1+1) 

G(J)  =  W(l) 

RHO(J)  =  GAMMA *RHO (JMl) 

15  TERM  =  PS I JMl 


INTERPOLATE 
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Fig.  G-2(e).  Control  Surface  Deflections  for  Pitch 
Rate  Command,  Flight  Condition  3--CSC  Mode  5 


O 
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Fig.  G-2(f).  Control  Surface  Deflections  for  Pitch 
Rate  Command,  Flight  Condition  3 — CSC  Mode  6 


Fig.  G-2(c).  Control  Surface  Deflections  for  Pitch 
Rate  Command,  Flight  Condition  3--CSC  Mode  3 


Fig.  G-2(d).  Control  Surface  Deflections  for  Pitch 
Rate  Command,  Flight  Condition  3 — CSC  Mode  4 
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Fig.  G-2(b).  Control  Surface  Deflections  for  Pitch 
Rate  Command,  Flight  Condition  3 — CSC  Mode  2 


Fig.  G-l(f).  Control  Surface  Deflections  for  Pitch 
Rate  Command,  Flight  Condition  2 — CSC  Mode  6 


Fig.  G-l  (d) .  Control  Surface  Deflections  for  Pitch 
Rate  Command,  Flight  Condition  2 — CSC  Mode  4 
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Pig.  G-l(a).  Control  Surface  Deflections  for  Pitch 
Rate  Command,  Flight  Condition  2 — CSC  Mode  1 


Fig.  G-l  (b) .  Control  Surface  Deflections  for  Pitch 
Rate  Command,  Flight  Condition  2 — CSC  Mode  2 


Appendix  G:  Control  Surface  Deflections 

In  this  appendix  the  control  surface  deflections 
due  to  q  and  p  commands  are  illustrated  for  flight  condi- 
tions  two,  three,  and  four. 
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P12 

DATA  VP12K/. 2418 ,-.2118, -6. 472, -5. 8 5/ 

DATA  VP12Z1/0. ,0. ,0. ,0/ 

DATA  VP12Z2/-. 06537, -.0066 97, -.01254, -.034 59/ 

DATA  VP12Z3/-. 2589, -1.861, -1.646, -.6861/ 

P22C 

•  DATA  VP22K/-4. 478, -17. 446, -51. 06,-14.168/ 

DATA  VP22Z1/0. ,0. ,0. ,0. / 

DATA  VP22Z2/ (-.205,. 853) , (-.219,1.607) , (-.3541,2.927) , 
+(-.4083,4.916) / 

DATA  VP22Z3/(-. 205  -  . 853)  ,  (- . 219 ,-l . 607 ) , 
+(-.3541,-2.927) , (-.4083,-4.916) / 

DATA  VP22P1/-.6835,-. 8265,-. 02719, -.03448/ 

DATA  VP22P2/- .1041, -.07795, -2. 6 97, -2. 171/ 

DATA  VP22P3/ (-.2741,1.909) , (-.211,1.953) , (- . 391 , 2. 962) , 
+(-.4996,3.129) / 

DATA  VP22P4/ (-.2741,-1. 909) , (-.211,-1.953) , 

+  (-.391,-2.962)  ,  (-.4996,-3.129) / 

P21 


DATA  VP21K/-4 .284,-13. 584,-51.72,-46.6/ 

DATA  VP21Z1/0. ,0.  ,0.  ,0.  / 

DATA  VP21Z2/(-. 3017, 1.562) ,  (-.2442,2.101) , 

+  (-.3749,3.578),  (-.3774,3.848)/ 

DATA  VP21Z3/(-. 3017,-1. 562) , (-.2442,-2.101) , 
+(-.3749,-3.578) , (-.3774,-3.848)/ 


U12  =  .25 
U21  =  .25 
Rl  =  10. 
R2  =  50. 
RETURN 
END 


..>1 
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Zdata 


SUBROUTINE  ZDATA 
■  DIMENSION  - 

DIMENSION  VPllK (4) ,VPllZl(4) ,VP11Z2(4) ,VP11Z3(4) 

DIMENSION  VP12K (4)  ,VP12Zl (4)  ,VP12Z2  (4)  ,VP12Z3 (4) 

COMPLEX  VP11P1 (4) ,VP11P2 (4) ,VP11P3 (4) ,VP11P4  (4) 

DIMENSION  VP21K(4)  ,VP21Zl  (4) 

COMPLEX  VP21Z2  (4)  ,VP21Z3(4) 

DIMENSION  VP22K  (4)  ,VP22Z1 (4) 

COMPLEX  VP22Z2(4) ,VP22Z3(4) 

DIMENSION  VP22P1 (4) ,VP22P2 (4) 

COMPLEX  VP22P3 (4) ,VP22P4 (4) 


-  COMMON  - 

COMMON  /R/  Rl ,R2 
COMMON  /U/  Ul2,U21 

COMMON  /VPllZ /VPllK, VP11Z1,VP11Z2,VP11Z3 
COMMON  /VP11P/VP11P1 ,VP11P2 ,VP11P3 , VP11P4 

COMMON  /VP12Z /VP12K,VP12Z1 ,VP12Z2 ,VP12Z3 

COMMON  /VP21Z /VP21K/VP21Z1 ,VP21Z2 ,VP21Z3 

COMMON  /VP22Z/VP22K,VP22Z1,VP22Z2,VP22Z3 
*  COMMON  /VP22P/VP22P1 ,VP22P2 ,VP22P3 ,VP22P4 


PLANT  DATA 
Pll 

DATA  VPllK/- 2 . 236,-5.862,-24.06,-32.9/ 

DATA  VP11Z1/0. ,0. ,0. ,0./ 

DATA  VP11Z2/-. 01822, -.01004 ,-.01263, -.02996/ 
DATA  VP11Z3/-. 4568, -.5502, -1.51, -1.097/ 

DATA  VP11P1/ (-1.3,  .0,),  (-2.028,0.)  ,  (-3.223,0.)  , 

+  (-.01516,  .02343)  / 

DATA  VP11P2/ (.  3633,0. )  ,  (1 . 167 , 0.),  (.  9645 , 0. )  , 

+ (-.01516,-. 02343) / 

DATA  VPllP3/(-. 07683, .2065) , (-. 006472 ,. 07803)  , 

+  (-.007553, .05384)  ,  (-.8012,6.592)/ 

DATA  VPllP4/(-. 07683, -.2065) , (-. 006472 ,-. 07803) , 
+  (-.007553,-. 05384) ,  (-.8012,-6.592)  / 
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IF (ISlA  .EQ.  0 )  D1A  = 
IF (ISlB  .EQ.  0)  DIB  = 
IF (IS2A  .EQ.  0)  D2A  = 
IF  (IS2B  .EQ.  0)  D2B  = 

SATURATION  - 

IF (D1A  .GT.  25.)  DlA 
IF (DlA  .LT.  -25.)  DlA 
17 (DIB  .GT.  25.)  DIB 
IF (DIB  .LT.  -25.)  DIB 

IF (D2A  .GT.  20.)  D2A 
IF (D2A  .LT.  20.)  D2A 
IF (D2B  .GT.  20.)  D2B 
IF (D2B  .LT.  -20)  D2B 


RETURN 


YP  (17)  =  R1  -  F1P*Y ( 17) 
FI  =  F1K*Y (17) 

YP (18)  a  R2  -  F2P*Y (18) 
F2  =  F2K*Y (18) 


G1 (S)  - 

El  =  FI  -  Cl 

YP  (19)  «  El 
G1A  =  Y (19) 

YP (20)  =  G1A  -  G1P1*Y (20) 

1  ■  YP (20)  +  G1Z1*Y (20 ) 

YP (21)  =  GIB  -  G1P2*Y (21) 

GlC  =  YP  (21)  +  G1Z2*Y (21) 

YP (22 )  =  Y (23) 

YP (23)  =  GlC  -  G1P3*Y (23)  -  G1P4*Y(22) 
GlD  =  YP  (23)  +  G1Z3*Y  (23)  +  G1Z4*Y(22) 
Gl  =  G1K*G1D 


G2  (S)  - 

E2  =  F2  -  C2 

YP  (24 )  =  E2 
G2A  =  Y  (24 ) 

YP  (25)  =  G2A  -  G2Pl*Y  (25) 

G2B  =  YP  (25)  +  G2Z1*Y (25) 

YP  (26)  =  G2B  -  G2P2*Y (26) 

G2C  =  YP  (26)  +  G2Z2*Y2) 

YP  (27)  =  Y  (28) 

YP  (28 )  =  G2C  -  G2P3*Y (28)  -  G2P4*Y<27) 
G2D  =  YP (28)  +  G2Z3*Y (28)  +  G2Z4*Y(27) 
G2  =  G2K*G2D 


DELTA  - 

G12  =  U12*G2 
G21  =  U21*Gl 


non  non  non  n  n  non  non  non 


Cll  =  P11K*D11B 

P12  (S)  - 

D12  =  D2A  +  D2B 
YP(5)  =  Y  (6) 

YP (6)  =  D12  -  P11P1*Y (6)  -  P11P2*Y  (5) 

D12A  =  Y(6)  +  P12Z1*Y (5) 

YP  (7)  =  Y  (8) 

YP (8)  =  D12A  -  P11P3*Y(8)  -  P11P4*Y(7) 
D12B  =  YP  (8)  +  P12Z2*Y(8)  +  P12Z3*Y(7) 

C12  =  P12K*D12B 

P21  (S)  - 

D21  =  D1A  -  DIB 
YP (9)  =  Y  (10) 

YP  (10 )  =  D21  -  P22Pl*Y  (10 )  -  P22P2*Y(9) 
D21A  =  Y  (10)  +  P21Z1*Y  (9) 

YP  (11 )  =  Y(12) 

YP  (12)  =  D21A  -  P22P3*Y  (12)  -  P22P4*Y(11) 

D21B  =»  YP  (12)  +  P21Z2*Y  (12)  +  P21Z3*Y(11) 

C21  =  P21K*D21B 

P22  (S)  - 

D22  =  D2A  -  D2B 

YP  (13 )  =  Y (14 ) 

YP  (14 )  =  D22  -  P22Pl*Y  (14)  -  P22P2*Y(13) 
D22A  =  Y  (14)  +  P22Zl*Y  (13) 

YP (15)  =  Y  (16) 

YP  (16)  =  D22A  -  P22P3*Y  (16 )  -  P22P*(5) 

D22B  =*  YP  (16)  +  P22Z2*Y  (16)  +  P22Z3*Y(15) 

C22  =  P22K*D22B 


C  - 

Cl  =  Cll  +  C12 
C2  =  C21  +  C22 


noon  n  n  on  n  no  oo  non 


Derv 


SUBROUTINE  DERV (T, Y , YP) 

RECONFIGURATION  SYSTEM 

DIMENSION  Y (28) ,YP (28) 

*************  COMMON  ************** 

COMMON  /P11Z/  P11K,P11Z1,P11Z2, P11Z3 
COMMON  /P11P/  PllPl,PllP2,PllP3,PllP4 

COMMON  /P12Z/  P12K,P12Z1,P12Z2,P12Z3 

COMMON  /P21Z /  P2lK,P21Zl,P21Z2,P21Z3 
COMMON  /P21P/  P21P1,P21P2,P21P3,P21P4 

COMMON  /P22Z /  P22K,P22Z1 ,P22Z2 ,P22Z3 
COMMON  /P22P/  P22P1 ,P22P2 ,P22P3 ,P22P4 

COMMON  /F/  F1K,F1P /F2K,F2P 

COMMON  /G1Z/G1K, G1Z1 , G1Z2 , G1Z 3 , G124 

COMMON  /GlP/GlPl ,G1P2 ,G1P3 ,G1P4  'j 

COMMON  /G2Z/G2K/G2Z1 , G2Z2 ,G2Z3 ,G2Z4 
COMMON  /G2P/G2P1 ,G2P2 ,G2P3 ,GP 

COMMON  /R/  R1,R2 
COMMON  /U/  U12,U21 
‘COMMON  /C/  Cl ,C2 
COMMON  /D/  DlA,DlB,D2A,D2B 
COMMON  /IS/  IS1A/IS1B , IS2A , IS2B 
COMMON  /M/  FI ,F2 ,G1 ,G2 , El , E2 , Dll , D12 ,D21 , D22 , 

Cll ,C12 ,C21 ,C22 


Pll  (S)  - 

Dll  =  DlA  +  DIB 
C 

YP (1)  -  Y (2) 

YP (2)  =  Dll  -  PllPl*Y  (2 )  -  P11P2*Y  (1) 
D11A  =  Y ( 2)  +  PllZl*Y  (1 ) 

C 

YP  (3)  =  Y{4) 

YP (4)  =  DllA  -  P11P3*Y  (4 )  -  P11P4*Y(3) 
DllB  *  YP  (4)  +  P11Z2*Y  (4)  +  P11Z3*Y(3) 
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DO  20  L  «  1 ,NEQN 
YPOUT(L)  =0.0 
20  YOUT(L)  =0.0 

DO  30  J  =  1,KI 
I  =  KIPl  -  J 
TEMP  2 

=  G  (I) 

TEMP 3  =  RHO(I) 

DO  25  L  =  1 ,NEQN 

YOUT(L)  =  YOUT(L)  +  TEMP2*PHI (L, I) 
YPOUT(L)  =  YPOUT(L)  +  TEMP3*PHI (L,I) 
30  CONTINUE 

DO  35  L  =  1 , NEQN 
35  YOUT(L)  =  Y  (L)  +  HI*YOUT(L) 

RETURN 

C 

END 
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Fig.  G-3(a).  Control  Surface  Deflections  for  Pitch 


Rate  Command,  Flight  Condition  4 — CSC  Mode  1 


O 

O 


Fig.  G-3(b).  Control  Surface  Deflections  for  Pitch 
Rate  Command,  Flight  Condition  4 — CSC  Mode  2 


Fig.  G-3(c).  Control  Surface  Deflections  for  Pitch 
Rate  Command,  Flight  Condition  4 — CSC  Mode  3 


O 

o 


Fig.  G-3(d).  Control  Surface  Deflections  for  Pitch 
Rate  Command,  Flight  Condition  4 — CSC  Mode  4 
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Fig.  G-3(e).  Control  Surface  Deflections  for  Pitch 
Rate  Command,  Flight  Condition  4 — CSC  Mode  5 


O 
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Fig.  G—  3(f).  Control  Surface  Deflections  for  Pitch 
Rate  Command,  Flight  Condition  4- -CSC  Mode  6 


o 
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Fig.  G-4 (b) .  Control  Surface  Deflections  for  Roll 
Rate  Command,  Flight  Condition  2 — CSC  Mode  2 
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Fig.  G-4 (d) .  Control  Surface  Deflections  for  Roll 
Rate  Command,  Flight  Condition  2 — CSC  Mode  4 


Fig.  G-4  (e) .  Control  Surface  Deflections  for  Roll 
Rate  Command,  Flight  Condition  2 — CSC  Mode  5 
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4 (f ) .  Control  Surface  Deflections  for  Roll 
Command,  Flight  Condition  2 — CSC  Mode  6 


Fig.  G-5  (a) .  Control  Surface  Deflections  for  Roll 
Rate  Command,  Flight  Condition  3 — CSC  Mode  1 
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Fig.  G-5  (b) .  Control  Surface  Deflections  for  Roll 
Rate  Command,  Flight  Condition  3 — CSC  Mode  2 
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Fig.  G-5(d).  Control  Surface  Deflections  for  Roll 
Rate  Command,  Flight  Condition  3 — CSC  Mode  4 


Fig.  G-5(e).  Control  Surface  Deflections  for  Roll 
Rate  Command,  Flight  Condition  3 — CSC  Mode  5 
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Fig.  G-5(f).  Control  Surface  Deflections  for  Roll 
Rate  Command,  Flight  Condition  3 — CSC  Mode  6 
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Fig.  G-6  (a) .  Control  Surface  Deflections  for  Roll 
Rate  Command,  Flight  Condition  4 — CSC  Mode  1 
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Fig.  G-6 (b) .  Control  Surface  Deflections  for  Roll 
Rate  Command,  Flight  Condition  4--CSC  Mode  2 


Fig.  G-6(c).  Control  Surface  Deflections  for  Roll 
Rate  Command,  Flight  Condition  4 — CSC  Mode  3 
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Fig.  G-6 (d) .  Control  Surface  Deflections  for  Roll 
Rate  Command,  Flight  Condition  4 — CSC  Mode  4 


Fig.  G-6  (e) .  Control  Surface  Deflections  for  Roll 
Rate  Command,  Flight  Condition  4 — CSC  Mode  5 
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Fig.  G-6  (f ) .  Control  Surface  Deflections  for  Roll 
Rate  Command,  Flight  Condition  4 — CSC  Mode  6 
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able  flight  control  system.  Compensators  are  synthesized 
to  control  pitch  rate  and  roll  rate  through  individually 
controlled  elevators  and  flaperons.  Robust  control  of 
these  variables  is  required  over  a  large  portion  of  the 
flight  envelope  despite  flight  control  surface  failures. 

Linearized  aerodynamic  data  are  used  to  develop 
the  aircraft  model  in  state-variable  format.  The  longi¬ 
tudinal  and  lateral-directional  equations  are  coupled  in 
the  control  matrix.  Individual  control  of  the  elevators 
and  flaperons  is  obtained  by  dividing  the  dimensionalized 
control  derivatives  for  a  control  surface  pair  in  half  and 
assigning  each  surface  of  the  pair  one-half  of  the  total 
derivative  value.  The  system  with  individually  controlled 
surfaces  represents  a  four  input-two  output  system  which  is 
transformed  into  an  equivalent  two  input-two  output  system 
for  each  control  surface  configuration  and  flight  condition. 
Quantitative  feedback  theory  is  then  applied  to  the  equiva¬ 
lent  systems.  ~ 

Fixed-compensators  and  prefilters  are  synthesized 
for  the  pitch  rate  and  roll  rate  channels,  and  a  digital 
simulation  is  conducted  at  four  points  in  the  flight 
envelope  with  up  to  two  simultaneous  surface  failures,  as 
well  as,  for  the  unimpaired  aircraft.  Results  of  the  simu¬ 
lation  show  that  the  compensators  provide  robust  control 
for  each  control  surface  configuration  at  each  flight  condi¬ 
tion.  Compensator  designs  are  also  synthesized  to  provide 
only  stability  when  three  of  the  four  surfaces  are  failed. 

The  theory  developed  by  Dr.  Horowitz  is  found  to 
be  effective  in  designing  for  robustness  despite  large 
uncertainty  arising  from  control  surface  failures,  and 
plant  parameter  variation.  Flight  control  research  using 
quantitative  feedback  theory  should  be  continued.  The 
fixed-compensator  design  approach  has  potential  to  signifi¬ 
cantly  reduce  the  complexity  of  a  reconf igurable  flight 
control  system.  A  computer-aided  design  package  for  quanti¬ 
tative  feedback  theory  should  be  developed  to  facilitate 
design  work. 
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